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PNSASTS 


Fortran (‘ormula TRANslation ) is a computer language that has been around 
for a@ quarter of a century: publication of this book coincides with fhe 
twenty-fifth anniversary of the first Fortran compiler. Fortran is still much 
used despite predictions throughout its life that more elegant languages 
would supersede it. But Fortran is still the only language in which it is 
possible = with care = to write really portable programs. A portable pro- 
gram is one that should work with little or no alteration on any computer 
that has Fortran: in other words on most computers of any power. Some 
firms employing less than twenty people have invested more than fifteen 
man-years in Fortran programs: in big companies such investment can be 
enormous. The potential portability of Fortran has enabled such firms to 
transfer their investment from an outdated computer to a modern one «= 
in some cases several times during the firm’s history. Changing to a 
more elegant programming language would be more expensive than 
sticking to Fortran # so Fortran is here to stay for a while. 


Phe Fortran described and illustrated in this book is that defined by the 
American National Standards Institute in 1966 @ Fortran 6 ; roughly fortran 
IV) but with allusions to the same Institute’s standard published in 1978 
and affectionately known as fortran 77. At present not every computer can 
understand the new language, but even if it were possible to write widely 
portable programs in full Fortran 77 it would not be practical fo present 
the whole language in a book such as this. Fortran 77 = of which 
Fortran 66 is a subset = is a QYA language. 


b| KX > Kx 


THICKNESS OF ANSI THICKNESS OF ANSI 
REPORT: FORTRAN 66 REPORT: FORTRAN 77 


“Phe most important aim of this book is to provide a reference manual 
for Fortran emphasizing the self discipline needed to achieve portable 
programs. In deciding what could safely be advocated for portability I 
referred to Compatible Fortran and the HECB Programming Instruction Manval 
= both cited in the bibliography. Apart from using swréceR and REAL declar- 
ations € forbidden by the wecs manual 2 I believe my presentation of 
Fortran 66 takes account of all other important recommendations made in 
those two books. I have described every facility in Fortran 66 but 
encouraged the reader not to use €or use only with caution ) those that 
might lead fo non-portable programs. 


SPhe second aim of this book is to introduce programming fo the person 
who has not tried to program a computer before, or has used only the 
ubiquitous language called sAs/c. No previous knowledge of computers or 
computing is assumed, but the novice may not be able to take the advice 
given to Alice: “Begin at the beginning ... and go on till you reach the 
end: then stop.” That is because this book, being in the form ofa 
reference manual, has forward references here and there. So the novice 


Vtit 


is invited to skip certain pages @ these are few and clearly marked ) until 
he or she has a rough grasp of programming in Fortran before returning to 
cope with the subtleties. 


phe third aim of this book is to illustrate a few tricks of the pro- 
grammer’s trade such as sorting numbers, plotting graphs, solving simul- 
taneous equations, finding routes through networks and decoding 
algebraic expressions. Several of the examples show how fo handle 
characters € letters, digits and symbols ) in programs that should neverthe- 
less be portable. 


‘he examples of Fortran programs in this book are written with the idea 
of clarity in mind; not efficiency. A program is seldom a static thing; 
more often it remains in a state of periodic Improvement, extension and 
correction, An intelligible program is easier to maintain than one that 
Sacrifices clarity for efficiency. Because programmers’ time gets ever 
More expensive = and electronic computation cheaper and cheaper «= the 
clearly written program makes more sense economically than the “clever” 
and efficient one. Even so the reader is sure to find pieces of program 
in this book that cry out for greater efficiency with little loss of clarity. 
By all means note any desired alterations in the margins of this book 
Wherever my programming offends your sense of style or thirst for 
efficiency. 


An Claiming the advantages of clarity some may protest there are few 
comment lines in my examples. That is because most comments are made 
in adjacent explanatory texts or in little clouds Clike“this so as to fit 
a complete example on one page. 


Reigate, Donald Alcock 
Surrey, U-K. May 1982 
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WT HE: CORCE PT ees are M” WITHOUT INVOLVING 


6-5 
A\ssume there is no computer to 

help solve Se problem connate Pine 

Q painter: how many pots o 

paint does he need to paint the kth | 
roof and wall of this water tank ? pee TAN 


‘Phe paint manufacturer says each pot has 
enough paint to cover an area of 236. 





femembering that the area of a circle is given by mr* ( where r is 
its radius) or %d?+4 (where d is its diameter) the painter can work 
out the area of the top of the tank: 


AREA OF TOP = 3.14 x 6.57 +4 = 33.1663 


[Remembering that the circumference of a circle is given by td 
(where d is its diameter as before) the pointer can work out the 
area of the wall of the tank: 


AREA OF WALL = CIRCUMFERENCE X HEIGHT 
= 3-14 x 6.5 x 27-0 = 551.070 


‘Whe area to be painted is the sum of the areas of top and wall: 


TOTAL AREA = 33.1663 + 551.070 = 584.2363 


Qnto this area must be divided the coverage of each pot of paint so 
as to give the number of pots needed: 


POTS = 584.2363 + 236.0 = 2.47558 


you cannot buy a fraction of a pot of paint, so the final answer 
must be rounded up to the next whole number. (Looking at it 
another way, take the /nfegra/ part of the value and add unity. ) 


NUMBER OF POTS = QE +1 = 3 
Ws 


dr the items of data had been such that Pots had worked out at, 
say, 3.00000 instead of 2.47658 then the solution would have been 
3+1= 4 pots. Although this would be wrong mathematically a 
prudent painter would be happier with the extra pot in case he 
spilled a few drops of paint. 


D/AM 
IMOW; suppose the painter wanted to set 
down this method of calculation so that he 
could more quickly calculate the number of 
pots of paint (having perhaps oa different | 


capacity) to paint a water fank of any 
diameter and any height. p= HEIGHT 


rAN possible list of instructions | 
is set out on the opposite page. oogt 








/[\ Draw three little boxes to contain the items of data needed to solve 
° the problem. Name these boxes DIAM, WEIGHT, COVRG: 
pbam{ sd meiGHT[ «id coveg 

also drow boxes to contain intermediate results of the calculation: 


eee 


Draw a little box to contain the number of pots needed; but draw 
° it differently to emphasize that its contents must be a whole 


number = an énteger. NPOTS 


Ci Read three items of data for ao particular case to be solved, and 
© put them in the boxes named DIAM, HEIGHT, COVRG. 


ms EXAMPLE <> 
DIAM HEIGHT covrc[ 236.0 | 


Work out the area of the top of the tank by the formula 3.14xd*z4 
where d is the diameter found in the box named pDIAm. Put the 
answer in the box named /7oP. 


rop (“35-1663] 


York out the orea of the vertical wail of the tank by the 
formula 3.144xdxh where d andh are the diameter and height 
found in the boxes named DIAM and HEIGHT respectively. Put the 
answer in the box named WALL. 
























+ 3.14 x 6.5% + 4-0 =e EXAMPLE = 





7 ARAMIEE, WALL|_ 551.070 |] + 3.14 x 6-6 * 27.0 









ZA\ad the areas found in boxes named 7oP and WALL. Divide this 
sum by the coverage of oa pot of paint. This is found in the box 
named Covrcg. Pur the answer in the box named /Po7s. 










EXAMPLE © (33, 1663 + 551.070) + 236.0 -> POTS 
Dow take the integral port of the number found in the box 


named ors; addi ; put the result in the integer box named 


NPOTS . 
i enmeate 2Ga7sss) +1 > wpors [3 T 


[rl Write out the solution; in other words wrile out fhe integer found 
© in the box named wpors. Give this note to the painter. 












= EXAMPLE 










Sror work: there is nothing more to do. 


This is the END of the list of instructions. 








Qn computer jargon such a list of instructions @ ex- 
cluding the embedded examples 2 is called a 
program; the person who writes programs is called a 
programmer. The obove program is written in 
English, but over the page if is translated into 
fortran which is a language a computer can 
understand and obey. The jargon for obey is 
execute.) The instructions are obeyed in sequence 
starting at the first and stopping on meeting 57opP. 





3 


ALG STP THON BE Ah iy la 


REAL DIAM, WEIGHT, COVRG, TOP, WALL, POTS AR 
INTEGER . NPOTS = 


DIAM, HEIGHT, COVRG Z, 


RITE 6,200) NPOTS 
@RMAT (1X GH YOU, NEED, I2 5H PoeTS, ) 


G=> 8S Gun9 ® BE 





lAlere is the program in Fortran. The letters A to /Y with their little arrows 
(not part of the Fortran) correspond to steps in the English program on 
the previous page. There ore many niceties embodied in the Fortran 
notation above so it worth working right through the program yet 
again =» pretending to be the computer. 


[Xeep the items of data the same as before so as to simplify com- 
parison with the program written in English. These items are written 
on a special form rather like the one used for the Fortran program 
above. Character positions are called columns for reasons ex- 
plained later. 





Start obeying the Fortron program at its first /sastruction € also called 
a statement ) and continue in sequence until you meet one which 
says STOP, 


To REAL DIAM», HEIGHT» COVRG», TOPs WALL», POTS} A 

Lhis tells the computer to name six “little boxes” with the names shown. 
For the English program we drew the boxes. The computer, however, uses 
storage locations each capable of storing a R&AL number (in other words 


@ number with a fractional part after a decimal point 2. The locations 
DIAM, HEIGHT, COVRG, TOP, WALL, PoTS are called REAL variables. 


i eeEN JORIS 


his fells the computer to name another “litle box” WPOTS. In the 
English program we drew this box differently from the boxes designed to 
hold REAL numbers. This box is designed to hold an /N7EGER. The location 
MPOTS is called an INTEGER variable. 


 ... | READ (5,100 DIAM, HEIGHT», COVRG 
0. 
| 1,00] |F ORMAT, 3F 10.0) 


“Jhis instructs the computer to read some waiting items of data and put 
them into variables D/AM, HEIGHT, CovRa. The FORMAT statement says there 
should be 3 waiting items, each in fixed point form (hence the 3F ) 
and that they reside in successive fields of 10 columns each @ hence the 

10.0). Chapter 10 deals with formats in detail; until then all examples 

use F10.0 for reading real numbers which shovid be written with decimal 
points. 




















WS @ 





Whe 100 (in READ(5; 100) ) associates the READ statement with its FORMAT 
statement, This number is arbitrary and the fwo statements do not have to 
follow one another. Indeed it is common practice to group FORMATs together. 


4 


Phe 5 (in READ(5,100) }) is a unit number. You can read from all sorts of peri- 
pheral devices each of which is associated with an unique unit number. Con- 
ventionally unit 5 denotes a card reader. All examples until Chapter {1° have 
5 as the unif number in READ statements. 

<x 


TOP = 3.14 * (DIAM#*2) / 4.0 
rg egal el 


Whis is as described for step D of the English program, but notice the use 
of * for multiply by, the use of **# for raise to a power, the use of / for 
divide by. The brackets are not necessary here (exponentiation is always 

done first }) but they clarify the instruction and do no harm. 








[For the equals sign don’t say “equals” @ this is not an equation in algebra ) 
say “becomes”. Thus: Top becomes 3.14 times DIAM to the power 2 divided by 4.0. 
Statements with an equals sign ore called assignment statements or assignments. 


ALL = 3.14 * DIAM .* HEIGHT <1 [z 
—— 4 eee — el 


This assignment says: WALL becomes 5.14 times DIAM times HEIGHT. It corresponds 
precisely to step E of the English program. 


TR OTE STOP A WALL DZ COUR nF 


Vhis assignment says: /073 becomes (TOP plus WALL) all clivided by COvRG. Notice 
the importance of the brackets. Without them Pors would become 70P plus 
WALL/COvRG which is ridiculous. In the absence of brackets, multiplications 
and divisions are done before additions and subtractions » just as you 
would expect from the rules of algebra. 

1g 


| WPOTS = INTCPOTS) + 1 


Whis assignment has oa function on the right-hand side. This function, 
INT( ), causes the computer to consult whatever real number is represented 
inside the brackets and take just its integral part. Thus the statement says: 
NPOTS becomes the tntegral part of PoTs plus one. (Fortran has many other 
functions as described in Chapter 6 .) 

peti 


“A 


OE 
| 200 FORMAT <1X GH YOU NEED, 32, 5H POTS 


Uhis instructs the computer to write out the value found in the integer variable 
NPOTS. The 1X in the FORMAT statement says to start a new line. The 9H says 
to print the following 9 characters “ you NEED” which includes the space before 
yOoU and the space before NEED. The I2 says fo print the integer @ found in 
variable NPOTS ) in a field of two columns, justified to the right. The 5H 
says to print the following 5 characters “ Fors” including the initial space. 
So the result should look like this: % YOU NEED 3 POTS 5 

















VW hy does ¥ signify characters ? This is history. Hermann Hollerith pioneered 
punched-card equipment in the 1890s and his initial, H, survives in the termin- 
dlogy of Fortran. Characters may be called Hollerith characters. 


‘Phe 200 Cin WRITE(6,200)y associates the wriTE with the FORMAT statement as 
explained for the READ statement above. Again this number is arbitrary and 
the FoRMAT statement need not follow the wei7e. The 6 denotes a unit number 
as described for the READ statement. Conventionally unité6 is associated with a 
line printer and is so used throughout the examples in this book. 

ic ee he et ee 
Whis tells the computer to stop obeying instructions « in other words to stop 
execution. 


Ne erg ES cats te gene ps oe Da no ee 


his simply marks the end of a deck of cards for a program like a lid on a 
box. It is ot an instruction: it is called the END line. 











eSS PAULIN iro'A corm A COMPUTER CAN READ 


Pplere is @ program and its data written on standard coding forms: 


C TOP + WALL / COVRG 
POTS = INTCPOTS + 1 


RITE 6,200) NPOTS 
FORMAT C1X%9, 9 YOu E,ED Z,2 SH. POTS. 
BO gs ee Pk aoa Ly. Hal ae ee 1 


FORTRAN DATA program: 
1 0 20) 





"Phe traditional way of encoding such material is to punch cards ata 
device called a keypunchs BX 





a HOPPER 
GPhis leaves you with two 
decks of cards: oa 

program deck [A 


PUNCHING 
STATION 






KEYBOARD 
LIKE A 
TYPEWRITER 






and a data deck. >>. 


ZAnother traditional way of encoding programs and data is to punch 
paper tapes ot a device called a teletypewriter, 
Tt iy re “This leaves you with two 
reels Of punched paper 
tope: QO program tape, 
a) ae 
wl; 
e% 

and a data si gia 


Pe 
ee xf 
. 


TYPEWRITER 
KEYBOARD 





EMERGING 
PUNCHED 
TAPE 


Nowadays the most common method of getting a program and data 
into a computer is to type the material at a visual display unit 


@ vbU for short >. 


DISPLAY 
SCREEN 


TYPEWRITER 





‘fhe information from this source is stored magnetically on the surface 
of a disk, The program constitutes a Fortran file and the data 
constitutes a data file. The disk may be a “hard” disk: 


K#— 1/4 INCHES >| 
i. ee <KJ EXCHANGEABLE 
———= DISK PACK 
————— 
lS (6 DISKS ON 
ONE SPINDLE ) 





pur DING WINDOW 


CUT IN ENVELOPE 


or a “floppy” disk: 
‘ 


a 
FLOPPY D/sk Ee 4 


or the information may be stored as files on a reel of magnetic tape 
or on a cassette tape: 





ZA REEL ABOUT THE 


SIZE OF A 
p\ DINNER PLATE Lo 
ee an 


SKI (NGNENC-TAPE 
DRIVE 





VZhatever the medium there are ways of correcting mistakes in typing; ways 
of deleting lines, inserting lines, changing individual characters within lines, 
and so on, This is called editing. And you end up with two sets 
Of information the computer can read: 

@ a program 

@ a set of data 


EXEC TOM Wie Reva yy aa 


Using the “traditional” forms of input for illustration, what appears to 
happen when you feed your program and data to the computer is this: 







You put the program tape in 
the tape reader @ or program deck 
into the card reader) and the 
computer reads your program 
into its memory. 


RESULTS [>> VYOO WD BPS 3 


7 
PROGRAM = TAPE 


TAPE: D) 







Then you put the dato tape into the 
tape reader @ or data deck into the 
card reader } and the computer starts 
to execute your program = faking its oe 
data from the reader and printing 

results on the line printer. 


IN MEMORY 


fut it doesn’t happen as simply as that. The computer can certainly 
understand and execute instructions, but only a program in binary form 
s something like this: 

0010111001011010 

Oo1010OLo1oOO1o1L 

L1O100L0110L 01101 

L1110010110000110 

etc. etc. 


Qne such program is called a fortran compiler. This program is able to 
read your Fortran program and translate it into Os and 1s = then hand 
control over to the first instruction of your newly-translated program. In 

computer jargon a compiler comptles a program from source form € for 
example from Fortran form ) to object form €0s and ts). So the simple 
two-stage process depicted above is more nearly as follows: 





COMPILATION 
IN SOURCE FORM CFROM TAPE TO Disk) 
(STOPS AT THE “END” LINE) 


YOUR 
PROGRAM 
NOW IN 
OBJECT 
FORM 


YOUR RESULTS 
YOU MEED 3 POTS 







INTO MEMORY 


civ) 


EXECUTION 


TELLING THE COMPUTER WHAT 
OPERA TANG SYSTEMS: RESOURCES TO USE AND WHEN 
Whe computer has fo be told in detail what compiler to use and what 
peripheral devices @ card reader, tape reader, magnetic tape drives, disk 
files, card punch, tape punch, line printer ) to use. This cannot be done 


just by throwing switches or pressing buttons. You have to instruct an 
operating system. 


An operating system is a clever program rather like a compiler. Its 
job is to make the connections referred fo opposite « for example 
ensuring the Fortran compiler on disk is read into the computer’s 
memory, or that results from a particular program are routed to the 
line printer. Just as you have to know the language called fortran to 
write the kind of instructions illustrated on page 6, you also have to 
know the language of the operating system to tell the computer what 
peripheral devices to use and what compiler to employ. 


Whe language of an operating system is called a job-control language 
q@/cL for short). Every computer has one. Even a modest microcomputer 
has a JCL ~ with a vocabulary perhaps no more extensive than the words 
LOAD and SAVE and RUN. But a large modern computer serves many users 
simultaneously = each wanting to employ a different compiler @ Fortran, 
COBOL, Pascal, APL éfc.) or be competing for the next turn oat one of the 

line printers. Such an installation has a daunting JCL with huge vocab- 
vlary and complicated grammar. You may have to “talk” fo such an 

operating system over a telephone line by typing commands in JCL at 

a visual display unif. 


Below is illustrated a “fraditional” program deck and data deck = 
interspersed with commands fo an operating system written in a typical 
( but imaginary p JCL. 






PROGRAM 
DECK 





4J0B 136 USER TFL 
io 1 puso oo 






SPhis book does not describe operating systems because there are so many 
of them and they are so different in nature from one make of computer 
to onother. There are few short cuts to learning the JCL of an 

operating system. You should persistently ask questions of those who seem 
fo know it all « study whotever manual is available even if it seems, at 
first, incomprehensible « and above all try out your efforts fearlessly. If 
the operating system then mokes a mess of everything _in the_ computer, 
losing everyone’s files, if is 3 









CHAPTER 1 


fod wry making the program about water tanks run on the computer 
to which you have access. Try it with several different sets of 
data. q This exercise might take longer than you expect. ) 


1.2 Hote the necessary job-control cards @ or what you have to type 


when “signing on” and submitting such a job from your terminal ) 
far future reference: 


Kt proceditre FOF thE veces snvisssvicewnecn vnoee, CRStAllatION 
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AOS TORE 


PUNCHED CARDS 
LINES 

LABELS 

PROGRAM UNITS * 
ORDER * 
EXERC/SES 


% SK/P ON FIRST READING 


HERMANN HOLLERITH PROVIDING SOME 
OF THE TERMINOLOGY OF FORTRAN 


You may not have to use punched cards. Nowadays if is common to 
use fortran from a terminal to a large computer = or on a personal 
computer where the means of communication is oa typewriter keyboard. But 
Fortran was originally designed with punched cards in mind « and some 
ideas and terminology derived from punched cards remain in the 
terminology of Fortran. 


So here is a punched card: 
Ty FoRMAT cariu. 09 


i FORT RAN STATEMEN T IDENTIFICATION 
FD OODOMOOOODOMOMOOOODODDODDDODOODDDDUDODDDODI@M DIN 0D000NN 


lasal PVM U NIMH EAB VIM MRE GUMS we Oe a ob ot es ee 
| 


rT eet RRREG! (POPP e eee eee eee eee enna eee eeeee Prajbadtyadd 
pai Be i dd) lB Si ay Ul ae By Sy Sy ay i So By Sy ee Sa A a ei SS oe De 22/22222222 


| 

t 

AS ATS UASAY ESA EA ON EOS SESSA BO 33/33333333 

NeacaldesaQuascrcaaaacaa cect cece aaleaaaegaa 
| 
l 
| 
| 


SSSSSSISSSSSSSSMSSSSSSMPSSSSSSSSSSSSSSSSSSSSS5S5S5555 SIS 555555555 
SHG G GG MMC CCFC CCG MGCCC SSC EGGEEGEGGGG EGCG EGG GGGE6EEEE PpbiGG GE EGES 
WTA TT TTT TTT ITTV TTT TTT T TTI TTTT TTT TT TITIAN 


| 
Le ne $9.95°5'9'3 $ 2 $9595S595 
aie Hl” gas tig . Parr es 45:8; 


prize sts ae eee fim) PUUD Bee ah be Sf 


ee a Be | 





POL Otivt 





[Zach charocter occupies one column of which there are eighty per 
card. 


“AA card is punched by typing on a keyboard much like that of a type- 
writer where o fouch on a key produces a corresponding pattern of 

holes in the column of a card. <A touch on the space bor produces a 
blank column. A fouch on a special key causes the current card to be 


ejected into the ovtput hopper and o new card made ready for punching 
in column 1. 


An Fortran statements the columns are employed as follows: 
@ blank cards are not allowed in the body of the program 


@ columns 1 to 5 are for the label if any @note labels 100 
and 200 in the introductory example p 


@ column 6 is normally left blank: if not blank or zero this 
signifies continuation of a statement from the previous card 


@ columns 7 to 72 are for the body of the Fortran statement 
columns 73 to 80 are ignored by the Fortran compiler 
@ column 1 is special. The letter C in columni indicates the 


whole card is a comment to be ignored by the Fortran 
compiler during compilation. 






Qt is usual to write Fortran on coding forms divided into eighty 
columns and with columns 1,6 and 72 emphasized in some way: 


DO NOT LEAVE 
SF] BLANK LINES 





Qn the examples in this book column 6 is indicated by two vertical 
lines. Columns elsewhere in Fortran statements are usually not marked. 
Indication is unnecessary because blank columns ore generally ignored 
inside Fortran statements in columns 7 to 72. 

‘here is Oo notable exception to this. In Hollerith items blank columns 
ore not sila after letter A: 





i et Oe ee ee ee ee 


[Sater AT _C1X>, on” 1 VOY WEED NEED, 1.21, 5H POTS ) 
hoo 
Ca 










Whe rules for punching @€ or typing ) data for a Fortran program are 
completely different from those set out above. So if is customary fo 
employ a differently printed card. 


DATA 


UEVEREEER CORRS On ee Cer nir si beer a. ie ae rerryer ery ian Ar TMS se 


$000000000|000000§ 000000000 00cl09005cc00N0000000 pi ieoagons 
VUVVELVEL VTP TATED EAT ATED Tuddaerairat 


statements. 
numbers 
Qa common 
on the punched card above, 
Chopter 10. 


22222222292228222222)2 28 2222222122222222222222222 
93320333333333383333/3 338 38333 3)339393333993133333333 
VUCPRP EUR UCCURUEES  CCCUPEEESY CEEECECCE ECE EE) 
$5 5B5555/55555555555555555555559595955 555555555 
COSMSSEEESECEEGECGESGGE GB S666 666666666 66666565! 
PPRERRRRRIGEE) SRERE PROPER Pee ePe eer iPeeeee: 
CECH SS esse SSSM SS RH RES SM Ass ss OR CSAS escsegease 


‘ : 


= LNLEO BIMGOOM  nsbtilhw 


ANG G00 809/590 99 9995585090989; 99989 a8 
pris nae 38) 


convention: note 


THESE peer 








2212222222222 
3413333133323 
444aehgggqas 
$5)/5555555559 
BEERES SEG SHS 
THIVITIIIIIG 
ees eseseses 





Data are punched in fields as dictated by their respective FORMAT 
In the introductory example FfFORMAT(3F 10.0) demanded three 
located in three successive fields of fen columns each. This is 
the positions of the vertical 
Formats are dealt with exhaustively in 


COUNT AS ZEROS. 
°T LEAVE CAVE BIANKS INSIDE I 


DEER SEI 


lines printed 





THE CONTENT OF A_ SINGLE PUNCHED CARD 
= OR /TS EQUIVALENT IN SOME OTHER MEDIUM 


Zn End line is not \ike any other statement in Fortran. It simply 
marks the end of a program or subprogram Q(we meet subprograms 
over the page). The letters £& then W then D ore punched anywhere 
in columns 7 to 72 leaving the rest of the card blank: 





Some Fortrans (including Fortran 77) permit the £wd line to imply S7oP 
in the main program and RETURN in a subprogram, but this is not 
allowed in Fortran 66. 


“AA continuation line has a character other than blank or zero in column 
6. There may be up to nineteen continuation lines following the initial 
line of a statement. The standard does not forbid it, but some 
Fortrans object to labels on continuation fines. Such labels can serve 
no useful purpose anyway. 


WALL = 3.14 


x DIAM 
“E CONTINUATI Ni 
& WEIGNT gh CONTIRUATION Lines) 


The use of 1,2 etc. as the non-blank character in column 6_ illustrates 
a tidy convention for continuation lines. 





Dlever split a Hollerith item @ such as 9H, YoU.NEED) between one line and 
a continuation line. This would cause a Hollerith item to be filled out 
with blanks, thus making the next continuation line begin with rubbish. 
Although Fortran does not specifically forbid it, do not split any of 
Fortran’s keywords @ such as FORMAT) across continuation fines either; 
some Fortrans object to the practice. 

S34 











MAT (1X, 9H YOU 
NEED, 14, SH PD > 


‘Whe comment line has letter C in columni1, fhen any comment € com- 
posed of characters from the standard set = see page 20) in subsequent 
columns. It is safer to leave column 2 blank because there are 
Fortrans that attach some significance to characters here. Comment lines 
may be inserted anywhere before the £nwo line except immediately before 
continuation lines. 


REAL DIAM, HEIGHT, COVRG, TOP 
IS BEFORE A CONTINUATION LINE 
WALL, POTS 
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FOR ANY EXECUTABLE STATEMENT AND ALL 
FORMAT STATEMENTS 


“A label may be written as digits anywhere in the first five columns: 





Although standard Fortran permits, it is safer not to spread out 
the digits of oa label or write a label with leading zeros: 





Standard Fortran provides no limit to the size of integer used asa 
label = and different Fortrans impose different limits. A common limit 
(suggested as the one to ensure ao portable program) is 32767 (2-1). 
Every label must be unique in its program or subprogram. 


Statements in Fortran are either executable @i.e. cause something to 
happen when obeyed } or non-executable (i.e. declaratory in nature 
and not subject to being “obeyed”. Referring to the introductory 
example, statements are marked £4 and wWoNEXx respectively in the 
reproduction below: 


REAL DIAM, HEIGHT, COVRG, TOP, WALL, POTS << HONEX) 
INTEGER NPOTS SLWONEX) 


READ (5,100) DIAM, HEIGHT, COVRG 

FORMAT (3F10-0) 
TOP = 3.14 * (DIAM ##2) / 4-0 

WALL = 3.14 ® DIAM * HEIGHT 


POTS = (TOP + WALL) / CPVRG 
NPOTS = INT(POTS) + 1 
WRITE (6,200) NPOTS 


FORMAT (1X, QH.YOUNEED, 12, 5H.POTS ) 


STOP 


END 





“\ny executable statement may be given a label. The only non- 
executable statement that may be labelled is the FORMAT statement = and 
that statement must be labelled to be of any use. € Standard Fortran 
does not specifically forbid labels on non-executable statements but 
these have been known to cause trouble in some otherwise standard 
Fortrans .) 


ovices fo programming should skip the rest of this chapter on first 
reading. 
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PROL:RAA THE MAIN PROGRAM & ITS SUBPROGRAMS 
: (SKIP THIS PAGE ON FIRST READING ) 


A Fortran program may be just a simple main program as illustrated in 
the introductory example. On the other hand the main program may 
refer to @or invoke ) any number of subprograms any number of 
times. 


pects 
Pee marato (4.8.0) = 













ANOTHER 
SUBROUTINE 
SUBPROGRAM 






FUNCTION 
SUBPROGRAM 
UNCTION AREA(P, Q,N) 





| ae main program may 
share COMMON storage INITIALIZE CERTAIN DATA 


@ see page 76) with one 
or more subprograms. Initial values may be set Tee Cis 
in such storage by a BLOCK DATA Subprogram . 


Subprograms invoked by the main program may, in turn, invoke 
other subprograms. The only restriction is that no main program or 
subprogram may invoke étse(f — either directly or indirectly. In other 
words recursion* is not allowed. 


a. SUBROUTINE A 
CALL B 


# it is possible fo use the 


Where ore three kinds of subprogram: principle of recursion in 
Fortran by eaitiog. bs up your 


@ SUBROUTINE @ Chapter 7 ) own stacks efc.. 
@ FUNCTION @ Chapter 7 )p eke rophy: Pin ae 
@ B10CK DATA @ Chapter 9 ) 4g 





very subprogram begins with a heading containing one of the above 
names and ends with an END line as illustrated. The main program 
(of which there must always be precisely one ) is distinguished by the 
absence of any special heading. Many Fortrans demand such a 
heading beginning with the word PROGRAM or MAIN but this should be 
considered part of the job-control language rather than a Fortran 
statement. 


‘Whe term program unit) is used to mean a main program or a 


subprogram. 
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OR DE:R OF STATEMENTS IN A PROGRAM UNIT 
(eK: @ SKIP THIS PAGE ON FIRST READING ) 


‘Whe order of statements in a program unit is defined by the following 
table: 


One SUBROUTINE or FUNCTION or BLOCK DATA statement 
Cor nothing at all if this is the main program )p 


REAL, INTEGER, DOUBLE PRECIS/ON, LOGICAL, COMPLEX 
@é.e.type statements)in any order among them- 
selves 


DIMENSION 
COMMENT statements EXTERNAL 


lines statements 

anywhere | COMMON Cintermingled among 
between statements COMMON and DIMENSION 
heading statements p 


and £ND 
line but EQUIVALENCE statements 


not 
between DATA statements @at least one in a BLOCK DATA sub- 
continuation Program ) 


lines » 
Statement functions 


Executable statements FORMAT 

Gat least one except ina | statements 

BLOCK DATA subprogram }) | (intermingled among the 
executable statements D 


END line @ obligatory in every case > 





his table @ devised by Colin Day ) prescribes an order more restrictive 
than that specified by Fortran 66, but an order more Jikely to result in 
a portable program. In Fortran 66 FORMAT statements are permitted the 
same freedom as CommeEnr lines. REAL, INTEGER, DOUBLE PRECISION, LOGICAL, 
COMPLEX, DIMENSION, COMMON, EQUIVALENCE ond EXTERNAL statements may be 
in any order among themselves. 
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EXERCISES CHAPTER 2 


Boll (ead what your Fortran manual says about COMMENT lines, 
continuation lines, \Sabels, END lines and the allowable order of 
statements. Note in the margins of your manual where there 
should be tighter restrictions so as to achieve portable programs. 


BoB Note below your particular “house rules” for writing: 


<= letter Z 
digit 1 
- letter 0 
digit 90 
= leter Z 
ss digit 7 


Qo Lt you are using punched cards, punch one with all forty-seven 
characters of the standard set and stick the card in the space 
below @ cut off the last thirty columns and it should fit }. The 
standard character set is shown overleaf = page 20. 
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UPS OF FoR 


CHARACTERS 
SYMBOLIC NAMES 

TYPES OF VARIABLE 
TYPES OF CONSTANT 
ARITHMETIC EXPRESSIONS 
LOG/CAL EXPRESSIONS * 
ASSIGNMENT 

LOANS QAN EXAMPLE ) 
EXERCISES 


* SKIP ON FIRST READING 





LETTERS , DIGITS & SYMBOLS 


"Phe Fortran 66 character set has forty-seven characters: 
@ the twenty six capital letters A to Z 
@ the ten digits O to 9 


@ a space or blank (as made by pressing the 
space bar of o typewriter keyboard once p 


@ the ten special characters tabulated below: 


eguals q or equals sign ) 

plus ( or plus sign ) 

minus ( or minus sign ) 

asterisk (or star ) 

slash @ or oblique or solidus ) 

left parenthesis (or left bracket » 
right parenthesis €or right bracket ) 
comma 

decimal point (or point or full stop ) 
currency symbol Cor dollar sign ) 
Although this is a standard character 
of Fortran there are some peripheral 
devices that cannot handle it © or 
print something different such as £ 


“here is implied by the order in which these characters are 
ere. 


presented 


Ml LCUNN*X I+ | 


ther characters ore permitted in COMMENT lines and in Hollerith items 
but they are best avoided if a program is to be fully portable. Some 
peripheral devices are not capable of handling every non-standard 
character that might be employed. 





FORMAT (8H%& &bo# <>) 
NON-STANDARD CHARACTERS LIKE @[1t!? etc 







SOME COMPUTERS 
DON’T HAVE ALL 
OF THESE 










‘WH 


[Fortran 77 adds the following two special characters: 


" apostrophe or quotation mark ) 
: colon 


Jy ower-case die. “small”) letters are not standard characters in 
Fortran 66 or Fortran 77 « but there are many sophisticated programs 
for processing words and the texts of books. It would be foolish to 
refuse to use (€or produce ) such programs on the grounds that they 
have to be written in non-standard Fortran. But they would not, of course, 
be fully portable programs. If standard Fortran seems to be letting us 
down on this point, be assured that other languages demand similar 
restrictions. Fortran @ with a few extensions ) is often the only language 
available for writing programs to handle words and letters. 
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PNT UTS INAS NAMES DEVISED BY THE PROGRAMMER 


Phe introductory example began: 


REAL DIAM, HEIGHT, COVRG, T®P, WALL, POTS 
NTEGER NPOTS 


and a little box was drawn for each variable: 


DIAM [__]__ HEIGHT[_]  Covrc[—__] 
gl es ak omen a. | ee 
wpors [I] 





[Zach variable is given a symbolic name @or just name) by the 
programmer as illustrated above. 


A symbolic name must begin with a letter and may contain anything 
from one to six letters and digits. @ Some Fortrans permit more than 
six letters and digits but for the sake of portability keep the limit 
to six. The limit in Fortran 77 remains at six. 2 


(nlere ore a few more symbolic names a programmer might invent: 


A 

N 

H2504 SKIP THE REST OF THIS 
M46 PAGE ON FIRST READING 
RESULT 





Symbolic names are for naming variables as olready illustrated. They 
are also needed for naming: 


Arrays q Chapter 5 ) 

Common blocks @ Chapter 8 ) 
Statement functions q Chapter 6 ) 
Functions @ Chapters 6&7 
Subrovtines @ Chapter 7 » 


A\lthough Fortran permits one name to be used for more than one thing 
(for example, naming a common block X and naming a variable X ) if is 
safer and less confusing not to do so. Likewise alihough there is no 
rule against using Fortran’s own keywords as symbolic names ¢€ for 
example using the word stop to name a variable ) it is safer not to do 
this either. Two good rules are: 


@ keep symbolic names of variables, arrays 
and statement functions unique in each 
program unit 


@ keep the names of functions, subroutines 
and common blocks unique over the com- 


plete Fortran program 
PERHAPS IT'S O-k. 
BUT I’M CONFUSED 


ye C STOP -£Q. IF +GO7T$) STOP Ligy 












GO7O 10 
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WARD (el WA Pee VN AU pre aloe rng PRECISION, 


he introductory example began: 


REAL DIAM, HEIGHT, COVRG, TOP, WALL, POTS 
INTEGER  NPOTS 


and a little box = a variable — was drawn for each item declared: 


pum (| omelanr []covrsQ KH 
i) eae (SSO (| Ca ep 


npors [i 


‘Whe size and precision of numbers fhat may be kept in these boxes 
varies from computer to computer. The magnitude of a fypical REAL 
may lie between —10% and +10°* and be represented to a precision of 
at least six significant decimal digits. A typical INTEGER may take any 
value from -32767 to +32767. On some computers these sizes and 
precisions are much greater. 


Six significant decimal digits may not offer enough precision in some 
calculations. For such cases Fortran provides DOUBLE PRECIS/ON variables 
which may have the same size range as ordinary REAL variables ( typ- 

ically t10% ) but ao precision of at least eleven significant decimal 


digits. 
D@UBLE PRECISION PI Aste tne 
ee te = 3.141592653589793D0 Fe arn 


TYPE 
PISS RS DOWALEE. PRECISION 


‘ihe simplest type in fortran is type LOGICAL. A _ logical variable may 
take only two values: ¢rue or false. 


LOGICAL OK, QUERY 


PK = .TRUE. 
QUERY = .FALSE. 


ok. |] query[ | 


Whe most complicated type in Fortran is type COMPLEX. Every complex 
variable consists of two parts: the real and the imaginary. Each part 
has the same range and precision as a REAL variable on the same 

computer. 


COMPLEX PPTEN Tae nae 
feel! POTEN = (123.45, —6.0) Add SEE PAGE 25 


POTEN TYPE 
es Gigm CD 
[Fortran 66 refers to storage units of which: 
@ ore is needed to store an JNWTEGER, 
a REAL, OQ LOGICAL variable 
@ two ore needed to store a DOUBLE 
PREC/SION Or COMPLEX variable 
but in practice different Fortrans employ computer “words” in entirely different 
proportions. Do not rely, for example, on a REAL variable and an /W7EGER variable 
occupying the same space: they often don’t. 
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A variables may have their types declared as illustrated below: 


REAL X,Y, 2, AA 
INTEGER I, IA, KOUNT 


D@UBLE PRECISION PI, DX4 





but it is not essential to declare INTEGERS or REALS. When the type of a 
variable is not explicitly declared it is implicitly declared as type 
INTEGER if its initial letter is: 


I,J,K,L,M,N 


otherwise it is implicitly declared as type REAL. This convention is 
sometimes called implicit typing. 


‘Whe forms of type statement are: 


REAL same, name,...,name 

INTEGER same, nameé,..., name 

DOUBLE PRECISION same, name,..., Name 
LOGICAL rame, name, ..., name 

COMPLEX name, name, ..., name 


where: 
mame is a symbolic name of a variable or of 
an array ( see Chapter 5). When a sym- 
bolic name refers to an array then the 
dimensions of that array may also be 
given as though the name were in a list 
folowing a DIMENSION statement 


When ao variable is named as above this does not imply that it has an 
initial value of zero; its content is completely undefined; there is no 
soying what the little box might originally contain. 


Some Fortrans (@ including Fortran 77 ) provide a statement called JMPLICIT 
using which the programmer may alter the range of initials denoting types 
of variable. This facility shovld not be used in a program intended to be 
portable. 


Qome Fortrans (including Fortran 77) provide two kinds of integer : short 
and long. A short integer is denoted INTEGER*2 because it consists of two 
“bytes” each of eight binary digits; its range is t2°-1. Similarly the long 
integer is denoted JN7TEGER#+4 and its range is t2%-1. A program that 
relied upon the long integer would not be fully portable. 


[Fortran 717 provides an extra type called type CHARACTER. If a program is 
to be portable, however, it is better to manipulate characters by storing them 
in variables of type INTEGER as discussed in Chapter 9. 


‘WPhus there are just five types of variable in Fortran 66. There are also 
the same five types of array € Chopter 5) each being an array of 
elements. In general what can. be done with a variable can also be done 
with an array element of like type. There are, however, some important 
exceptions to this rule and attention is drawn to such exceptions wherever 
they occur. 
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REAL, INTEGER, DOUBLE PRECISION, 
LOGICAL, HOLLERITH, COMPLEX 
Phe introductory example had the line: 


TOP = 3.14 * (DIAM ¥*2) / 4.0 
2 4 


where the arrowed items are constants of type REAL. REAL constants are 
written with decimal points. Thus a RFAL constant with a value of two 
may be written 2.0 or even as 2 but never as 2 because 2 is an 
INTEGER constant as explained below. A REAL constant of o half may be 
written as 0.5 or just .5 without the leading zero. 





“A REAL constant may also be written in exponent form as illustrated below: 


[|7OP = 0.31461 *( DIAM#*2) / 400.08-2 
He ry 


where letter £ says “ times ten to the power of...” then must come an 
integer to specify the power. The example above shows 0.3/4£1 which 
means 0.314x10! and 400.0£-2 which means 400.0x10-% . One million couid 
be written in many ways including 1.0&6, 10.0£F#5, 10.E5, 0&5. 
Notice there need not be a decimal point in a REAL constant if written 
in exponent form € 10&5 ). 





SPhe introductory example had the lines: 


NPOTS = INT(POTS) + 1 Stee 
WRITE (6, 200) WNPOTS CONSTANT 
be 


Vas 
where the arrowed items are constants of type INTEGER. INTEGER constants 
consist only of digits. There should be no decimal point and no letter £. 
An INTEGER constant of thirty thousand is written 30000 and not 3&4 which 
is one of the forms of a REAL constant. 













A DOUBLE PRECISION constant has the same form as a R£AL constant written 
in exponent form except that letter D replaces letter £. 


DOUBLE PRECISION DIAM, TOP 
TOP = 3.14/592653589793D0 * (DIAM*#2)/4.0D0 
A 


z} 
DOUBLE 
PRECISION 
CONSTANT. 


Go far constants of type JNTEGER, REAL Gnd DOUBLE PRECISION have been 
defined without mention of a sign @+ or-) as a prefix. Without the 
sign these constants are said to be unsigned : with a preceding + or—- 
they are said to be signed. If there is no mention of signed or unsigned 
then the sign is optional. For example “An integer constant of two” may be 
written 2 of +2 ; “An integer constant of minus two” should, of course, be 
written -—2 . 














DOUBLE 
PREC/S/ON 
CONSTANT. 






Shere are only two constants of type LOGICAL. These ore written .7RUE. 
and .fALSE. as illustrated below: 


LOGICAL  @K, QUERY 


OK = .TRUE,. 
N59 LOGICAL 
QUERY. = FALSE. Soe iSischurs) 
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lermann Hollerith pioneered punched-card machines in the 1890s. His 
initial, H, survives in Fortran 66 to denote a sequence of characters as 
illustrated in the introductory example: 


| |FoRMAT (1x, 9H YOUNEED, I2, SH POTS ) 


& eerie” wees 





The unsigned integer before letter H specifies the number of characters in 
the literal immediately following the 4. This number includes spaces. In 
the above example there is a space before and ofter You and a space 
before POTS. In a FORMAT statement these items are called Hollerith 
literals. There are also, in Fortran, Hollerith constants which have 
precisely the same form as Hollerith literals: 


but may be stored in variables of type INTEGER. In fact Fortran 66 
permits Hollerith constants to be stored in variables of any type and 
Specifies no limit to their length = length being a property set by the 
particular computer being used. To make a program fully portable, 
however, confine Hollerith constants to INTEGER variables and do not store 
more than two characters per variable. Storage of characters is intro- 
duced in Chapter 9, and an exomple of fheir manipulation in Chapter 12 . The 
examples should demonstrate that the restrictions advocated above do not 
prevent the programmer from manipulating characters effectively. 


We anticipate Chapter 9 with the following example by which the Hollerith 
constant 2HAB (Qz.e. the letters AB) may be stored in the integer variable 
named K : 


[ INTEGER K 


DATA K / 2HAB/ 
CJ = 





‘ 2 aie ; 
and this may m= be achieved by assignment: 
—/— 


pi 





Some Fortrans allow Hollerith literals and constants to be written between 
quotation marks «= dispensing with the count and letter H: 


200 FORMAT (1X, ' yOU NEED’, I2, ' PQTS’ ) 


which is certainly more elegant @ and is the form preferred by Fortran 77) 
but shovid not be used if a program is to be fully portable. There is no 
quotation mark in the Fortran 66 choracter set. 


@onstants of type COMPLEX ore written as two REAL constants separated 
by Q@ comma and enclosed between parentheses: 





where the 2.7 constitutes the real part and the -0.9 constitutes the 
imaginary part. Mathematicians would write the above complex 
constant as 2.7-0.9i where the i stands for w-=i € the square root 
of minus one}. Non-mathematicians need not bother with type ComPLex. 
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REAL, INTEGER, DOUBLE 
PRECISION, COMPLEX 


ADRES ERE 


Phe introductory example showed: 










TOP = 3.14 * CDIAM*&2)/ 4.0 
WALL = 3.14% DIAM * HEIGHT 
POTS = (TOP + WALL) / COVRG 


Where in each case the expression on the right of the equals sign is 
evaluated and the result is assigned to the variable nominated on the 
left. 








[2xpressions may consist of variables, constants and other terms* bound 
together with operators. The arithmetic operators in Fortran are: 


4 for addition <= or as a prefix 
to denote a positive quantity @eg. A=+B) 


= for subtraction« or as a prefix 
fo denote a negated quantity @e.g. Z=-7 D 


* for multiplication 
7. for division 


* x for exponentiation 


“An expression is evaluated generally from left to right in three sweeps. 
Exponentiation is done in the first sweep; multiplications and divisions 
@with equal precedence 2 in the second sweep; additions and subtractions 
with equal precedence) in the third sweep. 


(Parentheses may be used to emphasize or change the natural precedence 
described above « and achieve what one expects from the rules of algebra. 
Notice that the parentheses in the first expression above are not necessary: 
DIAM#*2 would automatically be taken first. But parentheses are vital in 
the third line: POTS = (TOP + WALL)/ COVRG . 


On general, all terms in an expression should be of the same type. 
There is a notable exception in the case of exponentiation illustrated 
in the introductory example: 


| |79P = 3.14 * (DIAMRE2) / 4.0 
a 


where DIAM is of type REAL and the power 2 is of type JNTEGER. Fortran 
66 allows an arithmetical term of any type @ ZNTEGER, REAL, DOUBLE 
PRECISION OF COMPLEX Q to have an INTEGER exponent ; the result of such ex- 
ponentiation having the same type as the term being exponentiated. 





Qn the above example it would be correct to have DIAM*#2.0 in place of 
DIAM*#2 but this might force the Fortran compiler to evaluate this term by 
generating logarithms instead of multiplying DIAM by DIAM. 


*the other terms are “function references” (illustrated opposite) and “array elements ”(ch.5 ). 
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A DOUBLE PRECISION term may be exponentiated using a R£AL exponent ora 
DOUBLE PRECISION exponent «= the result being DOUBLE PRECISION in each case: 


a 





DBL**I~"+ DBL**¥R + DSLevD 


Fortran 66 also allows combination of REAL and DOUBLE PRECISION terms in other 
ways «= the result being DOUBLE PRECISION. Fortran 66 also allows combination 
Of REAL and COMPLEX terms = the result being COMPLEX. Examples follow. 










REAL 
DOUBLE PRECISION DBL, D 
COMPLEX CMPLX, C 
BL = R*¥D+D/R 

CMPLX = R*#C4C/R 


A\lthough many Fortran compilers permit other mixtures of type their details 
and interpretations differ. So for the sake of portability do not take 
advantage of “mixed mode” facilities apart from those just described. 


“Phe introductory example showed: 
pO BED) 
( 


where INT() is a Fortran function for converting a REAL quantity into an 
INTEGER by truncation @ for example IN7T(3.995) yields an integral value of 3 Dp. 
Thus both terms in the expression above ore of type INTEGER. Fortran pro- 
vides a range of such functions for converting expressions and parts of 
expressions from One type fo another. These functions @ defined in 
Chapter 6 ) make if unnecessary to rely on non-standard “ mixed mode” 
facilities. 






There ore some pitfalls to be avoided when writing expressions. In dividing 
integers any fractional part is lost. In the example below Z and J would 
be assigned values of 2 and —2 respectively: 


J = -7/3 


and there is no implied multiplication. The expression below should have 
an asterisk between the two terms: 


| [JA = ce-ae-F) } 
3 


and it is not permitted to have two operators together without a bracket 
intervening. The expression below should read B*(-c) or -Bec: 


¥ 


Sorne expressions might be ambiguous to the reader: 
le = B/c/D 


P = QaeRees a 


Xx = -Y¥«x*Z 





There is a world of difference between 8/(4/2)=4 and (8/)/2=1. Fortrane6 
Should treat the first example as A=((B/c)/D) but there might be a compiler 
somewhere that would treat if the other way. The second example illustrates a 
forbidden form, but either (@Q#*R)**s or @=#(Re¥S) would be allowed. Fortran 
should treat the third example as xX=(-(y##zZ)) but there might be a compiler that 
would try evaluating X as (-Y)s*z and get into trouble. When in doubt use brackets ! 
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eT ATOANG LEN STINGS asiie sits pace On resr READING) 


“An expression of type OG/CAL is one that has the Boolean value ¢rve or false. 
The main use of logical expressions in Fortran is in the logical IF statement 
described in Chapter 4 = hence the invitation to skip this double page on 
first reading. 


Whe simplest logical expressions are the logical constants: 


. TRUE. 
FALSE. 


but a logical expression may be more complicated = consisting of 
relational expressions bound together using logical operators. 


A relational expression consists of two arithmetic expressions bound to- 
gether by one of the following six relational operators: 


-LT. meaning “less than” 

-LE. meaning “less than or equal to” 
-EQ. meaning “equal to ” 

-NE. meaning “not equal to” 

-GT. meaning “greater than” 

-GE. meaning “greater than or equal to” 


An example of a relational expression involving type RéAL is: 


A*B/C .LE. X+45.0 


which takes the Boolean value trve if the numerical value of A*B/c turns 
out to be less than or equal to the numerical value of X+5.0 . If the value 
of A*B/c turns ovt not to be less than or equal to that of xX+5.0 the 
expression takes the Boolean value false. 


An example of a relational expression involving type JNTEGER is: 


Shere is no such thing as a relational expression involving type COMPLEX. 


Ot is wrong to compare an expression of type JN7EGER with an expression 
of any other type: 


: 


however it is admissible jin Fortran 66 fo compare an expression of type 
REAL with an expression of type DOUBLE PRECISION: 


but such practice can be dangerous because you may not be able to tell 
in which mode the comparison is made. Is X converted to double precision 
before comparing with 1.5D6 ? Or is 14.5D6 converted to standard precision 
before being compared with X ? In such cases it is best fo use one of 
the intrinsic functions @ Chapter 6) to make certain about the mode of 

comparison: 







DBL(X) .£ 


CONVERTS X 


3 
g8 4 
iS 


PRECISION 
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She simplest logical expressions @ as said before) are the logical 

constants .7RVE. and .FALSE. = but the expressions may be much more 
complicated, consisting of relational expressions @ defined opposite ) 
bound together with logical operators. The three logical operators are: 


-OR. = for “logical disjunction ” 
-AND. for “logical conjunction ” 
-NOT. for “logical negation ” 


[2xamples of logical expressions are: 
«TRUE. 

















A RELATIONAL EXPRESSION 














- NOT. . TRUE. ON ITS OWN 
I*xZ2 .£Q J RELATIONAL EXPRESSIONS 
I EQ. J OR. I .EQ. K “UA LOGICAL OPERATOR 





An the fourth example immediately above, the logical expression takes 
the Boolean value é¢rve if I equais J or if Z equals K or if I,J,x« 
are all equal to one another. In other words the .or. is inclusive 
rather than exclusive. 


SPhe binding strength of .wo7. is stronger than that of .AWD. which, in 
turn, is stronger than that of .oR. ; thus: 


-NOT. BLACK .AND. WHITE 
means (N07. BLACK) .AND. WHITE rather than .NOT. (BLACK .AND. WH/TE) . 


Brackets may be used in the construction of logical expressions @ and 
to add clarity to expressions that do not really need brackets ) just as 
with arithmetic expressions: 


CX .GE. 1.0) .AND. (X .LE. 2.0) 


This expression tokes the value ¢rue if the value stored in variable X lies 
between 1.0 and 2.0 , 


Qn the absence of brackets Fortran deals with relational expressions first, 
then applies the logical operators € first .wo7. then .AWD. then.OR ) to the 
resulting Boolean values. Make sure you don’t leave a relational 
expression “dangling” : 


X .GE. {0 .AND. 


at LE. 2.0 
DOES NE INTEND 
AN X HERE? 


| 4 LOGICAL EXPRESSION 
pal) <ead_8u7 Nor FORTRAN, 
a 
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STATEMENTS WITH AN EQUALS S/GN 
(EVERY ASSIGNMENT JS AN EXECUTABLE STATEMENT ) 
“Whe introductory example showed several assignments: 


TOP = 3.14 *( DIAM*¥#2)/4.0 <x, 
WALL = 3.14 * DIAM* HEIGHT 
POTS = (TOP + WALL)/COVRG << 
NPOTS = INT(PQTS) + 1 


where the arithmetic expression on the right of the equals sign is evaluated 
and the result assigned to the variable nominated on the left « obliterating 
any previous value found there. 
















“She first three assignments above are of type REAL: the result of a REAL 
expression being assigned to a REAL variable. The last assignment above is 
of type /NTEGER: the result of an /NTEGER expression being assigned to an 
INTEGER variable. The piece of program below illustrates DOUBLE PREC/S/ON, 
COMPLEX and LOGICAL assignments: 


DOUBLE PRECISION D, DBL 
COMPLEX C, CMPLX 


DOUBLE PRECISION 
ic a pe 


BL= 1D6x D 


eMPLX = (1-5, ~1.0) xC 










[for the sake of clarity it is best to ensure the expression on the right of 
the equals sign is of the same type as the variable* on the left of it. 
However, Fortran 66 does permit a change of type across the equals sign 
in the case of J/NTEGER, REAL Ond DOUBLE PRECISION. The allowable forms of 
assignment are: 


INTEGER variable” gz | INTEGER expression 
REAL variable * wm me < REAL expression 
DOUBLE PRECISION variable*|\% | povBLe PRECISION expression 


COMPLEX variable * = COMPLEX expression 


LOGICAL variable* = LOGICAL expression 


Qn the first of the three forms depicted above the Fortran compiler automatically 
converts the result of the expression (if necessary ) to the type of variable 
to which the value is assigned. Thus it would be permissible to write: 


|__| | wrers = POTS + 1.0 
Uinreac®) Eee) Vea) 


[wears = qwr(rovs)+i 
[inate Ranoce) Wneaee 


Warning / Do not assume the computer would preserve the full accuracy of 
the product in an assignment such as: 


ce ae eae aa TATE 


where D is Q@ DOUBLE PRECISION variable and RA and RB are REALS. But you may 
ensure full accuracy by changing the form of the assignment to: 


| | [a= DBL (kA) * DBL(RB) Cie caonRGroR DOD) — | 


% “variable” is here taken to include “array element”. 


in place of: 
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ROARS AN EXAMPLE TO ILLUSTRATE ARITHMETIC ASSIGNMENTS 
‘Whe monthly repayment, $c, on a mortgage loan of $s at P% 
compound interest over N years is given by: 
SR(i+R)" 
wl(i +R)" =1] 
R = P + 100 


where: 


(Alere is a set of data for a program designed to compute C: 





A\na here is the program itself: 


READ(5, 100) $,N,P 
R= P/ 100.0 


C = S*R*(1.0+R)*#N /(12.04((1.04R) #4#N — 1.0)) 
WRITE (6,200) C 

TOP 
FORMAT (F 10.0, I10, Fi0.0) 
FORMAT (1X, 19H MONTHLY BURDEN = #, F6.2) 
END 





Ana here is the resulting output: 
° 


° Py 
o MONTHLY BURDEW = $255.45 
° 





Gome innocent looking assignments can cause unexpected trouble: 


A= B/c 
X= YueZ 


If variable C holds zero = or a number near to zero = there would probably 
be an error message printed and execution would cease. Variable X would be 
set to 1.0 if Z were zero and Y greater than zero @e.9. X=6.5%«0 ) but if 
Y contained a value of zero or less the program would probably fail. 
A negative value cannot be raised to a non-integral power = this is a 
mathematical rule having nothing to do with Fortran specifically. 


2.1 
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CHAPTER 3 


Bog Read what your manual says about the set of allowable 
characters and strike out those not in the standard set of 47. 


PS ec Le your manual permits names to consist of more than six 
characters, make oa marginal note that no more than six should 
ever be used. 


309 Note in the table below the size and range of each type of 
variable on the computer to which you have access 


[type ___—— number of bits | range precision _ 
mene ee ee 
2 SG REE” (SP aeen Gaeaies 
CTT ES a, PE: Sees 
ISR ie OE S| EERE Se, 
SEE EAE RE 








COMPLEX 


“« ” 

3 ob FPortran is an acronym for “formula translation”. Take a 
formula used in your own field of expertise. Write a small 
program (like the one on the previous page ) to read data; 
evaluate your chosen formula; print the result. 
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CONES WihhlN 
DN ARORA ONT 


SIMPLE LOOPS 

SHAPES (OR STRUCTURES ) 
LOGICAL IF 
UNCONDITIONAL TRANSFER 
COMPUTED GO TO 
CONTINUE 

THE DO LOOP 
ARITHMETIC IF 

ASSIGNED GO TO 

AREAS OF SHAPES QAN EXAMPLE ) 
EXERCISES 


INTRODUCING THE GO TO AND 
LOGICAL IF STATEMENTS 


ere is an example of a program that does nothing but read an integer, 
print it, then repeat the process until there are no more integers to 
read « when an error message from the computer would be printed to 
indicate’ an unsatisfied READ statement. 


READ (5,100) INT 


WRITE (6, 200) INT be 
GO TO 10 INFINITE 


Loop 
STOP > a 
FORMAT (110) 
FORMAT (1X, I10) 
END 


Notice the Go 70 statement which causes control to pass to the statement 
having the nominated label. This particular 60 70 would prevent the STOP 
statement ever being obeyed. 





Ls the first item of data were a count of the number of integers following: 


SAYS TWAT 3 
INTEGERS FOLLOW 
IN THIS, CASE 





then the program could be rewritten as follows: 


READ (5, 100) K®UAT 
10] |IF (KOUNT .£q. 0) GO TO 20 
READ (5,100) INT 
WRITE (6, 200) INT 
KQOUNT = KQOUNT - 1 





and this program wovid stop properly even if the count of integers were zero. 


AXlternatively the program could be written like this: 


READ (5,100) K@UNT 

READ (5,100) INT 

WRITE (6,200) INT 

KOUNT = KOUNT — 1 

IF (KOUNT .NE. 0) GQ TO 10 
STOP 


but this program would fail if the count of the number of following integers 
were zero. Statement {0 cannot be avoided so there must be al least one 
integer to be read and printed. 





ZAnother way to tackle the problem is to omit 

the count of integers from the data but ferminate 
the list of integers with a zero @ assuming, in 

this case, that the program is designed to 


handle only non-zero integers ). The program 
could be written like this: CmIsey 


10||READ (5, 100) INT 
F (INT .£Q. 0) GO TO 20 
WRITE (6, 200) INT 
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ST/CK TO THESE SHAPES OTHERWISE 
OR “STRUCTURES” YOU RISK GETTING INTO A MESS 


She last three programs opposite display the shapes @or structures) defined 
ge the following flow charts: 


THE 

| statements | HYBRID | statements _| 

Loop ars Loop b> <Gonin> 
ae =r F statements | LOOP ee nar 


fda to these loops the following two shapes: 





Di? oo 





Qsing these shapes @even omitting the hybrid | 
loop » it is possible to program anything C=R 
capable of being programmed. And the 

program would be properly structured. 


Whe box called statements may be empty. For example JF..7HEN..£LSE may 
be drawn thus: or thus: 


THEN it ELSE 
| seataeneats 






eS 


i 


Gonversely any box called statements may contain very many Fortran 
statements to be obeyed = including those comprising the shapes illustrated. 
In other words all these shapes may be nested. Here is an illustration of 
nesting: 





LOGICAL nF THE MOST USEFUL CONTROL STATEMENT 


Fortran offers some special methods of control @ considered later ) but 
the shapes previously sketched may all be achieved using the LogIcAL IF 
statement. Its form is: 


IF ( logical expression) statement 
where: 


logical expression is an expression whose value is 
either true or false as described 
On pages 26.and 29. 


Statement need not be the 60 7o statement 
previously illustrated. It may be 
any executable statement other 
than another L06/CAL /F or A DO 
Statement @ see page 40). 


ere ore some examples of LOGICAL IF statements: 











IF (0K) GO TO 66 
F CI .NE. J) STOP 
IF CA®B/C .LE. X#5.0) AmAt 1.0 


IF ((X.GT.¥) .OR. (CX .LT. SORT(Z))) THT 





When such a statement is obeyed the logical expression is first 
evaluated. If the result is the Boolean valve true then the statement 
incorporated in the LOGICAL IF statement is obeyed. If the Boolean value 
turns out to be false then control passes straight to the statement 
immediately following the LOG/CAL IF. 


A subtle point arises when the logical expression is evaluated. In the 
statement labelled 66 above the computer might compare values of X and 
Y; discover that X was greater than Y ; and so obey the statement Z=-J 
without bothering to compute the square root of Z. That’s all right. 
But if the function were one devised by the programmer @ Chapter 7) and 
which changed the valve held in J, what then? The result would probably 
vary from one Fortran to another. Never rely on such an expression 
being fully evaluated. 


pposite is Q@ program designed to illustrate some nested “shapes” 
previously defined. The program reads a number of integers to discover 
the maximum and minimum. The data are written after a count of 
the number of integers to follow « precisely as in an earlier example. 





Rlere are the data: 


COUNT OF INTEGERS 
FOLLOWING 
CAT LEAST ONE) 





Prlere is the program. It will cope with any number of integers but no 
integer may exceed a million in absolute value: 













INTEGER MAX, MIN, INT, KOUNT 
READ (5, 100) K@UNT 
MAX = -1.0E6 <EG NO ITEM OF DATA MAY 

MIN = 41.0E6 sty FALL OUTSIDE THE RANGE 


t 1,000,000 IN THIS 
READ (5, 100) INT PROGR. 




















? 
OGRAM 


IF (INT .GT. MAX) MAX = INT <xe 

IF (INT «LE. MIN) MIN = INT og ue 
K@UNT = K@UNT - 1 

IF (KOUNT .NE. 0) G@ TO 10 


WRITE (6, 200) MAX, MIN 
100| |FORMAT ( I 10) OF 10 COLUMNS 


200} FORMAT (1X, 4H,MAX, I7, ‘a MIN, I?) 
















&\nd the output would be: 


MAX 217 MIN -73 
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NNO MONG ub Bas Suse 7” 


“Thus far this chapter has introduced the GO 7o statement by example 
only. The form of this statement is: 


GO : TO label 


label is the label of an executable statement 
in the same program unit as the Go 70. 
This item should consist only of digits 
(the name of an integer variable is not 
allowed ). 


where: 


7 hen this statement is obeyed control is passed to the statement bear- 
ing the nominated label: 





She introductory example illustrated the S70P statement. Its form is: 


STOP octal number 


where: 
octal number is optional. If included it should be 
limited in length to 5 octal digits. 
qOctal digits are those from O to7, 
but Fortran 77 permits decimal digits 
from 0 to 9 p 


When this statement is obeyed execution of the program ends normally. It 
depends on the computer installation what happens to the octal number, 
but that number should somehow be conveyed to the person trying to 
run the program ~ perhaps by being sent to unit 6 which is convention- 
ally the line printer. There may be any number of S70P statements ina 
program; the programmer might be interested to know which of them 
caused his program to stop. 


[otice that the s7oP statement is the standard way to cause a program 
to stop execution. The END line has an entirely different purpose as 
explained on page 14. The CALL EX/T statement ~ available in many 
Fortroans «= is not a statement defined by Fortran 66 or fortran 77. 


[Fortran 66 and Fortran 77 both recognize a PAUSE statement: the word 
PAUSE optionally followed by an integer @ octal in Fortrané6 ). This statement 
may have been useful in the days when the programmer had the computer 
all to himself. He could note the octal number displayed by the computer and 
decide whether or not fo cause the program to start execution again = 
carrying on from the statement immediately following the PAUSE. But on 
large modern computers the PAUSE statement is either ignored or dealt 
with according to rules local to the installation. So the PAUSE  state- 
ment is best avoided if programs are to be made fully portable. 
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COMPY TED (ee) aT.) “CASE 36 OR oun” 


[Fortran provides a statement called the computed Go To which causes 
control to pass fo one of several labelled statements according to the 
value currently stored in an integer variable. The form is: 


GO TO (label, label, label, ... ,label), variable 
GD 


label is the label of some executable statement 
in the same program unit as the computed 
GO To. There should be one or more of 
these inside the parentheses. Each may 
consist only of digits; the name of an 
integer variable is not allowed. 


where: 


variable is the name of oa variable of type JW7EGER. 
@An array element is not allowed here.) 


When this statement is obeyed the value stored in the integer variable is 
consulted. If this value is unity, control passes fo the statement labelled 
with the first label written between the parentheses; if the valve is two, 
control passes to the statement labelled with the second label between the 
parentheses = and so on. 


GO TO (10, 10, 20, 10), IMP 


A=xX 
GO 76 30 


A= Y 
CONTINUE 


In this example control would pass to statement 10 if the value stored in 
JMP were ior 2 or 4, buf to statement 20 if the value in JMP were 3. 





What if the value in JMP turned out to be negative or greater than 4 ? 
The answer is that different Fortrans do different things. It is best fo 
protect the computed GO 70 as illustrated below. 


a IF CCIMP .GT. &) -OR. (IMP .LT.1)) STOP OFF es 
G® 7T@ (10, 10, 20, 10) IMP PAae 


CONTINUE AN EXECUTABLE STATEMENT ALWAYS LABELLED 


Phe CONTINUE statement is a “do nothing” statement of the form: 
CONTINUE 


‘Phis statement is classed as an executable statement even though it causes 
no work fo be done. It is a useful statement = when labelled ~ to act 
as the terminal statement of a DO loop @ page 40) or as a point ina 
program where several paths converge, as in the example above at 
label 30. 


“Although Fortran 66 does not forbid it, a@ CONTINUE statement without a 
label is pointless « and probably points to a mistake in the program. 
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Snll2 LOOP “REPEAT UNTIL” 


[Fortran provides a special form of repeat until structure called the DO 
loop. Its form is: 


DO label control = initial, terminal, increment 
DO label control = initial, terminal 


or: 


where: 
label is the label of some executable statement 
physically following the DO. This item must 
consist only of digits. 


control is the name of an integer varioble @ not an 
array element ). 


initial is either an unsigned integer or the name 
of an integer variable @not an array element ) 
which must hold a value greater than zero 
when the Do is executed. 


terminal is of the same form as initial and must hold 
a value greater than or equal to that of initial 
when the DO is executed. 


increment is of the same form as initial and terminal 
and must have a value greater than zero when 
the DO is executed. If omitted, increment is 
assumed to be unity. 


Inlere is part of the “max. & min.” program again, but using the DO loop: 


READ (5,100) kOUNT 
MAX = -1.0E6 

MIN = +4#1.0E6 

DO 10 I= 1, KOUNT 

READ (5,100) INT 

IF (INT .GT. MAX) MAX =INT 
F (INT .LE. MIN) MIN @ INT 
6, 200) MAX, MIN 









On meeting the DO the computer 
assigns the énitial value to control 
din this case the integer variable I 
is assigned a value of 1). Then the 
computer obeys subsequent instruct- 
ions down to and incivding the 
one with the nominated label @ in 
this case the statement labelled 10). 















Shen the control variable is incremented by the valve held in increment 
qin this case unity by default D. If the augmented control variable now 
holds a value greater than that specified by terminal the loop is satisfied 
and control passes to the statement immediately following the labelled one; 
otherwise the loop is executed again. The word DdO is short for “ditto”. 


Rote carefully the restrictions on the components of a 300. It is wrong to 
count backwards: 


| sae 10 =r = KouNT, 1, -1 


and it is wrong to count from zero: 


[— [be 00 t-0,9 med 


although there are some Fortrans that permit such abuses. Many Fortrans 
Obey any DO /oop they encounter at least once « even when the controlling 
parameters are in conflict: 


a |, 
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SPhe label must be the label of an executable statement other than: 

(i) arithmetic IF ( page 42, Gi) RETURN € page 6p, (ii) STOP, (iv) PAUSE , (v) GOTO 
or (vi) another DO statement. The terminating statement may, however, be a 
logical IF @as already illustrated ) provided that this does not incorporate 
any of the forbidden statements (i) to (vi) set out above. 


4 woLl= eR Ow SE Ly™ ee Soe RO Baron) 





LAany programmers moke it a rule to finish every DO loop at a CONTINUE : 


Dé 10 L=1, KOUNT 

READ (5,100) INT 

IF (INT .GT. MAX) MAX=INT 
TF (INT «LT. MIN) MIN =INT 





Mever tamper with controlling parameters inside a loop: 


, 
WRITE (6, 300) K,L 


FORMAT (1X,1I3, 10H TIMES 5 =, I3) 





fAne it is safe to jump out of a loop before it has run its course * in which 
case the control variable retains its current value: 


Dd 10 I t, 100 
READ (5, 400) X 
ZF (X .£Q. 0.0) G6 TO 20, ow 


WRITE (6, 500) X THE NUMBER OF NUMBERS 
READ — HAVING FOU! 
ahaa ZERO AMONG THE DATA 


WRITE (6,600) I 





ut never assume anything about the value of a control variable once the loop 
has run its course. In the example below the printed value of XK might turn 
out to be 4; perhaps 3; perhaps some arbitrary value. 





Some Fortrans allow a jump back into the middie of a loop = the statements 
executed whilst being out of fhe loop constituting an “extended range”. But 
some Fortrans @ including Fortran 77 ) object to this practice so it is best 

not to employ extended ranges. Never jump into the middie of a Da loop/ 


Dd 10 Z= 1, 3 
GO TP 20 
CONTINUE 
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NOW LARGELY SUPERSEDED BY 
THE LOGICAL IF STATEMENT 
Fortran provides a unique three-way switch of the form: 
IF (expression) label, label, label 
where: 


expression is an arithmetic expression of type INTEGER 
or REAL or DOUBLE PRECISION 


label is the label of an executable statement in 
the same program unit as the /F statement 


When this statement is obeyed the expression is evaluated: the result will 
be negative or precisely “zero” or positive. Control then passes to the 
statement having the first, second or third label respectively according to 
this result. | 


Terr (x+y) §=f0, 10, 20 
10) |Z = J <4 OBEYED NEXT IF 
GO TO 30 Pe seo 
20 <E-§ OBEYED NEXT IF 
30 (X4#¥) > 0.0 


‘Whe “zero” referred to above is the kind of zero appropriate to the 
expression being evaluated. Thus the zero would be: 


0.0 for a REAL expression @ e.g. (x+y) D 
fs) for an INTEGER expression 6.9. (M+N) ) 
0.0D0 for Q DOUBLE PRECIS/ON expression @ é.g. (DX+1.0D0) ) 


Whe arithmetic IF was once the only IF statement in Fortran. Its use has 
been largely superseded by the logicaé IF which mokes a program easier 
fo read and understand. 


A four-way switch can be achieved using a logical IF statement incorp- 
orating an arithmetic IF statement: 


IF (X+tY .G% 100.0) IF (xt¥) 10, 20, 30 
I=WN 


OBEYED NEAT IF 
(XtY) © 100.0 





but this is the sort of “clever” programming best avoided. 


SOME TRAFFIC & )} 
OTHER TRAFFIC S 
MORE TRAFFIC = 
LESS TRAFFIC $ 





“CASE C OF...” 
BUT IT 1S BEST NOT TO USE IT 


[Fortran 66 provides an alternative to the computed GO TO mechanism. This 
involves special assignments of the form: 


ASSIGN label TO variable 
where: 


label is the label of an executable statement in 
the same program unit as the ASS/GN 
statement. This item consists only of digits. 


variable is the name of a variable @ not an array 
element ) of type JNTEGER 


This assignment is used in conjunction with a special GO 70 statement of 
the form: 


GO TO variable , (label, label, ..., label ) 
where : 


variable is the name of an integer variable to which 
a valid label has qwhen the Goro is 
obeyed ) already been assigned by an 
ASSIGN statement. 


label is as defined for the ASSIGN statement. The 
parentheses surround a list of labels which 
must include that which has already been 
assigned to variable when the G0 70 is 
obeyed . 


Were is an example of the assigned GOTO in use: 


INTEGER LABL 

ASSIGN 10 T@ LABL 

ASSIGN 30 76 LABL 

IF (IZ .EQ. J) ASSIGN 20 T@ LABL 


Aehess OF EXtcurabie sianemens) 
| __—*{|6e 70 LABL , (10, 20,30) 


When the ASS/GN statement is obeyed the integer variable is assigned some 
“value” @a pattern of binary digits ) corresponding to the label of a state- 
ment. ASSIGN 10 TO LABL is {notk the same thing as 4ABL=10. So LABL 
should not be referred to save by ASS/cN statements. When GO TO LABL 
is obeyed contro! is transferred to the statement bearing the label last 
assigned to LABL. This must be one of those in fhe bracketed list. 


Different Fortrans have different restrictions and extensions to the 
assigned GO TO mechanism so this form of control is best not used 
in programs intended to be portable. The computed GO 70 described 
on page 39 provides a safer and better mechanism for dealing with 
the shape “case C of...” ., 
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AN EXAMPLE TO ILLUSTRATE 
CONTROL STATEMENTS 


[ere is Q program designed to calculate areas of rectangles, 
triangles and circles. 





Whe data comprise two lines @recordsp) for each shape. The first line 
is given as 1 or 2 or 3 according to the shape required: the second 
line contains two or three or one dimension(s) appropriate to the 
Shape selected. Thus the data could look like this: 
<4) [TRIANGLE 

<HA,BE&C 
<I RECTANGLE 

HS BAD 
<LY ANOTHER TRIANGLE 

HA,BEC 
<—N CIRCLE 

ND 
<EJ END OF RUN 














ZA flow chart for the program is shown below: 


READ 1,2 OR3 
ACCORDING TO 
SHAPE SELECTED 











READ A, BE&C 
COMPUTE 
S=Z(AtB +c) 


COM PUTE 
z =,/S(S-AXS-8X3-C) 













COMPUTE 
z= BD 
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“Whe program itself could be written thus: 


REAL A,8,C,D, zZ,5 << sercmchsow STATEMENTS 
INTEGER ITYPE Se i Pe) 
READ (5, 100) ITYPE 

IF CCITYPE .LT. 1) .@R. (ITYPE .GT. 3)) STOP 
GO TH (10, 20,30), ITYPE 
READ (5, 200) B,D 
Z=B*eD 

G6 TO 40 

































RTC) 13, A 
parece FUNCTIO 
(CHAPTER 6) 





—_——_ 










<t-f INTEGER /N FIRST 
COLUMN OF CARD 






FORMAT ( 2F 10.0) 
FORMAT (3F 10-0) 
FORMAT (F 10.0) 
(1X, 9H AREA IS , F#.2) 







Because there are only three possibilities specified in the statement 
GO TO (10,20,30), ITYPE , this line could be replaced by the old- 
fashioned arithmetic IF as follows: 


| |e CITYPE-2) 10, 20, 30 


(Rinally the output @using the data opposite ) would look like this: 





QYUL the contro! statements introduced in this chapter = viz: 


IF (logical expression) statement 

60 To label 

STOP 

Go To (label, ...,label), variable 

CONTINUE 

bo label control = initial, terminal, increment 
IF (arithmetic expression) label , label, label 
ASSIGN label TO variable 

Go TO yariable, (label, ...,label) 


ore executable statements and may therefore carry labels. 
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EXIERCIS'ES) CHAPTER 4 


4.4 A pair of simultaneous equations: 
ax 


+ bY = p 
ex + dy = 
may be solved using Cramer’s rule as follows: 
p Al ie p 
. J¢.4 _ le 4 
x = ab and y= art 
ke d ps d 


where the vertical bars denote determinants which may be 
evaluated as follows: 
& b 
cd 
Write a program to read the coefficients a,b,c,d and the 
right-hand side p,q then print the solution for X and Y. 
(Make sure you check the denominator is not zero before dividing: 
print the message “wo so.uvTion” if it is.) Then loop back fo read 
@ new right-hand side.’ 


= axd —- cxd 


4.2 Write Q program to generate and print the Fibonacci series: 


1 , 1, 2,3, 3, 8, 13, al, 34 y ++ 
where each successive term is the sum of the previous two terms. 
@ These numbers represent the numbers of rabbits generated by an 
initial pair of rabbits. The numbers grow rapidly because when 
each generation breeds so does its parents’ generation, its grand- 
parents’ generation, its great grand... Computer programs for 
printing Fibonacci numbers must be proliferating even faster, 
because this boring example seems to appear in many text books on 
programming.) Make sure your program stops before the term 
next to be printed exceeds the capacity of an integer variable, 
One way to do this is to check that Aalf the previous term, plus 
half the term before that, does not exceed half the capacity of a 
variable. 


4.3 \rite Q program { just as boring) to evaluate and print the 
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factorial of an integer read as data. As an example of a 
factorial, factorial five = written 5/ = is given by: 

SJ = 5x4 x3 *2*1 
Again make sure the capacity of an integer variable is not 
exceeded. 


A\\ternatively, make the program print a table of factorials 
Starting at 1/ 


2) 


TYPES OF ARRAY 
SUBSCRIPTS 
RIPPLE SORT (AN EXAMPLE ) 


EXERCISES 


WARES CLF LL = Fetcision, cocicar. comPeex 


Zn array is an arrangement of elements in one, two or three dimensions. 
The following: os Ks 

REAL A, BB, C¢ FONG 
DIMENSION  A(#) , BB(5,3) , C(3,2,3) 


specifies three arrays of type REAL depicted as follows: 





i, 2, 
(i 
ae Sa 
cca, [1 7-68 | 
1 2 CAYER 3) 
1 a 
ca 
eae meee 
i, 2, (LAYER 2) ) 
= 8 ane er, 
ce eee 
CG, 


[Zach array consisis of a number of “little boxes” called array 
elements. Except where specifically mentioned to the contrary, each 
array element may be used in the manner of a variable of like type. 
The element is written as the name of the array followed by one, two 
or three subscripts in brackets: 


5.38 
-6°-72 


7.68 
BB(4, 3) 





Phe DIMENSION statement has the form: 


DIMENSION array, array, + , array 


where: 
array is the symbolic name of an array followed 
by ifs dimensions in parentheses. For each 
array there may be one, two or three 
dimensions separated by commas. Except in 
the case of subprograms (Chapter 7) each 
dimension may conist only of digits. 


Whe DIMENSION statement — being a specification statement = is non- 
executable and therefore should not be labelled. D/MENS/ION statements 
should follow the type stataments as defined by the table on pagei7. 


Whe DIMENSION statement may be omitted if arrays have their dimen- 
sions specified by type statements: 


REAL A(7#), BB(5,3) , C(3,2,3) & Dimensrow 
INTEGER IA(3,2) <A omitred 


Zs in the case of variables, if the type of an array is not specified it 
becomes implicitly specified by the initial letter of the array’s name. 


I,J,K,L,M,N 
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A\rrays whose names begin with letters Z,J/,X,4,M,N are implicitly 
specified as type INTEGER: thosewith other initials as type REAL. So the 
previous two statements could be reduced to one: 


| | formenszan A(7), BB(5,3), C(3,2,3), IA(3,2) 


Where is a third way of specifying the dimensions of oan array = by 
the COMMON statement as illustrated below: 













INTEGER IA 
REAL A, BB, C 
OMMON A(#), BB(5;,3) , C(3,2,3), IA(3,2) 


where the use of COMMON is explained in Chapter 8. 









Dt is wrong to specify dimensions more than once for any array: 


REAL RED(6, 6) aI 
IMENSI@N  RED(6,6) 


“Arrays are stored by columns. For example, in the integer array: 


elements are stored in the order: 
JA(1,1), IA(2,1), IA(3,1), IA(1,2), IA(2,2), IA(3,2) 


ee 
ne THE FIRST SUBSCRIPT 
AG eT y.WARIES FASTEST 
LEE , 











_in other words: 


1) 2) 





SASS 


Gn many Fortrans no check is made on the bounds of an array, 
thus if is possible, for example, to write IxX=IA(4,1) to pick up fhe 
value stored in IA(1,2). This is allowable in Fortran 66 but not in 
Fortran 77. It may even be possible , using some compilers, to pick 
up the same value by writing IX= JA(4). for portable programs, 
however, such tricks should never be employed. The number and 
range of every subscript should conform to the dimensionality and 
size ranges originally specified. 


ZAlthough some Fortrans set all elements of each array to zero at 
start of execution most do not. The content of every array element 
= just like the content of every variable « is undefined at start 
of execution @ there is no saying what it might contain } unless pre- 
set by a DATA statement. DATA statements are described in Chapter 9. 


Ae 
LN 


CLYaLY 
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ONLY SEVEN FORMS ARE PERMITTED 
Cinteger constent « integer variable t integer constant ) 


[Zach subscript of an array element may be written as an_ integer 
constant as already illustrated: 


A(2) = 5.38 
BB(4,3) = ~-6.72 





Maat) 


A subscript may also be written @within limits ) as an expression 
of type JNTEGER : 


<I (WTEGER CONSTANT) 
<a (reser variaiet ) 


DS orecer VARIABLE’ # INTEGER CONSTANT ) 
<EWINTEGER CONSTANT ¥ INTEGER VARIABLE * ) 


Sewresee CONSTANT ¥* INTEGER VAR(ABLE' t 
INTEGER CONSTANT ) 


tyre an ARRAY ELEMENT 


AA ithough many Fortrans € including Fortran 77 ) permit more complicated 
expressions as subscripts there are only seven forms permitted by 
Fortran 66 « and all seven forms are illustrated above. Even an 
array element written A(1+#2*k) 93 or A(kK¥*2)>™& is non-standard 
and likely to prevent a program being fully portable. 





A(2*K -1) 


“An array element may be used in the manner of a variable except 
Where specifically noted. Thus it is wrong € page 40) {fo write: 


DO 10 ZZ = ZA(1,1) , TIA(1, 2) 





but the intended effect may easily be obtained at the cost of extra 
assignments: 


WJ = TA(1, 1) 


K = IA (1,2) 
10 I=J,K 





[lere are a few examples of the use of array elements. First oa 
piece of program to read values into a one-dimensional array 

@ often called a vector ». Each value is punched in the first ten 
columns of a card: 


REAL A(7) 

DO 10 T= 1, 8 EO Rate a ne 
READ I. SINGLE STATEMENT: 
pre ati nen ly READ(S, 100) A 


PAGE 9% 
FORMAT (F10.0) . Uy 





and here is a piece of program to “clear” a two-dimensional array: 


REAL BB(5, 3) 
De 20 ICO6LUM = 1, 3 
DO 20 TR@W = 1,5 


BB (IR@W, ICOLUM) = 0-0 
CONTINUE 





SO 


[nlere is Q piece of program to print a one-dimensional array @€ a 
vector ) as a column of values down the left-hand side of the output 


page: 


PO SO Fe hE? ES THIS LOOP MAY BE 
WRITE (6, 200) A(ZI) REPLACED BY THE 


SINGLE STATEMENT: 
ONTINUE WRITE (6) 200)A 


FORMAT (1X, F10.2) a mce 96D 





Qn the next piece of program oa vector is scanned to find the location of 
the first negative value ¢ if any ): 


DS 40 I= 1,7 
IF (AC) .LT. 0.0) 6@ T™] 50 
CONTINUE 


WRITE (6, 300) I 


FORMAT (1X, 13H LOCATIGN IS , £2) 




















Whe following example shows a matrix multiplication. Each element of 
array C(,) is made to contain the sum of the products obtained when 
each row of array AC(,) is multiplied by the corresponding column of 


arroy B(,). €.9. CONTENTS OF C(3,2) 1S THE INNER 
PRODUCT OF ROW OF ACs) AND 








2) 3) 


A(t, 8(1, 
SS & EE] 


AG, THE NUMBER OF COLUMNS OF Cae 
A(4, Muar aS THE NUMBER OF 
OF Cae ) 
[peg A(10,10), B(10,10) , C(410,10) 
eee ROWSA , COLSA, COLS8 


ROMA *'4 Se] SET UP THE 


COLSA = 2 DIMENSIONS 
COLSB =3 


Dd 10 J= 1, COLSB <1 
Dd 10 IL 1, ROWSA 
C(I,J) = 0.0 


CONTINUE 


DO 20 K=1, COLSB 
Dé 20 J = f, ROWSA 
Dé 20 I = {, C@LSA 
C(J,K) = CCT, K) + A(T, I) * BCI, K) 
CONTINUE 





he above piece of program could be made more elegant and efficient by 
omitting the “DO 10” loops altogether and replacing the three “DO 20”loops with: 


De 20 K= 1, COLSB 


DO 20 J= 1, ROWSA 
X= 000 3 Si itis) 


Dé 30 I= 1, COLSA 
X = X+#A(J,L) * BCI, k) 
C(I, K) = X 
CONTINUE 





but ts it as clear to you ? 
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[3 AN EXAMPLE TO /LLUSTRATE 
GPP S SUBSCRIPTED VAR/ABLES 

Gorting numbers into ascending order is simple in concept but remark- 
obly difficult to organize when there are large volumes of data. The 
example below uses one of the simplest techniques of all = the ripple 


sort =~ which is adequate for small volumes of data € a hundred or so 
numbers 9 stored as a vector. 


pe oe “Array A() = a column vector = is to be 

A(3) [10-2 sorted into ascending order; heaviest 

ACA) numbers sinking to the bottom. You can 

A(5) (3-5 reverse this order by reversing the condition 
of the logical IF statement. 


Start with an “index” I pointing to row 1; then advance J row by row. At 
every advance look oat the number JZ is pointing to « and also at the 

number one row ahead of ZI. If the former is greater than the latter 
swop the two numbers. 


[nlaving finished one “sweep” of I sweep again = but stop one row short 
of the previous sweep because the heaviest number must already have 
sunk to the bottom. 


@ontinve Sweeping = each sweep a row shorter than the previous one = 
until there is @ Whole sweep without a single swop in it or the length 
of sweep is reduced to nothing. 


(ere is the whole process: % shows where a swop has just occurred: 










HEAVIEST NUMBER 
HAS SUNK TO 
THE BOTTOM 


A(2) 
Ce 


Whe program opposite is designed to sort ao vector A(), having N rows. 
The vector is filled with numbers by the technique already illustrated; the 
results ore printed down the left margin of the ovtput page by the 
technique already illustrated. 


[Kogical variable EXIT is set trve at the start of each sweep but is changed 
to false if a swop has fo be made. Thus if there is a full sweep without any 
swops, control jumps to label 40 for an early exit. 


Oi 


2 


(lere is the full program: 


REAL A(100) , SWOP 
INTEGER I, KARDS, NSWEEP, LAST, LENGTH 
LOGICAL EXIT 


READ (5,100) KARDS 
IF ((KARDS LT. 2) .O0R. (KARDS .GT. 100)) STOP 1 


LAST = KARDS — 1 

Dé 20 NSWEEP. = 1, LAST <= 

EXIT = .TRUE. 

LENGTH = KARDS — NSWEEP <iF SHORTER LENGIN EACH SWEEP 
DO 30 I= 1, LENGTH 
IF ( A(I) .LE. A(It1)) G6 76 30 
SWOP = A(Z) 


A(1) = A(Z#4) pet 
A(I#1) = SWOP 
EXIT = .FALSE. 
CONTINUE 
IF (EXIT) GO T6 40 
® 50 I= 1, KARDS 
WRITE (6,300) A(1) oD, 
ROUNDED TO Two 


FORMAT (F10.0) a DECIMAL PLACES 


FORMAT (1X, F10.2) 








NSA NUMBER OF 
1 CARDS FOLLOWING 





SPhe output on unit 6 would look like this: 
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EXERCISES CHAPTER 5 


5.1 


5.2 


5.3 
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sing the early examples of this chapter as models write oa 
complete program for reading two arrays; performing a matrix 
multiplication; printing the result. Try the program on the 
following data: 


SALES PRODUCTS 
a) ae » 










SALES 
PERSON 





ArRay AC(,) 


NUMBERS OF ITEMS SOLD 





Array B(,) 


PRICES & COMMISSIONS LIST 





and the result should be the valves in the following table: 


SALES 
bond AMOUNT COMMISSION 
i Array CC) 


tG, 
C2,) 38.50 | 7.10 | = ArRAy A(,) x ARRAY B(,) 


c(3,]_ 30-00 _| 7. 


SUMS & COMMISSIONS EARNED 





Phere is Qa modification of the ripple sort called the shuttle sort 
which is explained in many books on programming. The method 
requires just one pass down the vector « but every time a 
swop is made, the number shifted up one place must be 
compared with the number now above it ™ and so on until this 
item has risen as far as it should go. Write a program to 
perform the shuttle sorf. 


“fhe enthusiastic reader may care fo explore three further 
methods of sorting explained in many text books on pro- 
gramming. These three methods are: 

@ the Shell sort 

® the monkey-puzzle sort 

® Quicksort 


of which the last two are mind-benders but worth the effort of 
Study. 


SV NAYS FONS ORS 


INTRINSIC FUNCTIONS 

BASIC EXTERNAL FUNCTIONS 
STATEMENT FUNCTIONS 

TRIANGLE @AN EXAMPLE ) 

ROUGH COMPARISON @AN EXAMPLE ) 
EXERCISES 


AN TRANS IC PUNE ULES “BUILT IN” TO FORTRAN 


The introductory example showed the line: 


BE ee ee 


where INT( ) is a reference to a function capable of converting its 
argument @in parentheses p from ao real value to an integer value by 
truncation. In other words all digits after the decimal point are thrown 
away. INT(3.999) yields the integer result 3 standing in place of 
INT (3.999) . 


‘Phere are thirty-one intrinsic functions in Fortran 66. The term “in- 
trinsic” implies these functions are usually built in @rather than linked ony 
fo the compiler. This, in turn, implies that you should not use any of 
their thirty-one names for any other purposes besides those defined oppo- 
site. In particular you shovld not use any of these names for functions 
of your own devising such as statement functions which are explained later. 


Phe form of a function reference (or invocation ) is: 


name (argument, argument, ... , argument ) 
where: 
name is one of the thirty one names defined in the 
table opposite. The type of each function 
(é.e. the type of value you get back from the 
function ) is indicated by the initial letter of 
the name of each function @ see initials below). 


argument is the name of a variable or array element, 
or is an expression of the type indicated by 
the initial letter of each symbolic argument 
in the table opposite : 


D signifies a function or argument 
of type DOUBLE PRECISION ; 

C signifies a function or argument 
of type COMPLEX ; 

L&M signify functions or arguments 
of type INTEGER ; 

e other initials signify functions 
of type REAL. 


A function reference may be used in the manner of o variable except 
where specifically forbidden. Thus if would be wrong € page 40) to 
wrife: 


| fog 40 T= INT(A4+B)*R INT (C¥D) 


but the intended effect could be had at the cost of extra assignments: 
le = INT(A+B) 


A = INTC(C#D) 
De 10 Law PFeik 





An argument of a function = whether intrinsic or otherwise «= may be 
an expression of much complexity: 


eae) DESCR = SORT (FLOAT (1B##2 ~4 #IA¥IC))/ FLOAT (241A) 
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UNSUION 


ABS (A) 
LABS (I) 
DABS ( D) 


AINT (A) 
INT (A) 
IDINT (D) 


AMOD(A1, A2) 
MOD (71, 12) 


AMAXo (11, 12, ..., In) 
AMAX1 (Af, A2,..., Aa) 
MAXO (J1, 2, ..., 10) 
MAX1 (Af, A2, ..., AM) 
DMAX1(D1, D2,..., D2) 


AMINO (11, 12, .-, £7) 
AMIN4 (AJ, AZ, -0, AM) 
MINO (J1,12,..., In) 

MIN1( Al, A2, ..., An) 
DMIN1i (D1, D2, .., DA) 


FLOAT (2) 
IFIX (A) 


SIGN ( A1, A2) 
ISIGN (71, I2) 
DSIGN (D1, D2) 


DIM (A1, A2) 
IDIM (11, J2) 


SNGL(D) 


DBL (A) 
REAL (c) 
AIMAG (c) 


CMPLX ( AR, AZ) 


CONJ (Cc) 


Dari wdon 


The absolute ¢ positive) value of an argument: 
ABS(-3-5) is 3.5 ; ABS(#+3-5) is 3.5 ; 
ZABS(-3) is 3 3; DABS(-3-5D0) is 3.5D0 


Vruncation: the sign of the argument is applied fo the 
largest integer Jess than the absolute value of that 
Orgument: AINT(-3.9) is -3.0 ; IN7(-3.9) is -3 3 IDINT(1.5D0) is 1 


Remainder when Al is divided by A2: = At —|A1+A2|xA2 
where the vertical bars contain an integer whose mag- 
nitude does not exceed A1l+A2 and whose sign is the 
same as the sign of Af+A2 : AMOD(S5.0,2.0)is 1.0; 
AMoD(-5.0, 2-0) is -1.0 @ similarly for Z2 and f2 ). 


Whe value of the largest argument: 
AMAX1 (-99.5, -16.4, 2-3) is 2.3 
AMAXO0(-99, -16, 2) is 2.0 @¢.e. REAL) 


There must be at least two arguments in each function. 


‘fhe valve of the smallest argument: 
AMIN1(-99-5, -16-4, 2.3) is -99.5 
AMINO(-99,-16,2) iS -99.0 (¢.e. REAL) 


Gonversion from type INTEGER to type REAL. 
Gonversion from type REAL to type INTEGER by 
truncation: the same effect as INT(A) above. 


Phe sign of the second argument @which should not be 
zero) applied to the absolute value of the first argument. 
SIGN (3.5, -1.2) is -3-5 5 SIGN(-3-5, -1-2) is -3.5 3 

€ similarly for r5szGN and DSIGN ). 


Yhe positive difference: DIM(A1,A2) is A1l- ABEND A2) 
and IDIM(I1,12) is I141-MINO(I1,12). 


Yhe most significant part of a double-precision 
argument expressed as a_ single-precision result 
= properly rounded, 


A single-precision argument extended to double precision. 


{Phe real and imaginary parts of a complex argument 
respectively: REALC(4.0, 2.0)) is 4-03 
AIMAG (( 4:0, 2.0)) is 2.0 


(xpress as a complex number with real part AR and 
imaginary part AZ: CMPLX(4-0, 2.0) is (4-0,2.0) 
Cin other words 4.0#2.0xW-1 «= a complex number ). 


Obtain the conjugate @ imaginary part reversed in sign) 


of a complex argument: Conw7( (4-0, 2.0)) is 
(4-0, -2-0) @in other words 4.0 -2.0xwW-1 ). 
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ABTS ECTERAAL GORSWIONS er rorrew 
@onsider this little program for printing a table of square roots: 


A = SQRT (FLOAT (I)) 
WRITE (6, 100) I, A 


FORMAT (1X, 14HSQUARE ROOT OF, 14, 3H IS, F10.4) 
END 





On the second line FLOAT() is an intrinsic function: SQRT( ) is a basic 
external function, Both kinds of function may be invoked in precisely the 
same way. The difference @ theoretically ) is that the programmer may 
override any basic external function by devising his own function @ see 
later 2 and giving it the same name. However, because some Fortrans 
do not stick to the letter of the standard it is safer mot to override 
basic external functions but treat them in the same way as intrinsic 
functions. Thus if you want to devise your own square-root function give 
it some name other than S@QRT : for example SQ@ROOT. 


Shere ore twenty-four basic external functions in Fortran 66 ond these 
are tabulated opposite. Many Fortrans, however, offer additional ones. 
For example there is no tangent function defined by Fortran 66 but many 
Fortrans offer 7AN() as a basic external function. Using such non- 
Standard functions invites problems of portability = don’t be tempted! 


Phe form of a function reference (oor invocation) is the same as that of 
an intrinsic function: 


name (argument, argument, ... ,argument ) 
where: 

mame is one of the twenty-four names defined in the 
fable opposite. The type of each function 
(i.e. the type of value you get back from the 
function ) is also tabulated because the initial 
letter of the function’s name does not always 
indicate type. 


argument is the name of a variable or array element, or 
is an expression of the type indicated by the 
initial letter of the symbolic arguments in the 
table opposite: 


A signifies an argument of type 
REAL ; 

D_ signifies an argument of type 
DOUBLE PRECISION 

Cc signifies an argument of type 
COMPLEX 


AAs in the case of intrinsic functions a function reference may be used in the manner 
of a variable except where specifically forbidden. Also the arguments may be com- 
plicated expressions, themselves containing function references. The following state- 
ment works out the fourth root of fhe absolute value of a REAL number stored in 
array element xX(2): 


eg = _S@RT(S@RT (ABS ( X(2)))) 
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FONSYION 


EXP(A) 
DEXP (D) 
CEXP(C) 


ALOG (A) 
DLOG(D) 
CLOG (Cc) 


ALOG 10 (A) 
DLOG10(D) 


SIN(A) 
DSIN (D) 
CSIN(C) 


COS(A) 
DCO$(2) 
c¢cos(c) 


ATAN(A) 
DATAN ( D) 


ATAN2(At1, A2) 


DATAN2( D1, D2) 


TANH (A) 


SQRT(A) 
DSQRT(D) 
CSQRT(C) 


DMOD( D1, D2) 


CABS(c) 


WS OF? 
FUNSWION 


REAL 


DOVBLE PRECISION 


COMPLEX 

REAL 

DOUBLE PRECISION 
COMPLEX 


REAL 


DOUBLE PRECISION 


REAL 
DOUBLE PRECISION 
COMPLEX 


REAL 


DOUBLE PRECISION 


COMPLEX 


REAL 
DOUBLE PRECISION 


REAL 
DOUBLE PRECISION 


REAL 


REAL 
DOUBLE PRECISION 
COMPLEX 


DOUBLE PRECISION 
(see also the 


DEF AMAT LON 


e raised to the power given by the argument: 
in other words the natural antilogarithm of the 
argument. £xP(0.0)is 1.0 ; EXP(1.0)is e. 


The natural ¢ base e) logarithm of an argument 
which must have a value greater than zero: 
ALOG(1.0) is 0.0 ; ALOG(e) is 1.0 


The common @ base 10) logarithm of an argument 
which must have a value greater than zero: 
ALOGy( 10.0) is 1.0 


The trigonometric sine of an argument in 


radians: * Nee | PL = 3-14169--. 
‘= a eal i Se SIN(-PI/6.0) is -0-5 
9 legos 42  \SIN(0-0) IS 0.0 
SIN (PI/2.0) is 1-0 


The trigonometric cosine of an argument expressed in 
radians: PI # 3.14159... 


= 1.0 
ea rn fOB a 7 COS(-PI/3.0) is 05 


ae ee €0S(0-0) is 1.0 


COS (PI/2.0) is 0.0 





The arctangent of an argument = the angle ¢ in 
radians ) whose tangent is... The range of the result 
is -1/2<angle¢ 1/2 


0 ol ae 


—1 
Laos 1/4 


| 
The arctangent of AL+A2 but the signs of Af and 
A2 are individually significant and A2 may be zero 
Qsimilarly for Di & D2). The range of the result 
is -W<angle¢ 1 


ATAN (ee) is %/2 
ATAN (0-0) is 0-0 
ATAN(-1-0) is -1/4, 


The hyperbolic tangent of an argument: 
CeA - e%)e(e4 + eA) = i.e. Simh(A)/cosh(A) » 


The square root of an argument which may 
not have a negative valve 


The remainder when Df is divided by D2: 
= Di —| Di +D2|xD2 where the vertical 


intrinsic functions bars contain an integer’ whose magnitude 
MOD(,) d& AMOD(,) does not exceed the magnitude of DL+D2 and 


REAL 
see also the 


intrinsic functions 


ABS( ), TABS() & 
DABS( ) ) 


Whose sign is the same as the sign of Di+D2. 


The REAL modulus of a complex argument. 
If ¢ is (AR,AZ) then CABS(Cc) is: 
SORT (AR ##2 + AI¥¥2) 


° 
t e=1+4,+ 


ae ie 
x 


| oe 


Ba ee 
+3 + 


Ey 


Led 


? 


{ radian (1% radians = 180°) 
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SAN UNSUIONS — reodeamnee” 


‘he following program is for computing the mass of a length of pipe, 
given its length, both diameters, and density of material: £,D2,D1, RHO. 


REAL £, D2, Di, RHO, A2, AL, W 
READ (5,100) E, D2, D1, RHO 


Al = 3.141593 * Di **2/4.0 y? 
A2 = 3.141593 * D2 #*2/ 4.0 

W = E * (A2—A1) * RHO 

WRITE (6,200) W ) 

STOP : ¢ 
FORMAT (4F 10-0) we 


FORMAT (1X, FHMASS IS, F10.2) 
END 










Whe lines beginning Af and AZ do much the same job and may be turned 
into a single statement function. The program could be recast thus: 


REAL E, D2, D1, RHO, W, X, AREA 
AREA (X) = 3-141593 * X *¥2/ 4-0 
READ (5,100) E, D2, D1, RHO 

W = E * ( AREA(D2) — AREA(D1) ) * RHO 


Ys. 
WRITE (6,200) W STATEMENT FUNCTION 
INVOKED TWICE 


[Rotice the statement function is defined immediately before the first 
executable statement. The statement function may be invoked @ some 
say “referenced” D) anywhere else within the program unit in which it is 
Gefined. But if one statement function invokes another the invoked one 
must be defined first: 


AREA(X) = 3.141593 * X ¥¥2/ 4-0 SH STATEMENT 

VOLUME (P,Q) = AREA(P) * Q ne 1K connec? 

WEIGHT (R,S,T) = VOLUME(R,S) ¥ T <eRORPE 

READ (5,100) E, D2,D1,RHO pg as 

W = WEIGHT (D2, E, RHD) -WEIGHT (D1,£,RHO) 

WRITE (6,200) W 

FUNCTION INVOKED 

TWICE 


Whe variable x in the first statement function is a dummy argument 
having no connection with any ofher X that might be used in the same 
program. We could just as well have defined the statement function: 

















‘DEFINGS A “STATEMENT FUNCTION’, 
CAULED AREAC ) 





<al 














































AREA (%) = 3.141593 * ri ** 2 / 4-0 


except that Fortran has no such character as £. The same goes for 
Gummy variables P/, 9,8, 5,7 in the piece of program above. 


ire = 3.141593 * X **2/4.0 
ay 
S DUMMY x! 
Se ARGUMENT * 
5 \ 












fhe form of definition of the statement function is shown below. This is 
not an executable statement and therefore should not be labelled. 


name (dummy, dummy, ... , dummy ) = expression 


where: 


name is the symbolic name chosen for the statement 
function. If shovid not be the name of an in- 
trinsic function, and it is safer that it should not be 
the name of a basic external function = or any 
of Fortran’s keywords. It should not be the 
name of any variable or array in the same pr- 
gram unit. The name may be declared ina type 
Statement, or the type of function becomes implicitly 
declared by the initial letter of name. 


dummy is the symbolic name chosen for @ dummy argv- 
ment. The types of dummy arguments may be 
specified in preceding ¢ype statements, otherwise 
their types are implicitly specified by initial 
letters @ I toN specify WTEGER; other initials REALY. 


expression is an arithmetic expression or logical expression 
involving all the dummy arguments. Array 
elements are not allowed in this expression. 
The type of expression must agree with the 
type declared explicitly or implicitly } for the 
function’s name. 


Phe statement function may be invoked in the manner of an intrinsic or 
basic external function as illustrated opposite. Here if is again: 


Foe W = E *(AREA(D2) — AREA(D1)) * RHO 
¥ re 


where the actval argument may be a complicated expression of appropriate 
type. But the actual argument may not be the name of an arroy, nor 
may it be an £x7ERNAL name € page 70). 





Fortran 6 has no tangent function although fortran 77 does. It is simple 
to provide a tangent function in the form of a statement function: 


and similarly for other trigonometric and hyperbolic functions missing 
from the definition of Fortran 66. For example: 


| [SZWH(X) = (EXP(X) —- EXP(-X)) / 2.0 


Fortran allows the expression in the definition to contain variables other 
than dummy variables, but this practice is not recommended because it 
may make a program difficulf to understand and errors hard to trace. 
In the example below A, 8,C and D are defined elsewhere in the program 
unit: their current values are used each time the statement function is 
invoked. A,8,C,D are called parameters. 


| | feuare(x) = ArxXee3 + BeXee2 + C#X + D 
t 


u t 
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AN EXAMPLE TO ILLUSTRATE INTRINSIC, BASIC 
EXTERNAL AND STATEMENT FUNCTIONS /N USE 
Given the lengths of two sides of a triangle and the included angle: 


ZN 
ere | 


Q 


a 


it is possible to derive expressions for the area, the other two angles, 
and the length of the third side. 


orea = 5 PQ sin 
g ; j - / Asinew 
iti sic on Vere 
-( Psinw 
A ange B = tan’ (> se— ) 


length of 
z third side, L = ,/P*+Q?- 2P@ cosm 
¢. 
Y 











[nlere is a set of data for the program. The values represent the two 
lengths, P and Q, followed by the included angle, “, expressed in 
degrees: 














REAL L 
DEGREE(R) = R¥ 180.0 [3.141593 <l 
RADIAN (D) = D*¥ 3.141593/ 180.0 <A 


STATEMENT FUNCTIONS 


degrees + fadians 
radians + degrees 













READ (5,100) P,Q, ALPHA 
C = RADIAN ( ALPHA) 

AREA =0-5% P *Q * SIN(C) 
A = INT (DEGREE (ATAN (Q ® SIN(C) / (P-Q* COS(C)) ))) 
B = INT (DEGREE (ATAN (P * SIN(C) / (Q@-P*COS(C))))) 
L = SQRTCP¥P + QQ - 2.0XP*Q*COS(C)) 





















yES! iT coup 
BE MADE 

SIMPLER, BUT 
THIS PROGRAM 
WAS CONTRIVED 
TO ILLUSTRATE 
FUNCTIONS 













WRITE (6,200) AREA, L 
WRITE (6,300) IA, IB 
STOP 

FORMAT ( 3F 10.0) 

@RMAT (1X, FHAREA IS, F10.2,17H OPPOSITE SIDE 15, F10.2) 
FORMAT (1X, 11HBASE ANGLES , I6, 4H AND, I6 ) 













and the output would appear as: 


AREA I5 
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ROUGI GORMRARISON sareneisdncrron™ 


SPests for approximate equality are rather messy in the body of a 
program. A _ statement function of type 40G/cAL such as: 


LOGICAL EQUAL 
EQUAL(X, Y, APPROX) = ABS(X-Y) .LE. APPROX 


can make such a program tidier. This statement could be invoked 
as follows: 


RUFFLY = 0.0014 
IF (.NOT. EQUAL (B*B, 4.O*A#C, RUFFLY))GO TP 10 


A quantity such as RUFFLY in the example above is better set ina 
DATA statement as defined in Chapter 9 but illustrated below: 


DATA RUFFLY [| 0.001 / 


where it may be easily found by the programmer and its value 
changed if the criterion for convergence is to be altered. 


“Phe logical function could, of course, be invoked with the criterion 
for approximate equality written as a constant: 


[] ZF ( EQUAL( P,Q, 0-001) ) STOP 





EXERCISES: CHAPTER 6 


6.4 Write a program fo compute and print the highest common 


factor of two integers read as data. Use Euclid’s algorithm: 


divide the larger integer, L, by the smaller, S 
if there is no remainder then S is the answer 
otherwise copy the content of S into L and 
put the remainder into S 

@ go back fo fhe first step. 


This method requires the use of the intrinsic function MOD(,) 


@.% Whe sum of the first integers is given by: 


N(N#1) 
2 





Make a_ statement function called ISuM(M) to deliver this result, 


6.3 ‘Whe area of o triangle with sides of length A,8,C is given 


by S(S-A)(S-BXS-C) where S is the semi-perimeter @ half of 
A+B+C ). Make a statement function called AREA(,,) using 
the three lengths as dummy arguments. @ As a statement 
function this will be long and messy. The function would be 
nicer wriften as a function subprogram. These are explained 
in the next chapter.) 


oA [Present the body of the “loans” program on page 31 as a 
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statement function. 


7/ 


AND 
BOUIN 
SOB ARAL 


FUNCTION SUBPROGRAMS 
SUBROUTINE SUBPROGRAMS 
EXTERNAL 

HORRORS 

AREAS OF POLYGONS (AN EXAMPLE D 
EXERC/SES 


PONSUION SUB RROERAWS arte Proseanmen 


he programming language called BAS/C provides a useful function called 
SGN( ) which returns a value of #1.0 or 0.0 or -1.0 depending upon 
whether the single argument proves to be positive or zero or negative 
respectively. A similar function may be included in o Fortran program by 
writing a@ function suvbprogram as follows: 


REAL FUNCTION SGN(X) 
IF (X .GT. 0.0) SGN = 1.0 A 
IF (X .£@. 0:0) SGN =0.0 espe deepen 


IF (X .LT. 0.0) SGN =-1.0 
RETURN 
END 





he function may then be invoked from another program unit as though it 
were an intrinsic function; actual arguments conforming in type and 
usage with corresponding dummy arguments: 





A function subprogram begins with a heading containing the word FUNCTION 
and ends with an £WD line to tell the Fortran compiler there are no more 
Statements to compile. The form of the heading is: 


FUNCTION name ( dummy, dummy, ...,dummy ) 
or: 

type FUNCTION name ( dummy, dummy, ..., dummy ) 
where: 


type may be INTEGER, REAL, DOUBLE PRECISION, COMPLEX, LOGICAL 
depending upon the type of result to be handed back to the 
program which invoked the function. If omitted @as in the 
first of the two forms defined above ) the type is implicitly 
declared by the initial Jetter of mame by the standard con- 
vention: I toN imply type /W7EGER ; other initials type REAL. 


name is the symbolic name given to the function by the programmer. 
This name must not appear in any ¢fype statement within 
the function subprogram; the means of declaring type has just 
been described. 


dummy is a symbolic name given to a dummy argument of any type. 
There must be at least one such argument in a function sub- 
program. Its type may be specified by a type statement in 
the subprogram or be implicitly declared by the initial letter 
of the dummy’s name. 


The dummy argument may represent oa variable oran array or 
an external subprogram @ all three forms are illustrated; the third 


on page 70). 
Names of dummy arguments may not be included in £@U/VALENCE 


Statements @ page 80) or in COMMON statements (page 7%) or in 
DATA statements @ page &) within the function subprogram. 


Somewhere in the subprogram there must be at least one executable state- 
ment of the form: 
RETURN 
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Whe RETURN Statement makes control return to the program which invoked 
the function. Also there must be at least one assignment having the 
name of the function on the left of the equals sign. This assignment 
determines what value is handed back to the program which invoked the 
function, 


he purpose of a function subprogram is to return a single valve in its 
place « just like an intrinsic function. The mechanics of Fortran 66 allow 
the programmer to write a function subprogram that: 

@ alters the values of its arguments 

@ changes values in COMMON @ page 76) 
but for the sake of portability these things should never be done. This is the 
province of the subroutine subprogram described next. 


[Mere is a different version of the SGN() function; it is of type /swrEGER but 
has an argument of type DOUBLE PRECISION : 


INTEGER FUNCTION IDSGN(X) se aaa 
DOUBLE PRECISION x —— AS A DOUBLE - 

seg i VARIABLE 
IF (X .67. 0.000) IDSGN = 1 


F (X .LT 0.0D0) IDSGN = -1 
RETURN 


REAL FUNCTION SUMSIX (VECTOR) 


DIMENSION  VECTOR(6) VN DUMMY ARGUMENT 
AS AN ARRAY 
SUMSIX = 0.0 Mame Cem 
I= 1,6 
SUMSIX + VECT@R(1) 


| \ozmenston x(6), Y(6) 
LES a! 





or it could be invoked in a more subtle way by using the name of an array 
element as the actual argument rather than just the name of an array: 





IAlere the actual argument specifies the first 
of six elements to be summed. Notice & 
that an actual argument of P(1) has 

the same effect as an actual argument 

of P, but that this is not true of Q(5)&Q. 


Ot is the responsibility of the person who 
invokes such a function to ensure that the © 
function works inside the bounds of the 
nominated array. Arrays with adjustable 
dimensions €@ these are permitted in functions ) 

are described on page 69. 


‘ 
=> 
2 
g 
3 
z 
> 
“a 
‘) 
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OOO UN = SUBS wvaron ane *cauen” 


frlere is an example of a subroutine designed to do the same job as that 
done by the Son() function illustrated earlier: 


SUBROUTINE SGNS (RESULT, EXPRSN) 
RESULT = 0.0 
IF (EXPRSN .LT. 0.0) RESULT = -1.0 


IF (EXPRSN .GT. 0.0) RESULT = +1.0 
RETURN 





this could be invoked (in the case of a subroutine it is usual to say “called”) 
from another program unit as follows. The point of return is the first execut- 
able statement following the cali; in this case at the line labelled 30: 


CALL SGNS(ANSWER, B*C) 
WRITE (6, 100) ANSWER <IP\ SOMOS SOK) 


A Subroutine subprogram has the same essential structure as a function 
Subprogram. The form of the heading is: 


SUBROUTINE same 
or: 

SUBROUTINE name (dummy, dummy, ... , dummy) 
where: 


name is the symbolic name given to the subroutine by the 
programmer. The initial letter has no significance. 


dummy is a symbolic name given to a dummy argument of any 
type. The type may be declared by a type statement 
inside the subprogram or be implicitly declared by the 
initial letter of the dummy’s name. Unlike a function 
subprogram a subroutine subprogram need not have 
any argument hence the first of the two forms above. 


A dummy argument may represent a variable, an array or 
an external subprogram (all three forms are illustrated; the 
third on page Top. 


Names of dummy arguments may not be included in £QU/VALENCE 
Statements @ page 80) or in COMMON statements @page 76) or in DATA 
statements ( page 8) within the subroutine subprogram. 


Somewhere in the subroutine subprogram there must be at least one executable 
statement of the form: 


RETURN 


to make control return to the program unit from which the subroutine was 
called. The point of return is the first executable statement after the CALL 
Statement as illustrated above. 


SPhe CALL statement is an _ executable statement of the form: 


CALL name 
or: 
CALL name Cactuval, actual, ..., actual ) 
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Where: 


name is the symbolic name of the subroutine subprogram 
to be called 


actual is an actual argument conforming in type and usage 
with the corresponding dummy argument of the 
subroutine subprogram being called. 


“Whe subrovtine subprogram may communicate via its arguments as 
already illustrated. But a subroutine may also communicate by referring 
to named common blocks and blank common. There are illustrations of 
this kind of communication in Chapter 12. 


\When Subroutine SGNS(,) @opposite) is called, the first actual argument can- 
not, of course, be a constant or expression; it should be the same of 
a variable or the name of an array element = in other words a “little 
box” in which fo carry back the result. Such an argument we can call 
an output argument because it delivers output from the subroutine. Conversely 
the second actual argument we can call an input argument, and this may 
be a constant or expression as well as the name of a variable or array 


element: 


Whe following subroutine is designed to add a given value to successive elements 
of a vector. The number of successive elements is specified as one of the 
Orguments. Notice how this number, as a dummy argument, is used as 0 


dimension of a vector. This usage is called ADGA DIARIES: 


SUBROUTINE ADDIN (VECTOR, NUMBER, VALUE) 


DIMENSION —_ VECTOR (NUMBER ) a 
DO 10 I= 1, NUMBER \=4 ‘pimensions 


VECTOR (I) = VECTOR(I) + VALUE 
















/wae 
CALL SGNS(ANSWER, 2.5 *A¥*B) 
CALL SGNS(VEC(2¥I+J), 2.5*A¥B) 












his subroutine may be called using the name of an array as the first actual 
argument : 










MAIN 
PROGRAM 






| «i cALe ©=ADDING(Y, 4, —3.0) 


or may be called using the name of an 
array element as the first argument: 


| | [cane ADDIN(V(3), 7, 4-0) ¢ 


and it is up to the person who writes the program that 
calls the subrovtine to ensure the subroutine does 
not work outside the bounds of the nominated array. 





Qn the above subroutine the first argument is an output argument and must 
therefore be a fame. The other two arguments are input arguments and may 
therefore be constants @as in the examples ) or expressions of appropriate type. 


“Phe use and abuse of arguments is further discussed on page 7. 
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SUBPROGRAMS WHOSE NAMES ARE USED AS ARGUMENTS OF 
OTHER SUBPROGRAMS ( SKIP THIS ON FIRST READING ) 


(Were are two function subprograms designed to return the cosecant and 
secant of an angle expressed in radians: 


REAL FUNCTIQ@N C@SEC(X) 
IF (SIN(X) .EQ. 0.0) STOP 1 
COSEC = 1.0/ SIN(X) 
RETURN 
END 























REAL FUNCTIQ@N SEC (X) 

F (COS(X) .E0. 0.0) STOP 2 
SEC = 1.0/ COS(X) 
RETURN 
END 















BASIC 
EXTERNAL 
FUNCTION 





BASIC EXTERNAL 
FUNCTION 


a 
a 
~ 


And here is a function subprogram with one dummy argument 
representing a function subprogram and the other representing an 
angle expressed in degrees: 


FUNCTION 
REAL ANGLE 
TRIG = RATIO (ANGLE * 3.141593 / 180.0) 











TRIG (RATIO , ANGLE) DUMMY 
ARGUMENT 
AS THE NAME 


OF A FUNCTION 







Phe TRIG(,) function could be invoked from another program unit as 
illustrated below. The EXTERNAL statement tells the fortran compiler that 
cos—C and S&C are not names of variables but ames of external sub- 


programs . 
ANGLE A 
iN DEGREES Ay 












|| JexreRnac 


eee = C * TRIG(COSEC, OF x ae (35: ik 
aah 
AS NAMES OF 
STERNAL FUNCTIONS 
AA fundamental principle of Fortran is that any program unit may be 
separately compiled. If you consider the two invocations of 7RIG above 
being compiled = with the Fortran compiler knowing nothing about sub- 
programs called Coséc or S&C = you will see the need for the preceding 


EXTERNAL statement. How else could the compiler know that Coséc and Séc 
were not just names of variables? 


COSEC , SEC 









“Phe EXTERNAL Statement has the form: 
EXTERNAL rame, rame, ...,name 


where: 
name is the name of an external subprogram @function or sub- 
routine } used in the current program unit as an argument 
of another function or subroutine. 


[Names of statement functions may not be declared FXTERNAL or be passed as 
arguments. The same applies to intrinsic functions. But although Fortran 66 
allows basic external functions to be declared £X7ERWAL = and their names to be 
used as arguments like CoSEC and SEC above * it is safer not to do so. Not all 
Fortrans agree what functions are intrinsic and what external. 


Whe EXTERNAL Statement is not an executable statement so should not be labelled. 
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HORRORS THE USE AND ABUSE OF LOCAL VARIABLES AND 
; ARGUMENTS @ SKIP THIS ON FIRST READING ) 


Qn Fortran 66 all local variables in a subprogram ¢ those not in 
COMMON storage s* page 76) become undefined as soon as the RETURN 
statement is obeyed. In other words a subprogram cannot remember 
what its local variables and local arrays contained from one invo- 
cafion to the next. In the following example: 


FUNCTION SUMATE (X) 


IF (X .EQ. 0.0) TOTAL =0.0 <kQumirauzarion) 


TOTAL = TOTAL + X 


‘TOTAL’ IS A 
SUMATE = TOTAL wy LOCAL VAR/ABLE 
RETURN 
END 


FORGOTTEN UPON, 
RETURN 





Successive invocations such as: 
|| A= Sumare(o.0) + SUMATE (50.0) + SUMATE( 25.0) 
Viraize) 
might well mot cause a value of 0.0+ 50.0 + 25.0 = 75.0 to be assigned 
to variable A. In many Fortrans these values would be retained, but 
reliance on such a feature makes oa program non-portable. 


A subprogram cannot be expected to return values via its arguments if 
one actual argument becomes associated with another. Consider: 








FUNCTION NORTY (J, kK) 


J=K#i1 4g 
RETURN 


END 








which is perfectly reasonable if invoked using two distinct actual arguments: 


LAST =2 
NOW = NORTY (NEXT, LAST) 


which would result in NOW becoming 2 and WEXT becoming 3. But if 
NORTY were invoked as follows: 


NO =2 en 
IDOUBT = NORTY (NO, NO) 


what would be the value assigned to I00UBT = 2 or 3 ? And what would 
be contained in variable No ? Different fortrans would probably deliver 
different values. 


QOhnintentional associations cause obscure bugs in Fortran programs. Be 
particularly careful about associating a dummy argument with an entity 
in COMMON storage € page 78). 


“AA subroutine subprogram @not a function subprogram) may, according 
to Fortran 66, have an actual argument which is a Hollerith constant 
@e.g. 4YCARD ) corresponding to a dummy input argument. However, the 
use of this feature is bound to make a program non-portable because 
the number of characters that can be stored in a variable differs 
from one computer to another and from one Fortran to another. 


CALL THINGY (IZ) 
CALL THINGY ( 4HCARD) 
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OF: ILLUSTRATING A FUNCTION SUB- 
PROGRAM AND ITS INVOCATION 


@onsider the diagram on the 
right. It is a simple matter 
to showt that the speckled 
area is given by: 


Aj = 4(NY -NN) eS 
= £ (2x3 - 2.5x1) 
= 1.75 





“Whe same formula may be used 
for computing the area on the 
left. But this area would turn 
out to be negative: 

NH Ai = 3G - XN) 

£ (32.5 - 5x4) 





= -6.25 
X4o%% 
Y / 
“Phe formula may be applied to E en 3% 


Sequential sides of a polygon = 
and the areas of triangles 
summed to give the area shown 
on the righf. 





ut if the polygon is closed, as 
shown on the left, the sum of all 
the triangular areas will equal 

the area of the polygon. 


SPhe bounded surface must be 
kept to the deft of each arrow. 
No check is made on silly data 
Where one edge crosses another 
as in a figure of eight. All points 
should have positive coordinates. 





t 
the area _ enclosed between the two vectors 
A ond B is given by $(A x B), where 











P FF F 
AxB = Q,Q@,0) = (a, bz - by an) F vs 
b; 62 0 ; 
unit t 
vectors ff 
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[Alere is a@ function subprogram devised to return the area of a poly- 
gon of W sides, where the coordinates of sequential vertices are 
Stored in vectors X(€) and Y(): 


REAL FUNCTION AREA(X, Y, WN) 


regi cso x(W) , Y(N) sw Diaicaratt pmtvavon ) 
D@ 10 L= 1,N CLOSES THE POLYGON BY 

J a LL JOINING FINAL VERTEX TO 
If (rhe: Fat Se eee 


AREA = AREA + 0-5 # (X(1)#Y(I) - X(T) #Y(1) ) 












DIMENSION P(50) , Q(50) 

READ (5,100) NUMBER 

IF ((NUMBER .GT. 50) .OR. (NUMBER .LT. 3)) STOP 1 
DO 10 I= 1, NUMBER 


(5, 200) P(r), Q(z) 


C@NTINVE 


A = AREA(P, Q, NUMBER) 


WRITE (6, 300) A 





WMotice that variable I and label 10 in the subprogram have nothing 
whatever to do with variable I and label £0 in the main program. 
Communication is via function name and arguments only, 


A set of data might be: 





A\nd the result would look like this: 
© 


© AREA IS 
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EXERCISES’ CHAPTER 7 


We Phe “Macaulay function” AMAC(X) returns the valve represented 
by its argument, X, if this value turns out to be positive; 
Otherwise the function returns a@ value of zero. Write AMAC() 
as a function subprogram. 


7/7. © \edork through exercises 6.2, 6.3, 6-4 but this time write function 
subprograms instead of statement functions. 


7.3 Recast the ripple-sort program on page 52 as o subroutine sub- 
program RPSORT(A,N) in which vector AC), with adjustable 
dimension N, is sorted. Write a main program to call this 
subroutine. @ See footnote.) 


Vo Ch Rework exercise 5.1 so that the matrix multiplication is by a 
call fo subroutine MATMUL(A,B,C,I,J,K) in which array A(I,J) 
is multiplied by array 8(J,k) to give array C(ZI,k). Write 
this subroutine and a main program to test it. @See footnote.) 


7.5 Write a “library” of subroutines to perform the fundamental 
operations of matrix arithmetic: 

copying from one array to another: MATCOP(A,8,I,J) 

addition and subtraction: MATADD(ISIGN, A, 8,C,1I,J) 

scalar multiplication: MATSCL(FACTOR, A, B,I,J) 

transposition: MATTRA(A,8,I,J) 

clearing a matrix: MATZER(A, I,J) 

creating an identity matrix: MATIDN(A,TI) 

matrix multiplication @ as exercise 7.4 above ) 

matrix inversion: MATINV(A,8,I) 


For an introduction to computing with matrices see: 
“Illustrating BASIC”, Donald Alcock, Cambridge University Press, 1977. 


Subroutine MATIDN(A,1) is done for you below, but you may want 
to devise a more elegant solution: 










SUBROUTINE MATIDN (A, I) 
DIMENSION A(I,I) 
DO 10 K=1,1L 





= ‘REPLACE ZERO 
DE IG SES sh DIAGONAL ELEMENT 
A(J,K) = 0-0 SPA WITH 1-0 
ZF (TS £0. K) A(T, K)= 1.0 1) 2) 3) 





CONTINUE A(t, 
AQ, 
AG, 


footnote: you may wish to add an argument of type 
LOGICAL which is set {true if the subroutine 
succeeds, or false if the actual arguments 
specify an impossible situation such asa 
negative number of elements. 
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SVAN Sulwav als 


COMMON 
COMMON @ CONTINUED ) 
STACKS © AN EXAMPLE ) 
EQUIVALENCE 

CHAINS @AN EXAMPLE ) 
EXERCISES 


(0 VN V\ey™ | A MEANS OF COMMUNICATION BETWEEN PROGRAM 
UNITS VIA SHARED VARIABLES AND ARRAYS 


‘Phe COMMON statement declares names of variables and arrays that may 
be used by any of all program units in a complete program. Consider the 
following main program: 





program: 


SUBROUTINE NEGATE 
COMMON B, A, Z(3,2), L,M, N(3) 
N(1) = -N(1) 
N(2) = -N(2) 
-N (3) 





“Phe elements of common storage declared in both program units bear a 

one-to-one correspondence « despite the completely different names being 

used. What the main program calls J the subprogram calls M, but it is 
nevertheless the same storage location in the computer. Notice that what 
the main program calls A the subprogram calls B; and what the main 

program calls B the subprogram calls A. 


The effect of fhis subroutine subprogram when called is to negate the first 
three elements of array K() on behalf of the main program. 


‘Whe order of elements in common storage is the order declared in the 
COMMON statement: 





and because the subroutine above is only concerned with the array N) 
at the very end of common storage it would be allowable to simplify 
the common statement so as to “step over” unwanted elements as below: 


SUBROUTINE NEGATE 
OMMON R(8) , I(2) , N(3) 
R 
‘STEP OVER? R(3) 


EIGHT 
hacen € , 
STEP OVER 
R@) Ls TWO INTEGERS 


REALS 
Some pitfalls await the programmer who assumes anything about the 
relative sizes of elements having different type. These potential 
pitfalls ore discussed later. 
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Phe opposite page illustrates blank €or unlabelled ) common. There is 
also named (Cor labelled }) ComMon as depicted below: 


INTEGER ON 
REAL P(3) 
DOUBLE PRECISION (2) 
COMMON / ALPHA / P,N 

COMMON / BETA / Q 
























SHARED BETWEEN MAIN PROGRAM 
AND SUBPROGRAM ‘SRA’ 


THESE S sies eats BE DIFFERENT 
AS LONG AS THE ARRANGEMENT OF 

REALS de INTEGERS ($8 KEPT THE aX: 
SUBROUTINE SRB 
DOUBLE PRECISION Q(2) 
COMMON / BETA / Q 
















Q(2) 






DOUBLE PRECIS/ON VECTOR 





SHARED BETWEEN MAIN PROGRAM 
AND SUBPROGRAM ‘SRB’ 





“Any number of named blocks may be specified. Each is accessible fo 
any program unit that declares the block’s name correctly and _ speci- 
fies a block of édentical length. @ Blank common need not be made the 
same length in every program unit: different lengths are illustrated opposite.) 


TPhe form of the common statement is: 


a a 
BLANK 
COMMON name, name, ... ,7ame {Ein 
or: 


COMMON /block/ name, name, ... ,name 


where: 
block is the symbolic name * unique over all program units = 
used to name a common block. This block may be 
referred to by any other program unit that declares 
Q@ common block of the same name and size. 


Omission of a block name @ leaving a pair of slashes 
with nothing between them ) specifies blank COMMON. In 
such a case the slashes may be omitted also, thus 
achieving the first of the two forms defined above. 


name is the symbolic name of a variable of any type or the 
name of an array of any type. If it is the name of an 
array the name itself may be followed by the dimensions 
of that array in parentheses = provided that these 
dimensions are not also given in a type or DIMENSION statement. 


[Fortran 66 allows more than one block to be declared in a Single comMMON 
statement, but this can be confusing especially in the use of commas: 


eae | COMMON / ALPHA/ P,N /BETA/ Q EG 0.K. BUT CONFUSING 
Pee a 


Dummy arguments may not be declared common, nor may an array in 
common be given adjustable dimensions: 


SUBROUTINE WRONG (A, I, J) 
DIMENSION ACI, J) 


COMMON A 





(continued ) 


77 


COM'AOM (cosmmweD ) RULES FOR ENSURING PORTABILITY 


[Fortran 66 refers to storage units of which one is occupied by a 
variable or array element of type /WTEGER, REAL or LOGICAL and of which 
two are occupied by a variable or array element of type DOUBLE PRECISION 
or CompLEX. So ideally it should be simple to map one common block 
upon another; in other words to associate items one with another. 















ones oa 1a) ie Senate a) 
EAL io, ot 
LOGICAL L(4) ge. 45) But there are few 
oe le (1) Fortrans in which 
DOUBLE PRECISION D(2) aO*~-Nn_.._] sie dedsen a 
COMPLEX C2) LD to (2) is 1deal wou 
“4) apply. If would be 
ee ee ees Dit) £3) remarkable if 
Di2 Ka) complex elements 
REAL oi \ Va C() and k(2) found 
piace weaetie ae caf Kk)" themselves in the 
cae ay ° C(2 same pair of 


COMPLEX K(2) 
LOGICAL M,N 
COMMON /GAMMA/T,S, E,K,M,N 





storage units. 


STORAGE UNITS 
/N COMMON 
BLOCK GAMMA 


{\n reality the fundamental storage unit is the computer word which has 
a different length on different computers. It may be 8, 16, 24, 32, 
36, 48 or 60 binary digits @bits) in length. On a typical computer 
having a 16-bit word an integer or logical variable occupies one word, 
a real variable occupies two, a double precision or complex variable 
occupies four words. The diagram above would |ook quite different 
if drawn with @ word as a_ storage unit. 


Sb make things more complicated, some 
compilers “map” multi-word items onto words | eomenon / Decra / Py Re 


of store in such a way that “dead” words ae para 
are left in the common block « making pereerrref BEAD) 
association difficult if not impossible. ey A 


ZAnd in some fortrans, logical 
arrays are “packed” info successive 


counen / TETAS 4(3) (0) words using only one bit per element, 


whilst each logical variable occupies a 


7) whole word. It ts therefore impossible 
| eeeaee Pee Fn 3 oe to associate logical variables with 
— logical array elements. 





[But it zs possible to use blank and named common in portable programs. 
The rules are these: 


@ always declare common variables and arrays in the order: 
DOUBLE PRECISION then COMPLEX then REAL then INTEGER then LOGICAL 


@ never associate items of different type 

@ never associate logical variables with logical arrays 

@ fo be quite safe from losing data in named common blocks 
(because of overlay problems on some computers ) include the 
definition of named blocks in the main program if any of 


their initial values are to be changed during execution. 
. (See Bibliography «= second book « on this subtle point. ) 
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AN EXAMPLE 70 ILLUSTRATE THE USE OF A 
NAMED COMMON BLOCK 


Sitacks are widely used in programming. The technique has its own 
vocabulary involving the word “push” @fo mean adding a number to a 
stack ) and “pop” (to mean taking a number off a stack D. The technique 
= used in anger on oe 129.0 aN be — 8 follows: 


~ A, 
ae >. 22 f *  S. 


Phe information pushed and popped may include integers representing 
labels of statements, letters in @ word, arithmetic symbols « in fact 
information of almost any kind. The picture above illustrates a stack 
of real numbers. 


flere ore three subroutines for maintaining a stack of real numbers in 
Q@ common block named STACK. The first subroutine should be called 
to clear the deck, but this subroutine could be replaced by a BLOCK 
DATA subprogram page 87) to initialize the value of JPOINT. 


SUBROUTINE CLEAR 
COMMON / STACK / A(10), IPOINT 


IP@INT = O 
RETURN 





Whe following subroutine is for pushing a real number onto the stack: 


SUBROUTINE PUSH ( EXPRN) 


COMMON / STACK / STK(10) , IPT EN 


IF (IPT .GE. 10) STOP 10 <E* srack CAUSES 


IPT = IPT +1 ERROR STOP 
STK (IPT) = EXPRN 

RETURN 

ND 


And the following subroutine is for returning the value popped from the stack: 


SUBROUTINE POP( TOP, OK) 
LOGICAL OK 

COMMON / STACK { S(10) , IPNT 
OK = IPNT .GT. O 





EMPTY STACK 


IF (.NOT. OK) RETURN 
TOP = S(IPNT) RETURN FALSE. 





All three routines refer to a common block named S7ACK. Names of variables in 
each subroutine have been made different so as to emphasize the independence of 
names = and the total dependence on order « within a common block. 


Uypical invocations of these subroutines are: 


CALL CLEAR ee ee) 
CALL PUSH(A/B¥*#2) ha ge 
Pfc Ate parc yasu, Any) <A vatu’ Fon 10P oF s7HK 
TOP 


-NOT._ ANY 
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A MEANS OF SHARING STORAGE SPACE, 
BUT Or FOR MAKING SYNONYMS 





@onsider the following function subprogram contrived to return the value of 
y in the cubic equation y=ax3+bx*+cx+d given the value of x and 
the four constants: 


REAL FUNCTION CUBIC(X, A,B,C,D) 
R= AxX x¥3 

S=#*R+ BeXt¥2 

T=S + C¥X 


U=T+OD 
CUBIC = U 
RETURN 
END 





“Phe intermediate variables R,5,7,U may be made to share one and the 
same storage location by inserting an £QU/VALENCE statement in the function 
subprogram as follows: 


[pre FUNCTION CUBIC(X, A,B,C, D) 
QUIVALENCE (R,5S,T,U) 
Phe form of the £QU/VALENCE statement is: 


EQUIVALENCE (list), (list), ... , Clist) 


where: 
list is a list of two or more items separated by commas, Each 
item in the list may be the nome of a variable or array 
element, but not a dummy argument. In the case of an 
array element the subscripts must be integers, and there 
must be as many subscripts as the array has dimensions. 


Despite the apparent simplicity of the E£@U/VALENCE statement there are 
difficulties not always appreciated by Fortran programmers «= difficulties 
that can easily make @ program non-portable. Consider: 


REAL OD, &, F 
EQUIVALENCE (D,€) 





Because F and D are of identical type, in most Fortrans the statement F=D 
would cause O value of 5.0 to be assigned to F. But Fortran 66 says that 
D @dhence also F ) would be wnaefined. Assignment to a variable or array 
element undefines its equivalenced room-mates. In other words the £QU/VALENCE 
Statement does not imply synonyms. This one ghastly fact is enough to 
make this statement a danger in programs intended fo be portable. 


When two array elements are equivalenced the other elements in each array 
fall into place as depicted below. Two- and three-dimensional arrays are 
strung out by columns: 


REAL A(6), B(4,2) Acayf Bat, 1) 

EQUIVALENCE ( A(3), B(2,1) ) > Agi ]] B2, 1) <a 
AW FT] BG, 1) 
AGS) [FT] B41) 
A@)|F ~~ {] B(1,2) 





EXTENSION Sees > 
RESULTS — Hey 
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By the device just described it is permissible @but not good practice ) to 
extend common in a forward direction: 


A(6), B(4,2) SEE SKETCH ON 
COMMON A Agr OPPOSITE PAGE 


EQUIVALENCE ( A(3), B(2,1) ) 





but it is not admissible to extend common “backwards” @ you can’t create an 


A(o) or A(-1) D: 
REAL A(6), C(3) x 
COMMON A 


QUIVALENCE (A(2) , C(3)) 


and it is obviously silly to equivalence one 
element of an array with another. This 
would imply “folding” the array over itself. 
This mistake could be made indirectly; 

























[fortran 66 permits a single subscript of unity @no other ) to imply the 
first element of a two- or three-dimensional array. Fortran 77, however, 


does not allow any such device: 
D(1, 1) i 


DIMENSION A(6), D(2, 2) 
EQUIVALENCE (A(1), D(1)) 
D(2,1) 
D(,2 


a tas D(2,2 














You can achieve what was intended above with greater clarity , and without 
falling foul of Fortran 77, as follows: 


| | fequrvacence (A(t), DU1,1) ) 


or even by equivalencing, say, A(3) with D(1,2). 


(Fortran 66 allows items of different type to be equivalenced. For 
example the arrays R&(,) and I() may share space as a result of the 
following £QUIVALENCE statement: 





Assuming a single computer word is used to store an integer, and a 
pair of words to store a real, arrays R(,) and I() might occupy 
the same space. In one well-tried and much-used program the 
programmer had taken advantage of this relationship, as a result of 
which the program would run on computers in which integers are 
stored in Single words and reals in double words. On moving this 
program to a new computer = in which an integer occupied the same 
space as a real = the portability problem was overcome by the 
drastic means of turning all real arrays and variables into double- 
precision arrays and variables. To avoid problems of portability, 
therefore, do not equivalence items of different type without forethought. 
of> 
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UAVS ILLUSTRATING A PROGRAMMING TECHNIQUE 
CALLED LIST PROCESSING 


“Phe manipulation | of chains is called dist processing = a fundamental 
tool of programming. Chains are essential for solving the problem on 


page 126. 


She simplest form of chain is illustrated below. It has a head and 
a vector of dinks. Corresponding to each link is a row of data: 
2 





DATUM (2, 
DATUM(3, | 52-8 | 17-9 | 
DATUMG, 
DATUM (5, 
DATUM (6, 





ZERO MARKS END OF CHAIN 


‘Phe following piece of program would cause dota linked by sucha 
chain to be printed: 


tI 
IF (N -EQ. 0) STOP 


WRITE (6, 200) DATUM(N,1) , DATUM(N, 2) 


STANDARD SHAPE FOR 
N = LINK(N) <I CLIMBING DowN TO 
NEXT LINK OF CHAIN 


GO T 10 
FORMAT (1X, 2F10.2) 





Quppose there were now some data in row 5 of array DATUM(,). The. 
new data could be linked to the head of the chain by the following 
piece of program in which variable I is first set to the desired row 


number ): 
DATUM(T,1) = 56-78 


DATUM(I,2) 21-73 


LINK(I) = ITHEAD weg STANDARD SHAPE FOR 
= LINKING RoW I 
IWEAD = IT SFR 70 HEAD OF CHAIN 


and now the picture is as follows: 


LINK(1)[|-— | 





1) 2) 
DATUM (1, 
DATUM (2, 
DATUM(3, 
DATUM (4, 
DATUM (5, 
DATUM (6, 


IHEAD 










“Phe last row linked may be unlinked @discarded ) as shown below. 
@For the sake of illustration the row of data is printed before being lost. D 


WRITE (6, 200) DATUM(IKEAD, 1) , DATUM(TIHEAD, 2) 


STANDARD SHAPE FOR 
IHEAD = LINK (IHEAD) Teta tte (Oe ) 
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The output would be: 


© 56.78 24.73 © 
Oo wat t) 
and the picture of arrays would once again be the first picture opposite 


(the one with INEAD containing 2 ) but with “garbage” in row 5 of 
LINK( ) and of DATUM(;). 


Notice these pieces of program cause the dast row linked to be the 
first discarded = like pushing and popping a stack. Indeed programmers 
Often use this chaining mechanism to organize stacks. 


VWhen using chains in more complicated ways it is useful fo incorporate 
the head of each chain into the vector of links. This wastes a little 
space in the associated array of data, but there become fewer special 
cases to deal with oat the heads of chains. 


LINK (1) 2 <Gneno DATUM(4, 


LINK(2),-3 | DATUM (2, 
LINK (3) DATUM (3, 
=—-— ——— oo 





A\s an example, suppose we want a subroutine for discarding a 
certain link part way down a chain. The link to be discarded is 
the one linking a given value @ VALv) to be found somewhere in 
the first column of array DATUM(,). First of all this value has to 
be found by running down the chain and peeping one link ahead. 
If the value cannot be found then the subroutine returns a logical 
argument set false. 


SUBROUTINE DSCARD (VALU, Ox) 
LOGICAL OK 
COMMON / CHAIN / DATUM(100, 2), LINK(100) 






NOW = 1 


NEXT = LINK (NOW) 


@K = NEXT -NE. O 


IF ( DATUM (NEXT, 1) .£Q@. VALU) G@ TO 20 


NOW = LINK (NOW) SKycuma bows) 


GO 7%) 10 
Then if the value is found the link may be discarded. 
LINK (NOW) = LINK (NEXT) 


END 
If you were to write this with the head of the chain in JHEAD you would 


have to make a special test for an empty chain before starting the 
search for VALU. 


[hist processing is extremely useful, and, fo enthusiasts, addictive. There is 
Q good introduction in fortran Techniques by Day @ see Bibliography ). 
This explains how to deal with “garbage” such as that left behind by the 


Qbove subrovtine every time a link is discarded. 
a 
-*) 
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EXIERCISES CHAPTER 8 


SPs | A\iter the “jibrary” of subroutines specified in Exercise 7.5 so 
that each subroutine communicates via a named block of 
comMMON rather than having names of arrays as dummy 
arguments. 


83.2 Whis is a challenge. Write a program to perform the ripple 
sort on numbers in vector AC) and print them in sorted order. 
But do not disturb vector A(). Instead link the elements of 
A() by a chain in Z(): 


A(2) | 6-5 
A(3) {13-9 | 
A(4) [10-2] 
A(5) 
A(6)|_ 3-5 | 


—~— 





then move the links using the logic of the ripple sort so as to 


end up like this: 
“line 
UNDISTURBED 





8.3 Spurn your program (above) into a subroutine that will sort on 
any column of a two-dimensional array: 


| | cae = mysorT (A, I, J, KEY) 


where I and J specify the dimensions of array A(ZI,J7) and KEY 
specifies which of the columns is to be used as the key for 
sorting. Test the subroutine by sorting on each column of a 
three-column array of numbers. Print the complete array = in 
the order determined by the chain = after each sort. 


Hint: put the vector of links into a named common block. If 
you write a subroutine fo do the initial linking: 


Peace LIWKER (LENGTH) Lorine Moe) | 


you will only need one call to LINKER = before the first call to 
MYSORT. Subsequent calls to MysoRT would then begin with the 
chain left behind after the previous call to Mysorr. 
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AT ALR ZAT ION 


DATA 

BLOCK DATA 

CHARACTERS 

STATE TABLE @AN EXAMPLE D 
EXERCISES 


LAN A STATEMENT FOR INITIALIZING VARIABLES AND ARRAY ELEMENTS 


Phe DATA statement is for setting initial values into local variables before 
execution starts. This is not an executable statement so: 

@ the DATA statement shovid not be labelled 

@ variables cannot be reset by a DATA statement during execution. 


flere is Q subroutine designed to return the surface area and volume ofa 
sphere given its radius as an input argument: 













SUBROUTINE SPHERE (RADIUS, AREA, VOLUME) 
DATA PI / 3.141593 / 

AREA = 4.0% PI *® RADIUS *¥*2 

VOLUME = 4.0* PI*® RADIUS *%3/ 3.0 
RETURN 
END 






LOCAL VARIABLE ‘PI' 
ts INITIALIZED TO 
3.141593 BEFORE 
EXECUTION ~ ¢.@. 
AT ‘COMP/LATION TIME’ 














LAore than one value may be set by a DATA statement. Here is a function to 
return the number of days in a month = given the month and year: 


INTEGER FUNCTION DAYS(MONTH, YEAR) 

INTEGER MONTH, YEAR, M(12), FEB, LEAP 

DATA M(1), M(2), M(3),M(4), M(5), M(6), M(7), M(B), M(39), M(10) 
ff M&, 2%, RH, 30, SM, 30, BM, 9, 34/, 

M(11), M(12)/ 30, 31/, FEB, LEAP/ 2,4 / 

DAYS = M(M@NTH) 

ZF((MONTH .EQ. FEB) .AND. (MOD(YEAR/ LEAP) -£@. 0)) DAYS=29 

RETURN 

END 

















Whe form of the DATA statement is: 
DATA names/constants/, names /constants/, ... , names/constants / 
where: 
names is a list of names of variables or array elements or both. 
Items are separated by commas. Items may not be dummy 
arguments, nor may subscripts of array elements consist of 
anything but digits. 


constants is a list of constants €Hollerith constants allowed ) separated by 
commas. Arithmetic constants may be signed + or-. Any constant 
may be preceded bya “multiplier” of the form: 
count * 
where count consists only of digits and specifies the number of 
times fhe subsequent constant is implied. 


‘Shere must be a one-to-one correspondence between names in the sames list 

and constants in the constants list. Each item in the names list must agree in 
type @ INTEGER, REAL, LOGICAL etc.) with the corresponding item in the constants list. 
This agreement does not apply to Hollerith constants which may, according 

fo Fortran 66, be stored in variables or array elements of any type. This 
topic is covered later in detail. 


LAany Fortrans allow you to write just the name of an array instead of 
laboriously listing its elements in sequence. For example the three-line DATA 
statement above is acceptable as follows: 


[Dara M/ 31, 26, 31, 30,31, 30, 2#31, 30, 31, 30,31/, FEB, LEAP/2, 4/| 


but this is, nevertheless, mon-standard Fortran 66. 


+ some Fortrans exclude Hollerith constants here. 
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A SUBPROGRAM FOR INITIALIZING VAR/ABLES 
AND ARRAY ELEMENTS IN NAMED COMMON BLOCKS 


Phe illustrations opposite show the initialization of variables and 
array elements (focal to the program unit in which the DATA statements 
appear. Fortran doesn’t allow items in blank common to be initialized 
under any circumstances. Items in named common blocks , however, may 
be initialized = but only in a BLOCK DATA subprogram. 


BLOCK DATA A 
REAL R, RAR BLOCK DATA 
INTEGER 1, TAR SUBPROGRAM 
DOUBLE PRECISION OD, DAR 

LOGICAL L, LAR 

COMPLEX C, CAR 


DIMENSION RAR(100), TAR(100), DAR(100), LAR(100), CAR(100) 
COMMON /IGTA/ R, RAR, I, IAR KS ONE ITEM MAY NOT 
COMMON / KAPPA/ D, DAR, C,CAR, L, LAR <td BE 'N MORE THAN 


OWE BLOCK 
EQUIVALENCE (IAR(1), J), (IAR(2), K), (TAR(3),M) 


DATA J,k,M/0,1,0/, RAR(1),RAR(2) /100.0, 0-0f 
DATA D,L,C/1.5D0, .TRUE, , (2-7, 3.6) / 
END 





A BLOCK DATA subprogram begins with the statement: 

BLOCK DATA 
and ends with an ewno line: 

END 
between which can come any number of fype statements @ REAL, INTEGER, 
DOUBLE PRECISION, LOGICAL, COMPLEX ) and any number of D/MENS/OoN statements, 
common Statements, £Qu/vALENCE statements and DATA statements. This order 
should be preserved; see page 17. No other kind of statement is allowed: 
the example above illustrates a// these non-executable statements. 


Fortran 11 (and several other Fortrans ) permit any number of 810Ck DATA 
subprograms = each having a name = plus one unnamed subprogram like the 
one above. Portable programs, however, should have no more than one 
unnamed BLOCK DATA subprogram. This subprogram may be used to inifialize 
any number of named Common blocks: the above S10cK DATA subprogram 
initializes Variables in COMMON blocks JOTA and KAPPA. 


At is important to specify the whole of a named common block even 
if you want to initialize just a few of ifs variables. Each named 
common block must be seen to have precisely the same length in every 
program unit in which it is declared: the lock DATA subprogam affords 
no exception. 


Phe following SL0cK DATA subprogram could be used to replace the sub- 
routine named CLEAR on page 79: 


BLOCK DATA 
COMMON / STACK f A(10) , IPOINT 


DATA IPOINT / Of 
END 





thus making it possible for the main program to do without the _initial- 
ization statement: CALL CLEAR. However, some programmers consider 
it a mark of bad style to initialize in this way any variable that 
might later be changed in value. 
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INITIALIZATION IN A BLOCK DATA SUBPROGRAM 
= INTRODUCING “FREE FORMAT” /NPUT = 


“According to Fortran 66 @not Fortran 77 which has a special type D 
characters may be stored in variables or array elements of any type. 












REAL R 
INTEGER z ae” 
DOUBLE PRECISTON D 

LOGICAL iL 

COMPLEX c 

DATA —_-R, I, D, L, C/ 5 * 3HABC 


but this practice leads to non-portable programs because different com- 
puters have words of different length. On a typical 16-bit computer 
the effect of the above DATA statement might be as depicted below. 
The indicates a blank or space character. 


R{AIBICWA TI D(AIBICVAVA774 “-|AIB] +~CIABIUWH7777 








“Phe safest fechnique for portable programs is to store just one character 
per integer variable , or per element of an integer array. Here is a 

BLOCK DATA subprogram to initialize the fortran 66 character set ina 
COMMON block named SIGMA: 


BLOCK DATA 
COMMON / SIGMA [ K(4?) 
DATA K(1), K(2), K(3), K(4), K(5), K(6), K(#), K(8), K(9), K(10) 
/ tHO, HL’, 142, 1H3,1H4, 1H5, 1H6, 1H7, 1HB,1H9 / 
DATA K(i1), K(12), K(13), KG), K(15), K(16), K(LW, K(18), K(19), K(20) 
/ IA, 1HB, {HC , IHD, 1HE , 1HF, 1HG, {HH , 1HI, 1HT/ 


DATA = K(21), K(22), K(23), K(24), K(25), K(26), K(23), K(28), K(29), K(30) 
/ 1HK, tL, 1HM, 1HN, HO, JHP, 1HQ, 14R , 14S, 1HT / 
DATA (31), K(32), K(33), K(34), K(35), K(36), K(3T)WA 
/ (HU, {HV, 1HW, LHX, JHY, HZ, tH */ 
DATA (38) , K(39), K(40), K(41), K(42), K(43), K(44),K(45), K(G6), K(4) 
J t=, 1H+, tH-, tHe, 1H/, IHC, 1H), 1H,, tH., AF / 
END 





Whenever a character is input it may be “looked up” in vector k(), and its 
subscript used for any necessary manipulation. There is a problem, however, 
in the “looking up”. Comparison of characters is nof trivial. Consider integer 
variables 71BW and JWW. Because the characters occupy the most signif- 
icant ends of I and J both variables hold enormous values. Furthermore 

the left-most bit is often used as a sign bit, so each integer value might be 
enormously positive or enormously negative. A statement involving IF(Z .£q. 7)... 
might cause the computer to subtract one value from the other in order to 
test for a zero result. If the sign bits were originally opposed this oper- 

ation might result in “integer overflow” and stop the program. The moral is 
to test the sign of each value and compare values only if the signs prove 
the same. Here is a logical function (used several times on subsequent pages ) 
for discovering if two integer variables or array elements hold the same 

character or not: 


LOGICAL FUNCTI@®N SAME(TI,; J) 
SAME «= .FALSE. 


IF((I .LT.O .AND. J.GE.0) .OR. (J .-LT.O .AND. I.GE.0)) RETURN 


CANE ee ERE xe true iF L&T ARE 
RETURN THE SAME, OTHERWISE 
END false 
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Using the subprograms opposite, here is a function to return a number 
@an “index” in the range 1 to 47 ) indicating what character occupies an 
integer variable. If this function returns zero it means the character 
Occupying variable NW is not in the Fortran 66 character set: 


INTEGER FUNCTION  INDEX(N) 

LOGICAL SAME 7 Saab 

COMMON / SIGMA / K(47) <I ~ Biock “31GmaA’ 
SHOWN OPPOS/TE 

D¢ 10 z= £, 47 

INDEX =I 


F (SAME (N, K(I))) RETURN 
CONTINUE Ames Feaetion 


INDEX = 0 = 
CHARACT: 


END 





‘fo illustrate a simple use of function INDEX() here is a program to 
read Q@ punched card €or line of data typed at a terminal 2 and print 
the values of unsigned integers punched into the card. The only charac- 
ters allowed are aigits and spaces. Each unsigned integer is assumed 

to be terminated by one or more spaces @ blanks). There is no con- 
straint on the particular fields in which the separate items are punched. 
This is called free format input. 


WWe anticipate page 112 which explains the A-Fformar in detail. In the 
example below the FORMAT (80A1) simply causes the character from each 
card column to be read and stored in the most significant end of the 
corresponding array element. In other words vector L(80) is made to 
contain an image of the punched card = ome character per element = 
as sketched below. The sketch assumes a computer that can store up 
to two characters per integer location. 


L£(80) , SPACE 
SPACE / 37/ 


IF (I .G6T. 80) STOP 

kK = INDEX (L(Z)) 

IF (Kk .£Q. SPACE) G® T 10 

IF ((K.LT.1) .OR. (K.GT. 10)) STOP 1 





= 0 


IF (I .GT. 80) GO Tb 30 

K = INDEX (L(Z)) 

IF (Kk .EQ. SPACE) GO TO 30 

IF ((k .LT. 1) . OR. (K .GT. 10)) STOP 2 
INTGER = 10% INTGER *# K - 1 

GO TO 20 

WRITE (6,200) INTGER 

FORMAT (1X, 10OHINTEGER IS, T10) 


Some input: 123 0046 9999 the output: [> INTEGER 1S 123 . 
© INTEGER {ts 46 o 
© INTEGER IS 9999 o 
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wiv INPUT OF ROMAN NUMERALS 70 
ILLUSTRATE USE OF THE DATA STATEMENT 


Assume all valid Roman numerals are composed of the following elements = 
never more than one from each box. @In fact classical Rome preferred to see 
Il rather than IV = the subtractive principle was a later development. ) 


D = 500 V:5 
{00 DC = 600 I=t VI=6 
CC=200 DCC=700 Il=2 VIL= 7 


CCC = 300 DCCC = 800 Ill=3 VIII=- 8 


Ht Ace CD = 400 CM= 900 V=4 IxX=9 





The logic of this program is contained in the following state table 
(or symbol-state table ): 


i a a os ee ee ee 
5 





(000% 02] 500% 03| 100% 09] 50%05 | 10810 | 5& o7 | 1% if | 
| 1000 & 02 | 500% 03] 100 & 09| 50%05 | 10% 10 | 5% 07 | 
é 





ro 

02 | | 1 & if | 
fos] error | érror [100% 09| s0&05 | 10K10 | 5a o7| 18 11 | 
Fos | error | error [1008 04| 50&05 | 10810 | sk o7] 18 u | 
fos | error —| error | error | 50406 | 10810 | 5807] 18 t1 | 
‘06 | error__| error | error | error | 10806 | 5&07| 1811 
mire 
06 | error _|error_| error | error 
roe | 800 & 05 | 300K 05| 100k on | 50R06 | 10810 | 54 08] 18 tt] 
r10_[ error Soe] 18 ut | 
ru] error | error_|_ error 





“STATE 








Sake the Roman number CIX as an example: begin with a value of zero. 
You are initially in state £4 @where the arrow is ) so look down from 
symbol C and find 100409 which says: “Add 100 to the value and 
change state to 09”. So add {100 to zero and move the arrow to 09. 
Now look down from symbol 1 and find 1& 41. So add f to the 
value @€ /004+f= 101) and move the arrow to state 1. Finally look 
down from symbol X and find 8 & 00. So add & fo the value 
(101 +8 = 109). The 00 means you've finished. Thus CIX is 109. 


Phe Program is designed to read one card = or row of data froma 
screen = containing Roman numerals separated by spaces: 


ROMAN NUMBER = 


the output: |° Roman numBer = 
ROMAN NUMBER = 


Some input: 


ROMAN NUMBER = 





[nlere is the program: 













REAL S(11,7) 
LOGICAL SAME 

INTEGER  C(#), L(80), STATE, SYMBOL, SPACE, PRINTR, KARD 
DATA SPACE, KARD, PRINTR / 1H , 5,6/ 

DATA C(t), C(2), C(3), C(4), C(5), CCE), CCF) 

/ HM , 1HD, 1HC , MHL , GHX , SHV, tHr / 


Elements of the state fable are packed in Ré&AL array S(11,7). Each element 
is packed as 100.0 times the first part plus the second part. A REAL orray 
is used because on many computers the size of an integer array element is 
limited to 32767. Error states are represented by -1.0 . 
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ere is the symbol-state table: 


Whe executable 








DATA S(1,1), S(1,2), $(1,3), S(1,4), 5(1,5), $(1,6), S(1, #) 
/ 100002., 50003., 10009., 5005., 1010., 507., Wt. / 
DATA S(2,1), S(2,2), 5(2,3), 5(2,4) , (2,5) » 5(2,6), 5(2,#) 
/ 100002., 50003., 10009., 5005., 1010., SO7., 111. / 
DATA S(3,1) , 5(3,2), $(3,3) » $(3,4) , 5(3,5) , 5(3,6) , $(3,%) 
-1. , il. , 10009., 5005., 1010., §507,, 111. / 
S(4,1) , 5(4,2), 5(4,3), 5(4,4), 5(4+5), 5(4,6), 5(4,9) 
-1. , 4. , 10004, 5005., 1010., 507.5 ft. / 
5(5,1), S(5,2), 5S(5,3), S(5,4), $(5,5), 5(5,6), 5(5,7) 
-{1. , -4. , -1. , §006., 1010., 507, {f1- / 
56,1), S(6,2), S(6,3), 5(6,4), $(6,5), $(6,6), S(6,7) 
“{. 5 he 5 he, dy , 1006, 507., 1. / 
S(#,1), S(#,2), S(#,3), S(#,4), S(7,5), S(7,6), S(#,#) 
he og: HW gy, Hho Hb ig “Ee g, SOB, difat 
5(8,1), 5(8,2), 5(8,3), 5(8,4), 5(8,5), $(8,6), 5(8,7) 
“hg Sh g Shag Steg —y Hh iy 108.7 
5(9,1), $(9,2), $(9,3), 5(9,4), S(9,5), 5(9,6), 5(9, 7) 
80005., 30005.,10004., 5006., 1010., 508., 111. / 
$(10,1), S5(10,2), $(10,3), 5(10,4), 5(10,5), $(10,6), S(10, 7) 
-1. -1. , 8007., 3007., 1006., 5086., 111. / 
S(it, 1), 5(11,2), $(11,3), $(11,4), $(1f, 5), 5(11, 6), S(1,#) 
“1. , ft. , tl. , ft. » 800., 300., 108. / 
part, below, uses function SAME(,) from the previous example: 
READ (KARD, 100) L READ ALL 80 CHARACTERS 
FORMAT (80A1) Ql as IN PREVIOUS EXAMPLE 
NEXT = 0 
WEXT = WEAT + 1 
IF (NEXT .GT. 80) STOP 
ZF (SAME ( L(NEXT), SPACE )) GO 76 10 


STATE = 1 

NUMBER = 0 <Hymiiauze ) 

NEXT = NEXT — 1 

NEXT = NEXT +14 

IF ((NEXT .GT. 80) .OR. (SAME(L(NEXT), SPACE))) G@ TO 50 
Dg 20 SYMBOL = 1,7 

IF (SAME ( L(NEXT), C(SYMBOL) )) GO TO 30 
CONTINUE 


idea 2 <4 wor M,D,C,L,X,V or I 


lenTRY = S(STATE, SYMBQL) ED 
IF (ENTRY .EQ. -1.0) STOP 3 FROM STATE 7. 
F CCNEXT .GT. 3) .AND. 


SAME ( SYMBOL, L(NEXT-1)) .AND. LF CHECK NO MORE THAW 
SAME (SYMBOL, L(NEXT-2)) . AND. 3 SUCCESSIVE SYMBOLS 


OF SAME KIND 
SAME (SYMBOL, L(NEXT-3)) ) STOP 4% 


PART = AINT (ENTRY / 100.0) 


NUMBER = NUMBER + IFIX(PART) 


STATE = IFIX (ENTRY — 100.0 PART) <*4SECOND PART OF ENTRY 
IF (STATE .NE. 0) GO TO 40 © NEW STATE 


, 110) 
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EXERCISES) CHAPTER 9 


Do 1 


9.8 


9.3 
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Shange the integer function DAYS(MONTH,YEAR) on page 86 so that 
it shares a named common block with a &8L0CK DATA subprogram. 
Write this subprogram so that it initializes the variables F&B and 
LEAP, and the array elements M(1)~to M(12). 


Recast the program on page 89 as a subrovtine for reading 
integers written in free format. Include a logical argument OK 
to be returned trve if the call succeeds; false otherwise. for 
example OK should be set false if the subroutine finds a letter 
among the digits ¢ possibly letter 0 typed instead of a zeroy. 
If the subroutine reaches column 80 this should indicate the end 
of an integer: a new card should be read automatically if the 
subroutine is called again. 


CALL IREAD(N, OK) 


WYNEAT INTEGER 
DELIVERED IN 
VARIABLE ‘N? 






Allow for optional plus or minus signs if you want this sub- 
routine to be really useful. 


Write Q@ converse program to the one on the previous page. The 
program shovid read unsigned decimal integers and print 
corresponding Roman numerals. Such a program doesn’t 
require a symbol - state table and shovid be easier to devise 
than one to convert Roman numerals to Arabic. DATA state- 
ments may be useful; for example: 


DATA K(1), K(2), K(3), Kh), K(5), K(6), K(?¥) 
yf 1HM, 1HD, JHC, FHL, THX, 1HV, IHF / 


but remember your program would not be fully portable if you 
Stored more than two characters in one integer variable: 


bara 7 aver 


Because we have not yet covered formatted output, print the 
results in a column down the left of the output page: 


eooo0oo0ooco090 


© 
© 
o 
© 
te] 
oO 
oO 
4] 





WRITE (6,200) k(t) 
200| |FORMAT (1X, Al) 


EDO 





ONTO 
OUTPUT, 


READ 

WRITE 

GENERAL 

I/0 LIST 

FORMAT 

FORMAT (CONTINUED ) 
RUN-TIME FORMAT 
GRAPH @ AN EXAMPLE ) 
DESCRIPTORS 

NUMBERS /N DATA 
FRUSTRATED OUTPUT 
DESCRIPTOR Fw.d 
DESCRIPTOR Ew.d 
DESCRIPTOR Dw.d 
DESCRIPTOR Gw.d 
SCALE FACTOR nP 
DESCRIPTOR Iw 
DESCRIPTOR Ly 
DESCRIPTOR Aw 
HOLLERITH LITERAL WH Ah...h 
BLANKS @SPACES) wK 
FREE FORMAT (AN EXAMPLE ) 
EXERCISES 


INPUT IN FORMATTED € READABLE ) FORM 
INPUT IN UNFORMATTED @ BINARY } FORM 
Whe introductory example showed: 
READ (5, 100) DIAM, HEIGHT, COVRG 
100| FORMAT _(3F 10.0) <3 Fieupis 0 counns PAM 


which caused numbers on a waiting card @or card image) to be read into 
variables named odD/Am, 








ONE CARD 
1S ONE 
"PHYSICAL 
RECORD’ 


‘Whe general form of the READ statement is: 


READ (unit , format) list 
i READ (unit ) list 


where: 

unit is either an integer constant € digits only ) or the name 
of an integer variable @ not an array element 2. unit 
denotes a connection to some peripheral device on which 
data are waiting to be read. The connection is made by 
command to the compuvter’s operating system =» perhaps bya 
PROGRAM command placed immediately before the main 
program. Any unit may be used for formatted or un- 
formatted input; not a mixture of both. Conventionally 
unit 5 denotes a card reader. 


format may be omitted altogether # in which case the RAD state- 
ment causes a complete wnformatted record to be read. 
This record shovid be waiting on a peripheral device such 
as a magnetic-tape deck. The waiting unformatted record 
= possibly spread over several physical records such as blocks 
of magnetic tape = should match the (ist identically. 


format, if present, may be the label @ digits only ) of a 
FORMAT statement @ page 98) describing the layout of a 
formatted record or sequence of such records. The pattern 
of items in the waiting data should match not only the 
list identically but also the description given by the 
FORMAT statement. 


format may also be the name of an integer array € con- 
yeniently an integer vector } containing the characters of 
the format description @ page 102). 


list (called the input-output list or 1/0 list 2 is a list of 
variables into which the waiting data are to be read. 
The example illustrates a simple 1/0 list: D/AM,HEIGHT, COVRG 
But there may also be array elements in the list, and 
their subscripts may be generated automatically. Page 96 
defines the I/o list in detail. 


[3very READ statement causes a fresh record to be started; for example 
the next card to be read. Any items not demanded by the 1/o list 
are Php to the program. You can’t retrieve subsequent items on the 
same “card for example. A subsequent RkéAD statement would cause the 
next card to go under the reading head. 
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OUTPUT IN FORMATTED @€ READABLE ) FORM 
OUTPUT IN UNFORMATTED €B/NARY ) FORM 
‘he introductory example showed: 


WRITE (6,200) NPOTSL” YS 
FORMAT (1X, QH Y@U NEED, 12, SH POTS) 


which caused the value stored in wNPOTS to be written as part of a single 
output record: ONE LINE 


NPors[ 3 | 1$ ONE 


"PHYSICAL 
RECORD’. 






YOU NEED 3 POTS 









VEIELD OF TWO CHARACTER 
POSITIONS; I2 








he general form of the wrir— statement is: 


ne WRITE (unit, format ) list 
WRITE ( unit) list <E{_UN FORMATTED) 
where : 


unit is as described opposite except that the peripheral device 
must be capable of output; for example a printer or tape 
deck but not a card reader. Conventionally unit 6 denotes 
the line printer. 


format may be omitted altogether = in which case the wasre state- 
ment causes a@ complete unformatted record to be written on 
the nominated unit which may be connected to a magnetic- 
tape deck, disk file or paper-ftape punch. The record to be 
written is as long as the dist demands = possibly spread 
over several physical records @ disk sectors or blocks of 
magnetic tape ). 


format, if present, may be the label (digits only » of a formar 
statement (page 98) describing, in conjunction with the J/o list, 
one or more physical records eg. printed lines; punched cards ). 


format may also be the name of an integer vector which 
Stores the format in character form @ page 102). 


list @1/0 list ) has the same form for output as for input. ist is 
fully described on page 96. 


CONCERNING READ AND WRITE STATEMENTS: 
As NsRAL 1. “READABLE” VERSUS “BINARY” 2. EMPTY I/O LISTS 
A formatted record = on printer or screen» is in character form readable by eye. 
An unformatted record is more compact than the corresponding formatted record but 
would appear as gibberish if it were possible to send one to the printer. 
Unformatted records are used, typically, to transfer intermediate results between 


Qrrays and disk « withovt conversion to and from character form, hence without 
losing speed and accuracy during transfers. The information stays in binary. 


An some circumstances the 1/0 list may be omitted from RéAD or waiTe Statements: 


READ (5, 100) <i SKIPS A COMPLETE PHYSICAL RECORD 
(6:9: THE NEXT CARD) IGNORING THE 
FORMAT (3F 10.0) NOMINATED FORMAT 


— al KAS {OVER fig Bmery RECORD 
WRITE (7) NOT ALLOWED witnour AN X/O “LIST 


WRITES MOLLERITH LITERALS IN THE FORM 
WRITE (6,200) 


Bur 
TERMINATES On, MEETING A DESCRIPTOR (seta As I4) 
200 | |FORMAT (18H THIS GETS PRINTED, 14 



















+ J3H THIS DOESN'T ) 
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DENOTING A STREAM OF /TEMS TO BE READ OR WRITTEN 
lV °) SST (CONSTANTS & EXPRESSIONS ARE “Ow ALLOWED ) 
She introductory example illustrated two simple J/o lists: 
ee eS ca) ah COVRG 
Whe general form of (ist is as follows: 


simple -list 

( simple-list ) 

(list, control = initial, terminal ) 

(list, control = initial, terminal, increment ) 


or a list of such lists separated by commas, where: 


simple-list consists of names of variables, names of arrays, names of 
array elements, or any combination of these. The names 
are separated by commas. Constants and expressions are 
not allowed in the 1/0 list. 


control \ these are as defined for the po loop € page 40) and have 

initial ( the same implications as in a DO loop. The range of 
terminal \ this implied DO loop is the preceding list contained 
increment J within the same pair of brackets. 


(ere are more examples of simple-lists: 


INTEGER — KARDS, LINPRN, TAPE1, TAPE2 
DATA _-KARDS, LINPRN, TAPE1, TAPE2/ 5, 6, 4,3/ 






READ (KARDS, 300) X,Y, A(2,1), A(2,2), A(1,I) 


WATE (TAPE2) MATRIX WSJ COMPLETE ARRAY NAMED 
‘MATRIX? TRANSMITTED 
BY COLUMNS 


Whe definition permits dist to be a Simple-list in brackets, but some Fortrans 
@including Fortran 77 ) object to unnecessary brackets: 


READ (KARDS, 300) (X,¥, A(2,1), A(2,2), A(I,7) ) 
WRITE (TAPE2) iG MATRIX 
AF Unnecessary 7 BUG 
BRACKETS Be 


[hlere is a_ statement illustrating the implied Do loop. This must be in brackets. 
| | |wezre (LINPRN, 400) (VECTOR(Z), I= 4, 10,2 ) 


oe VECTOR(4), VECTOR(6), VECTOR(8), VECTOR (10) 


And here is a@ READ statement illustrating a two-deep nest of implied Do 
loops and a weliTe statement illustrating a three-deep nest. The Fortran 66 
Standard does not forbid it, but Some Fortrans do not allow a_ nest 
deeper than three. 


READ (TAPE1) (C A(I,J), I =1,60), T= 1, 100) sb ne 
WRITE (TAPE2) ((( B(I,J,K), T= 1,20), J=1,30), K=1,50)| <ea RECORDS 


[Notice the recursive forms of the third and fourth definitions above; 
list is defined in terms of dist, thus permitting nested loops. These 
definitions account for the pattern of brackets and commas in the 
examples immediately above. 
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“All four of the forms of the I/0 list defined opposite have been illustrated 
individually, but notice that dist may be a list of lists: 


WRITE (LINPRN, 500) (A(1), 81), C(I), T=1, 10), X,Y,z 
WRITE (LINPRN, 600) ( VEC(Z), T= 1,6), ( COL(J), T= 1,3), P,Q 


“Phe items controlling the implied po loop do not have to be _ integer 
constants; they may be names of integer variables. Page 50 describes 
the allowable forms of subscript involving integer variables. 


| |wRITE (TAPEZ) (v(2#I +3), I= J,K,L) 

f. Wa e/- 
(2#I +3) VALUES CURRENT 
MAXIMUM COMPLEXITY WHEN THIS WRITE 


STATEMENT IS 
OBEYED 


















Or should be allowable to read subscripts and use them to address 
elements of an array on the assumption that fortran works strictly 
from left to right along the 1/0 list: 


READ (TAPE1) I,J, A(Z, J), B(J, IZ) 
READ (TAPE2) K,L, (V(I), T=K,1L) 


but do not try to change parameters of an implied DO loop whilst it is 
still looping: 
| Mee | [reap (rAPE1) CK, L, V1), 1=K,L) OR 


Vou are not allowed to transmit an array with adjustable dimensions. 
using only its name (in the manner of array MA7R/X illustrated opposite): 


SUBROUTINE CQPYC(A, I,J) 
DIMENSION A(I,J) 
READ (TAPE1) A 


WRITE (TAPE2) A 7th, 
RETURN 
END 





but this subroutine could be corrected by changing the RE&AD and weriTE 
instructions as follows: 


READ (TAPE1) (CA(K,L), K=1,1f), L=1,J7) 
WRITE (TAPE2) ((A(K,L), K=1,I), L=4,7) 


Gaany computers impose a limit on the length of an unformatted 
record, so find ovt what limit applies before attempting to transmit 
long unformatted records like those illustrated opposite. With formatted 
records it is a mistake to transmit more characters than the chosen 
medium can hold @e.g. 80 columns on a card; 132 characters on a 
printer ). But the stream of items generated by an I/0 list may be 
organized into a sequence of manageable records by FoRmAT statements 
as later described. 





FOR y~\ /\ DESCRIBING LAYOUT /IN “READABLE ” FORM: 
G EXHAUSTS THE 1/0 LIST =» BY “RESCAN” IF NECESSARY ) 
Phe introductory example showed: 


READ (5,100) DIAM, HEIGHT, COVRG 
goalie (3F 10.0) INPUT RECORD 
i (6, 200) NPOTS 


ORMAT (1X, 9H YOU NEED, [2, 5H POTS ) 


where fhe lines labelled 100 and 200 illustrate formaT statements. These 
may be interspersed among the executable statements but must always 
be labelled. 


(8 xecution of a formatted READ or wRITE statement causes a new record to 
be started: the subsequent input or output in that record is “driven” by 
the nominated format. On input all items are in “readable” @ character) 
form and have fo be converted to binary form. On output all items 
ore in binary form and have to be converted to character form. These 
conversions € from character to binary form and vice versa) are controlled by 
descriptors such as the 3F10.0 and I2 in the examples above. Descriptors 
are dealt with exhaustively later in this chapter. 


She form of the FORMAT statement is: 


FORMAT ( description ) 


where GQescription may take any of the following forms: 

@ a list of descriptors separated by commas 

@ as above but with one or more commas replaced by a slash 
or group of slashes each slash starting a new record. @ The 
slash or group of slashes may also be inserted at the end or 
the start of a list where there is no comma to replace, but 
Such practice is not recommended. ) 

@ as above but with one or more descriptors replaced by 
descriptions in brackets. The left bracket may be preceded 
by an integer saying how many times the content is implied 
@ the repeat count ). Omission of this integer implies once. 
@This definition is recursive; Gescription is defined in terms of 
description. » 


and where descriptor is any of the following Gall of which are described 
in detail at the end of this chapter D: 


nPtF w.d = fixed-point format 
nPtEw.d =  £-format 
nPtDw.d *s double-precision format 
nPtGwed = generalized format 
tlw = integer format 
tLw = logical format 
tAw ~ character format 
wHAhh...h =  Hollerith literal 
wx = spaces ¢ blanks ) 


and where: 
is an optional scale factor, 10”, described on page 110. 
Omission implies m=0; a scale factor of unity. 


(¢) is an optional integer saying how many times the 
descriptor is implied: omission implies once. 
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is an integer specifying the number of character positions 
s the field width= of the item fo be read or written. This 
integer must be greater than zero. Also w must be greater 
than a, where: 


is an integer specifying the number of digits assumed after 
the decimal point if no point is written in the data = 
or the number of decimal places to print on output. 

This integer must be included even if it is zero. 


each A € of which there are w in the field } is a Hollerith 
character * preferably from the fortran 66 character set, 


© 8 © 


Qn the simplest formats there is a one-to-one correspondence between items in 
the 1/0 list and the various descriptors: 


READ (5,100) DIAM, HEIGHT, COVRG 
100) |FORMAT (F10.0, F10.0, 10.0) 


but wherever a sequence of identical descriptors occurs the repeat count, t, may 
be used as illustrated by 3F10.0 opposite. Similarly the following two formats, 
300 and 400, are equivalent: 


300| |FORMAT (214, 3F 10.0, I¢) 
400) |FORMAT (14,14G,F10-0, F10.0, F10.0,I4 ) 


and making use of the third form of description defined opposite, the following 
two formats are equivalent. The format labelled 500 is “nested”: 


a ee (2F 10.2, 2(2F10.0, I%)) 

600} |FORMAT (F10.2, F10.2, F10.0,F10.0,14,F10.0, F10.0,I4) 
In fact these two formats are not quite equivalent. If the 1/0 list supplied more 
than eight items there would be a difference in behaviour as explained below. 


Descriptions may not be nested more than three deep: 


| 700] [FORMAT (\1X, 3(F10.2, Ih, 2GF 10.0))) 

Ce C—~ ——~ 

Woctmnt) Were 2) aD), 
Phe 1/0 list should provide enough items to satisfy descriptors in the nominated 
format. Each record is closed as soon as the last element in the I/o list 
has found its descriptor. In this example the word £x7TRA wouid not be 
added to the output record: 


ee (6, 800) V1,v2 
OW 


FORMAT (1X, 3F10.2, 6H EXTRA ) 
But if the 1/0 list provides more items than the format can deal with then 
the format is automatically rescanned as depicted below: 


Leal WRITE (6, 900) (VECTOR (I), L= 1,30) 
FORMAT (1X, FI10.2 Ep 
—— NEW RECORD 


where each rescan initiates a mew record. Thus on a line printer all thirty 
values above would appear as a column down the left-hand side of the 
ovtput page. In a nested format the rescan @and start of new record ) 
begins at fhe left bracket matching the dast but one right bracket = taking 
account of the repeat count if there is one: 


Pe] 






NOT ENOUGH VARIABLES 70 
ZER_ SATISFY 3F 10-2 




















LAST BUT OWE 
WRITE (6, 1000) (VECTOR(I), I = 1,30) pZRRIGHT BRACKET, 
FORMAT (1X, F10.0, 2(F 10.1, F10.2), F10-4 


RESCAN FOR NEW RECORD earmreee 









Ccontinued D 


99 


FORAY (aaspnvgo ) conrroiline Sue? LINE PRINTER 


[3very time format control meets a slash in place of a comma a new 
record is begun. A group of slashes thus produces blank records on 
output s or skips over complete records on input. Thus: 


DATA IVEC(1), IVEC(2), IVEC(3), IVEC(G), IVEC(5), IVEC(6) 

/ L , 2 , 3 , 4 , 5 ’ 6 7 
WRITE (6,100) (IVEC(J), T= 1,6) 
ORMAT (14,14 / 14,14¢,14// 14 /” ) 

















Should produce: 


Se SECOND RECORD 
<i) BLANK THIRD RECORD 


6 <I FOURTH RECORD 
<EBLANK FIFTH RECORD ?? 






A group of k leading or trailing slashes shovid produce & blank records. 
A group of k slashes replacing a comma should produce k-1 blank 
records. Some Fortrans, however, treat leading and trailing slashes 
differently so it is best to avoid leading and trailing slashes for the 
sake of portability. 


Lr an item of data being input fails to match its associated descriptor a 
warning message is usually printed. In some Fortrans this does not 
cause the run to stop but does leave the input variable undefined 

¢ containing rubbish ». 


REAL A 


READ (5, 200) A ~<&& MISMATCH LEAVES ‘A’ 
FORMAT (14) tb, FulL OF RUBBISH 





On output a mismatch usually causes the associated output field to be 
filled with asterisks or question marks: 


Same 


» 9H MISMATCH) 





300| |FORMAT (1X, I4@ 


O *%8% MISMATCH © 





[nlowever there is not necessarily a mismatch when a variable of one type 
is read or written using a format descriptor of another type. Both 

real and compjex numbers may be transmitted using -, £ and G de- 

scriptors. And because Fortran 66 has no variable of type CWARACTER 

¢ Fortran 77 does have this ) Hollerith data have to be stored in the guise 
of integers: 


| 
















INTEGER IVEC(36) 
READ (5, 400) IVEC 
FORMAT (36A2) 



















rec (1) [MIA JWls SKETCH 
IVEC (2) HG! ASSUMES A COMPUTER 
tvec (3) [YE TWAT STORES 


TWo CHARACTERS 
PER INTEGER 
LOCATION 





MASQUERADE 


1 


IveéC Ca AA < BLANKS 


100 


A complex number requires two descriptors in the format: 


COMPLEX C 
READ (5,500) C 
WRITE (6,500) C pAITy 
500] \FORMAT (1X, F9.0, F10.0 


SRL es S 4278 3.2 i 
Co 1} 24 10) ° () 


De an output channel is connected to a line printer @ conventionally unit 
6 is so connected ) the first Character is not printed. The first character 
position is left blank; the character itself is used to contro| the line- 
feed mechanism of the printer as follows: 


FIRST CHARACTER OF RECORD CARRIAGE CONTROL 


<E{TYPICAL OF A FIRST Record) | Skip to the first line of a new 
page before printing the record 


0 <KEze%0_) skip fwo lines before printing 
blank print this record on the next line 


+ SIE eT Ber ar) print record on the same line 
Se (i.e. overprint previous record ) 
other characters non-standard; but in some Fortrans 
Sy. treated as blanks 


[Previous examples of output formats begin with 1X to ensure each record 
begins with a blank. Another convention is to replace the /X with 1H . 


[zoo |rommar ix 7 = [200] roenar Cin. 7 
CLANK) 


Lr a Formar statement consists only of X and H descriptors make sure 
the associated wrir—E statement Aas no I/o list: 


DATA MTAPE /4/ ae 

WRITE (MTAPE, 100) I 

FORMAT (6H TRAMP, 2X_) 
—— 


















































Fortran would keep rescanning the format « hence keep writing “tramp, tramp, 
tramp” on the magnetic tape = vainly trying to find a descriptor for Z. 
This could go on until the tape ran ovt. More destructive is: 


DATA  LPRINT /6/ 
200 


WRITE (LPRINT, 200) I ie 

FORMAT (432 (1H+)_. ) 
which could cut through the paper, mash the ribbon, even damage the 
printer. @ I have not tested these routines and trust my readers won't. ) 
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ASSEMBLY VIA DATA 
RON = WA FOR VAL On READ STATEMENTS 


&\s already explained; format in READ(unit, format) or WRITE unit, format) 
may be the name of an integer array containing the format in 
character form. The keyword FORMAT is not itself included 
among the characters in this array. 


(RJollerith literals are not permitted in such formats because the number 
of trailing blanks in each array element depends on the make of 
computer being used. However it may be possible € but ot standard ) 
to get Qway with Hollerith Jiterals one or two characters long. 


SPhe introductory example illustrated two formats: 
100\ |FORMAT (3F10.0) 
200| |FORMAT (1X, 9H YOU NEED, I[2, 5H con 
aol eas 
OGD 
Apart from the two Hollerith literals these formats could have been 
assembled and used as follows: 





INTEGER IA(4), IB(4) 
DATA IA(1), IA(2), IA(3), IA(4) 
/ 2H(3, 2HF1, 2HO., 2HO) / 


DATA IB(1), IB(2), IB(3),I8(4) 


AD 2H 
[[eead (5, IA) DIAM, HEIGHT, COVRG j nN A 
ALL BLANKS /GNORED: 
TWE SKETCHES ASSUME 
|__| [WRITE (6, IB) NPOTS A COMPUTER THAT STORES 
FOUR CHARACTERS PER 
INTEGER LOCATION 





®r these formats could have been read as additional data: 


READ (5, 300) IA 
READ (5, 300) IB 
300) |FORMAT (_A2_) 


RESCAN FOR EACH 
PAIR OF CHARACTERS 
THE 'A2* MEANS TWO ) 












HOLLERITH CHARACTERS 





Phe ingenious programmer can make the layout of results dependent 
on the magnitude of values to be printed. If is only necessary to 
manipulate the integers into which the Fortran 66 character set has 
been stored » for example as on page 88. Remember you are not 
allowed fo assign Hollerith constants: 


but you can achieve the same effect by assigning integer values: 


DATA IX / 2H1. / 
IA(3) = IX 
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GAP 


AN EXAMPLE 70 ILLUSTRATE A USE OF 
RUN-TIME FORMATS FOR PLOTTING 


Phe following program is designed to plot a crude graph of y=cos(x) 
by printing asterisks as points on the curve. The scales have been chosen 
so that the graph would fit a modest 72-column printer @ or VDU). 





INTEGER ROW(61), STAR, BLANK, IPLACE, ANGLE, PRINTR 
REAL RADIAN, COSINE 


Le 
DATA STAR, BLANK, PRINTR / 1H#, 1H’, 6 / 
FE 


USED 
WRITE (PRINTR, 100) 


FORMAT (1H1, 41HGRAPH OF COSINE (X) FOR X=O0 TO 360 DEGREES) 


DO 10 T= 1, 61 
ROW(I) = BLANK 


DO 20 J*1, 25 


ANGLE = 15 * (J-L) Gy 
FLOAT (ANGLE) * 3. 141593/180.0 
COS (RADIAN) 


* 31.0 <SEAISCALE ROvGHLY TO 
= STAR CENTRE OF PAGE 
WRITE ( PRINTR, 200) ANGLE, ROW 


FORMAT (1H , I3 , 61A1 ) 
= BLANK 


She output should look like this: 
GRAPH OF COSINE(X) FOR X=O TO 360 DEGREES 


0 
15 
30 
45 
60 
15 
30 

105 
120 
135 
150 
165 
180 
195 
2410 
225 
240 
255 
270 
285 
300 
315 
330 
345 
360 
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[aye R A |? 1,0 RS RECAP/TULATION = AND SUMMARY 
Whe introductory example showed: 


READ (5,100) DIAM, HEIGHT, COVRG 
100) |FORMAT (3F 10.0) 


THREE CONSECUTIVE FIELDS 
OF 10 COLUMNS EACH 

















WRITE (6,200) NPOTS 
FORMAT (1X, 9H YOU NEED, I2, SH POTS ) 


eS (ay \ oe 

HOLLERITH INTEGER IN \ HOLLERITH 
LITERAL iN NEXT NEXT FIELD & LITERAL IN NEXT 
FIELD OF 9 COLUMNS) OF 2COLUMNS) FIELD OF 5 COLUMNS 








BLANK FOR 
CARRIAGE 
CONTROL 










YOU NEED 3 POTS 





The programs @and fragments of programs ) in earlier chapters used only: 


@ 1X to make the printer start a new line 

@ 10.0 for reading real numbers; F410.2 for writing 
to two decimal places; F10.4 to four places etc. 

@ 12 for writing integers in a field of two columns; 
Z10 in a field of ten columns etc. 

@ Hollerith literals for annotating results 

@ Az for reading a pair of characters into a 
variable; 8041 for reading a whole line of 
Characters into successive elements of an 
integer array. 


Qn the following pages all available descriptors @ plus the scale factor 
”P ) are explained in detail. The notation already introduced on 
page 98s is again used to summarize all descriptors as follows : 


nPtFw.d 


«= fixed-point format 
NPtEw.d = £- format 
”2PtDw.d = double-precision format 
NPt[GW.d = generalized format 
tlw # integer format 
tLw = logical format 
tAw *« character format 
WHAh...4 = Hollerith literal 
wx = spaces ¢ blanks p 
where: P is an optional scale factor; Z says how many times 
the descriptor is implied @ omission implies once ) ; Ww gives 
the field width; @ gives the number of decimal places 
implied if no point is written; h represents a_ Hollerith 


character =~ preferably from the Fortran 66 character sef. 
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FORM OF NUMBERS CONVERTED 
BY F,E,D,G DESCRIPTORS 


Data processed by f,£,D or G descriptors may be written in various 
ways. The essential thing is to confine each number to the 
field specified by w in its corresponding descriptor. 


Phe simplest form allowed is the ordinary decimal number «signed 
or unsigned: 


2.2,+,5 64,+,0 #22.5 -0.7 
A trailing or leading decimal point is permissible but not nice: 
64. -.7 


Numbers may also be written in exponent form: 
1, 5&6 or - 1. 5E%+6 for -1500000; 1.0£-2 for 0.001 


— ee ee ee ee 





The & may be omitted if the exponent is signed: 
- 1.5.4.6 1-0-2 





Furthermore it is permissible to write D in place of £& whether the 
number is destined to be a double-precision variable or not. 


PANS in the field count as zeros. Thus: 


: i Ha a implies 0.1020,3040,..00 
where fhe leading zeros are ignored and the trailing zeros do no 
harm. A _ fotally blank field is treated as zero by F,£&,D,G@ and 
I descriptors alike. 





L\t is safest ALIAS fo include a decimal point in numbers 
(especially those in exponent form ) to be processed by F,£,D,G 
descriptors. The decimal point in the data overrides an implied 
position of a decimal point located by a@ @ fw.d, Ew.d, etc. ). 
Not all Fortrans agree what £7.3 implies if the corresponding item 
of data is, say: Ra 


OUR OWN 
520E+6 or 52046 me ewe” 


but with 0.520€6 there is no doubt of the interpretation: 0.520x 10°. 


FRUGURATED COUR BMPedtPe 


“Phe value of w € in Fw.d, Ew.d,...,wx ) should specify a field 
wide enough to contain the number or message to be printed. 
If the field is too narrow, different Fortrans print different 
things = usually causing information to be lost = but few 
Fortrans cause the program to stop running. A _ typical outcome 
is a field full of asterisks or question marks: a frustrating 
result to receive: 


DATA VALUE / -12345.6 / 
100 













WRITE (6, 100) VALUE ae 
FORMAT (1X, 10H ANSWER IS, F6.1) 
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D3ZS9CRUP WOR F a ts 
wed 


“Phe item in the Wo list corresponding to this item should be of type 
REAL . The item may also be of type comptex provided that it 
corresponds to a pair of such descriptors. 


READ (5, 100) R,C 
100| |FORMAT (F10.0, 2F8.0) 





WRU T Cs The binary value in the computer is converted to charac- 
ter form and rounded to d decimal places. If the value is negative 
the number is prefixed with a minus sign. The number is output 
with as many preceding blanks as necessary to fit the number into 
a field of w character positions « right justified: 











A(4q) 
At) , M2) . A(3), A(4) 
/ -2368.4 , #12.0, -17.90767, 37.5£-2 / 
WRITE (6, 200) (A(I), LT=1,4) 
FORMAT (1X, F10.2.) 



















NOTE: iN FORMAT F6-0 THE 
VALUE —2368+-4 WOULD BE 
PRINTED AS 2368. 

(DECIMAL POINT IN SIXTH COLUMN) 









RSADINS 3 The number in the next w character positions of the 
card @ or card image ) is converted to a REAL value and stored 
in the corresponding variable or array element nominated in the 
1/0 list. If the number being read has no decimal point then fhe 
last ad digits are assumed to represent the fractional part. In 
other words there is an implied decimal point just before the last d 
digits. If you punch a decimal point in the number, however, 
this overrides the implied decimal point. Always do this / 


REAL BC?) 
(, 300) (a(D.,. ££, 7) 





B(1) ]0.53621 | 
B(2)| -16-40 | 
B(3)[ -0-1 | 


B(4) 


B(S) 









IMPLIED DECIMAL POINT 
JF NOT OVERRIDDEN 





ProR E FLOATING-POINT FORMAT 
w.d 


‘fhe item in the 1/0 list corresponding to this descriptor should be of 
type REAL. The item may also be of type compzex provided that it 
corresponds to a pair of such descriptors: 


REAL OR 
COMPLEX C 

READ (5, 100) R,C 
100\\FORMAT (E10.0, 2E8.0) 














WR PYNS 3 The binary value in the computer is converted to character 
form as _ follows: 


digits Et ee 

- digits E- ee 
where digits represents the ad digits specified; ee represents a two- 
digit decimal exponent. Valves are rounded to @ decimal places and 
the resulting number is right justified in a field of w character 
positions. 


Working from right to left the resulting number includes the two digits 
of the exponent, a sign for the exponent, fhe letter £, the d digits 
of the value, a decimal point, a minus sign or a blank in front of 
the value, enough blanks to complete the field width w. Some For- 
trans replace the £ with a D3; others with a blank. Yet others 
replace the E with an extra digit to the exponent. Some Fortrans 
precede the decimal point with a zero; others with the first significant 
digit €e.g. 4£.23&+00 instead of .123E+01 ). Some Fortrans precede 
positive values with plus signs instead of blanks. Such variations are 
permitted by the Fortran 66 standard. But if you make w greater 
than d@+6_ the field should be wide enough to contain the result 
whatever the style of output. 


REAL  A(4) 
DATA A(1), A(2), A(3), A(4) 
/ 46.573, -58796-36, 3F7-5E-2, 0.0066 / 


WRITE (6, 200) (ACI), T= 1,4) 
200} |FORMAT (1X, £12.6,) 
— 
















(a) *7T6E5730E+02 
O -.587964E+05 
ie} +375 000E+t00 
oO -680000E-02 


sooo 


CARRIAGE 
CONTROL 


RIZADANG, : The £ descriptor behaves as the Ff descriptor, but 
remember always to include a decimal point when writing an item 
of data. @ One famous Fortran may multiply the number in the data 
by 104 if the number before the — is integral / D 
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>) ae (e) re D FLOATING -POINT FORMAT FOR 
W.d@ POUBLE-PRECISION VALUES 


he item in the 1/0 list corresponding to this descriptor must be of 
type DOUBLE PRECISION. ( Some Fortrans make £ and dD descriptors 
interchangeable but the standard is more restrictive . D 


DOUBLE PRECISION DBL 
READ (5,100) DBL 
100||\FORMAT ( D12.0) 





VAROTONS ° The binary value in the computer is converted to 
character form as follows: 


- digits Dt ee 

- digits D-ee 
where digits represents the d digits specified; ee represents a two- 
digit exponent. Values are rounded to @ decimal places and the 
resulting number is right justified in a field of w character 
positions. 


SPhe variations in style of output described on the previous page 
@in particular the replacement of &£ by D and vice versa ) may apply 
to the D descriptor as well. These variations are all admissible in 
the Fortran 66 standard. Again, if you make w greater than d+6 
the field should be wide enough to contain the result whatever the 
style of output. 


DOUBLE PRECISIGN DP(4) 
DATA DP(1), DP(2), DP (3), DP(4) 
/ 76-573D0 , -587963.66D0, 37.5D-2, 68D-4 / 


WRITE (6,200) (DP(I), I= ft, 4) 
FORMAT (__ 1X, D12.6_) 
———————— 












-765730D+02 
-—+-587964D+t05 
*375000D+00 
*680000D-02 









o0o9g9 
9000 





A 


colt 
CARRIAGE 
CONTROL 






RSADACYA)? The D descriptor behaves as the F and £ descriptors, 
but remember the corresponding item in the 1/0 list must be of type 
DOUBLE PRECISION. 
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PKA Gua Ceres? formar 


She item in the 1/0 list corresponding to a 6G descriptor should be 
of type REAL. The item may also be of type comPiéx provided that 
if corresponds to a pair of such descriptors. 


REAL R 
COMPLEX Cc 

READ (5, 100) R,C 
FORMAT (G10.0, 268.0) 














Q 


° The binary value in the computer is converted to one of 
two character forms: 


@ if the absolute value lies between 0.1 and 
10% then in fixed-point form with four blanks 
to the right @ reducing the effective width of 
field to w-4 character positions ). The best 
use of the remaining field is made as illus- 
trated below 


@ otherwise in &-format right justified in the 
full field width of w character positions. 


REAL A(#) 
DATA A(1), A(2), A(3), A(&) 

/ 0.0666, 0.666, 6,66, 66.6 / 
DATA A(5) , A(6) , A(?) 

/ 666.0, 6660-0, 66600-0 / 
WRITE (6, 200) (ACI), T= 1,7) 


»-6660E-01 
0.6660 
6-660 
66.60 
666-0 
6660. 
-6660E +05 


o 
oO 
Oo 
oO 
© 
o 
oO 


CARRIAGE 
CONTROL 





RISAVDOLY A ¢ The G descriptor behaves as the F,£ and D descriptors. 
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P IN CONJUNCTION WITH F,E,D,G 
i DESCRIPTORS € DANGEROUS /) 


The nm is an integer constant which may be negative, zero or 
positive. The scale factor is 10”. 


° In F-format the number printed is 10” times the 


value _ stored: 


DATA B/1.5/ 
WRITE (6, 100) B, B 
100] |FORMAT (1X, 3PF8.2, -2PF8.2) 





Bn E-format and p-format the decimal point is shifted » places 
to the right @ when n7 is positive ) and the exponent is reduced by 
n. In other words the number printed is sot scaled at all but 
simply altered in appearance : 


9/-1.5/ 
WRITE (6, 200) 9,9 
FORMAT (1X, £8.2 












1PE8.2 ) 












© 








CARRIAGE 
CONTROL 


Qn G-format the scale factor has no effect on numbers printed in 
the style of F-format @ range 0.1 to 104 }. For numbers outside this 
range the effect of the scale factor is the same as that just 
described for the £- format. 


DIVA 2 When you write an item of data with no exponent 
then the value of that item is divided by 10”. When you write an 
item of data with an exponent then the scale factor has no 
effect on the value stored: 


READ (5, 300) (A(Z), I= 1,3) 
300) |FORMAT (2PF6.2, F6.2, F6.2) 
AF 


















Whe value of 7 is assumed to be zero when not specified at all 
@ 10°=1) but once a value for 7 is specified it remains in force 
for all other descriptors in the same FORMAT’ statement € even for 
descriptors in deeper nested brackets and in all rescans ) until 


re-set by another nP. 
ZERO NULLIFIES THE 
Jar 
FORMAT (2PF6.2, OPF6.2 ) 
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I ONLY INTEGERS MAY BE TRANS- 
Ww MITTED USING THIS DESCRIPTOR 


° Integer values to be printed are right justified in a 
field of w character positions: 


WRITE (6, 100) K,K 
ORMAT (1X, 2123) 





READING = Items are assumed to be right justified within a field 
width of w = blanks being freated as zeros. Negative values 
Shovid be preceded by a minus sign. Plus signs in front of positive 
values are allowable but not necessary. 


INTEGER K, L, mM 
READ (5, 200) K, L, mM 
FORMAT (3F3) 


22-2 2 ee ae 





Rial |i? ONLY LOGICAL VALVES MAY BE 
w TRANSMITTED WITH THIS DESCRIPTOR 


WRAL ° In a field of w character positions the letter 7 or F 
is right justified. 7 signifies the Boolean value ¢rve and F signifies 
false . 


M,N /.7RUE. , .FALSE. / 
WRITE (6, 300) M,W 
FORMAT (1x, Lt, L2) 


CARRIAGE 
CONTROL 


PUA * Inside a field of w character positions the first 
encounter of the letter 7 or F signifies the Boolean value frue or false 
respectively. The Boolean value is assigned to the corresponding 
logical variable nominated in the I/O list: 


LOGICAL M,N 
READ (5,400) M,N 
400\ |FORMAT (22410) 
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Sirol A FOR DEALING WITH 
Ww CHARACTERS 


Onlike Fortran 77, Fortran 66 has no variable of type CHARACTER . 
Instead characters may be stored in variables of any type = and it 
depends on the computer how many characters each type of variable 
can hold. On oa typical 16-bit computer an integer variable may 
hold two characters, a real variable four, a double-precision or complex 
variable eight. For the sake of generality we speak of a variable 
@hence also an array element } holding g characters in the illustrations 
below. In the examples g is assumed fo be 4, and a blank @or 
space ) is depicted as ‘ 


WRIA @ If the field width w is greater than or equal to g then 
characters are right justified in the field and enough blanks added 


to fill the field: 


DATA V(t), W(2), W3) , V4), VS) 
/ GHFOUR, 1H , 4HCHAR, 4HACTE, 2HRS / 


WRITE (6,100) (v(I), T= 1,5) 
200| |\FORMAT (1X, I6, 4I4) 






















If the field width is less than g then the left-most characters ore sent 
to fill the field, those to the right being lost: 


WRITE (6,200) (V(t), IL=1,5) 
200) |\FORMAT (1X, 512) 


year 






CARRIAGE 
CONTROL 







RRISADOIG: If the field width w is greater than g then only the 
right-most g Characters are stored in the associated variable, the 
foremost w-g characters being skipped over: 


READ (5,300) (CV(I), I= f,3) 
300} |FORMAT (378) 


FOUR CHARACTERS PER LOCATION 





If the field width is less than or equal to g then all w characters 
Ore left justified in the variable, the remainder of the variable 
being filled with blanks. 


READ (5,400) (v(Z), I= 1, 8) 
400||FORMAT (813) 


FOUR CHARACTERS PER LOCATION 
vay [Foo 


CBC! — vin [RIBBTY7 
t t F fi} 
a ' 0 





















HOLE RUA HT LITERAL wHlhh...h — CAprions 


WFNS 2 The w characters to the right of letter 4 define a 
field of width w and are sent to fill that field precisely. Blanks 
within this field are significant characters. The set of characters 
that may be used in a Hollerith literal is more extensive than that 
defined as the Fortran 66 character set and more extensive than the 
Fortran 77 set @ page 20) but for portable programs it is best not to 
use non- standard characters. 


WRITE (6, 100) 
100| \FORMAT (1X, 26H THIS TS A HOLLERITH LITERAL) 


© THIS IS A HOLLERITH LITERAL oO 


AF CARRIAGE CONTROL 


ADONS : According to the Fortran 66 standard the w characters 
of the Hollerith field on the input medium should replace the 
corresponding w characters appearing in the formar statement as 
illustrated here: 


READ (5, 200) 


200| |FORMAT (1X, GH«xxx ) 
WRITE (6, 200) 


© 
© 
o 
oe 
o 





ABCD 


Se008 


a ED 


Although this is standard Fortran 66 the practice is not recommended 
to those wanting to write fully portable programs. There are some 
q otherwise standard ) Fortrans that do not offer this facility of 
replacement. Try it on your own machine. 


4 FOR PRINTING SPACES 
Ww OR SKIPPING DATA 
WPI: A total of w blanks are sent to the output record 


on each wX encountered. See page 101 about carriage control. 


WRITE (6, 300) 
FORMAT (1X, 3HHOW, 1X, 3HNOW, 1X, SHBROWN, 3X, 3HCOW) 


O ,HOW NOW BROWN cow o 


Wangan cornea) 


°@ A total of w characters @including blanks ) on the 
input record are skipped over on each wX encountered: 


READ (5, 400) J, K 
400| |FORMAT (4X, Iqy, 8X, I3) 
NON 1234 SENSE 2468 AND MORE 


oP An ft 14 J 3 
















THE 
WORDS ON 
THE CARD 
ARE 
SKIPPED 
OVER 





x [Lia] 
k[ 468 | 










q 
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FORMAN: AN EXAMPLE TO ILLUSTRATE THE 
AVO/DANCE OF TROUBLESOME DESCRIPTORS 


Anput in Fortran is based on the concept of the punched card but nowadays 
very many (Cif not most) computer installations offering Fortran do not use 
punched-card equipment. It can be awkward trying to ensure that 
numbers typed at o VDU or on punched tape = appear in the correct 
fields. The solution to this problem adopted at many installations is to 
provide a special format descriptor or special keAD statement which 
simply reads the next number whatever field it occupies. The disadvantage 
of this solution is that it brings non-standard features into Fortran «= and 
these non-standard features vary in specification from one computer to 
Qnother. 


‘Phe subrovtine described here illustrates a way of reading numbers 
in free format. The subroutine 1's nevertheless written in standard 
fortran 66 and shovid be fully portable. 


Qn each call the subroutine reads the next waiting number as a REAL. 
The waiting number may be prefixed with a plus or minus sign and 

may contain a decimal point. This subroutine does not allow numbers 
written in exponent form, nor does it permit a leading or trailing 


decimal point: 
- 1.2783 ae --+5 aH. 36. 


Numbers should be written separated by one or more spaces as follows: 
-1270 -0.5 38.0 Gor just 38 D 


She subroutine @ which musi be initialized: see note opposite may be 
called as illustrated below: 


LOGICAL OK 


REAL V 
CALL INPUT(YV, OK) 


if ox returns true the item that was waiting is returned by Vv. If 
ox returns false it means the waiting item was wrongly formed. 
However, the pointer is then restored to where it was before the call « 
thus making it possible to attack the item with another subroutine € for 
instance a routine for reading keywords ). 





She subroutine weur seeks the next item. This subroutine reads a 
new card (€or card image) if there are no more items in the 
current image. 


fhe logic of this subroutine is based on the following state table: 
SYMBOL 


: : id : : : END OF 
ree 
ACTION 1: ACTION 3: ICTION Gi 
pt ACCUMULATE | NEXT STATE=2 NEGATE THE | exir raLse \weaT STATE =1 | REAP fen. 
NEXT STATE= NEXT STATE =2 AT STATE <1 
ACTION {: 
2 1 eaiahel EXIT FALSE | ExiT FALSE | EXIT FALSE | EXIT FALSE | ExiT FALSE 
NEXT STATE=3 
ACTION f: 
NEXT STATE=3 
ACTION 2: 
4 aS i ExiT FALSE | EXIT FALSE | EXIT FALSE | EXIT FALSE | EXIT FALSE 
NEXT STATE = 5 
ACCUMULATE | ExiT FALSE | EXIT FALSE | EXIT FALSE \ EXIT TRUE | EX/T TRUE 
FRACTION « 
NEXT STATE= 5 
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SUBROUTINE INPUT (VALUE, STATUS) 
LOGICAL STATUS, SAME 


INTEGER ROW(80), L(%4),7(5,6), READER, SYMBGL, STATE 


ENTRY, K, ACTION, FIRST , NEXT 

REAL VALUE, DIVISR, SIGN, DIGIT sey 

COMMON / POINT / NEXT we 

DATA READER / 5/ 

DATA L(1), L(2), £(3), L(4), £05), LC), 4(7), (8), £09), L (10) 
J {HO, tht, 1H2, 1H3, 14, 1H5, 1H6, 147, 148, 149 / 

DATA L(t1), L£(12), £(13), L014) a 

/ tht, tH-, tH. , tH 

DATA T(1,4) 5 74,2), THL,3), TC4,4), T019S)), TCLy6) 

















/ 3, 72, 32, 60, Wl; 48 / 
DATA T(2,1) , 7(2,2), T(2,3), 7(2,4), T(2, 5), 7 (2,6) 
13, 60, 60, 60, 60, 60 / 
DATA 7(3,1) , 7(3,2), 73,3), 7(3,4), T(3, 5), T(3,6) 
/ £3. 5 60, 60, ee, 50, 50 / 
DATA T(4,1) 5 74,2), 7(4,3), T(4,4), T4,5), 74,6) 
60 


P ; 60, 60, 60, 60 / 
DATA (5,1), 75,2), T05,3), T(5,4), 7(5,5), 75,6) 















/ 25 , 60, 60, 60, 50 , 50 / 
VALUE = 0. WF sympee 
DIVISR = 1.0 TABLE 
STATUS = .FALSE. ‘NEXT’ MUST BE INITIALIZED BEFORE 
FIRST = NEXT Ph a FIRST CALL TO ‘INPUT’ 

= e.g. 





BLOCK DATA 
COMMON / POINT / NEXT 
DATA ExT /81/ 

ND 









F (SME(K, L(I))) G6 76 20 
ONTINUE RE ee 

CHARAC WO 
co 8's > Cie) 


IF (LZ -GT. 10) SYMBOL = I-9 
ENTRY = T(STATE , SYMBOL) 
= ENTRY / 10 
STATE = ENTRY — LOX ACTION 
GO TO (1,2,3,4, 5,6,7), ACTION 

























VALUE = 10.0* VALUE + DIGIT 
7 





10.0% VALUE + DIGIT 
10-0 DIVISR 


READ (READER, 100) ROW 
(80A1) 


SIGN & VALUE /DIVISR 


NEXT = FIRST <i] FALSE EXIT 
RETURN (POINTER RESTORED) 


END 
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RCASIES CHAPTER 10 


10.0 1 Work through some of your earlier programs replacing the 
rudimentary output with neatly formatted results. In particular 
let input values feature among the results. For example the 
results from Exercise 9.3 should appear something like this: 


© ENCODING ROMAN NUMERALS 


1492 = MCDXCII 


0 
fe] 
© 





10.2 Write Q program to plot several functions simultaneously. For 
example y=sin(x), y=cos(x) and y<tan(x) using a different 
symbol for points on each curve: for example * and + and. 


79.3 [$xtend subroutine sNPUT @page 115) to cope with numbers 
written in exponent form. This involves an extension to the 
state table. 


10.4 U/rite a subroutine to read a “keyword” and match it 
against a list of allowable keywords stored in the program. 
For example let the list be RECTANGLE, TRIANGLE, CIRCLE and let 
the first ¢wo characters suffice for a match: RE, 7R, cr. If 
the subroutine discovers the word RECTANGLE @ or RECT OFT RE or 
even REk ) it should return 1. If the routine finds the word 
TRIANG it Should return 2, and so on. If the routine fails 
to find a match in the list for the item just read then the 
routine shovid return zero €or set a logical argument fase > 
and set its pointer back where it was when the routine was 
called ¢2.e. follow the principle used in subroutine sMPuT ). 


Anclude this subroutine = together with subroutine INPUT *» in 
fhe example program on page 45. The data for this program 
could be recast as follows: 


TRIANGLE 15.4 16.8 24.95 
RECTANGLE 73-67 10 
TR 3 4 5 


CIRC 15.9 ASSUMING “NO” /$ ADDED 
NOM@RE TO THE LIST OF KEYWORDS 
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FOIUSS 


FORMATTED FILES 
UNFORMATTED FILES 
ENDFILE, REWIND, BACKSPACE 
EXERCISES 


SOME CONCEPTS 
Le) AND TERMINOLOGY 


LA formatted file consists of physical records: each 1/o list speci- 
fies a logical record. 


ZA deck of cards punched in readable @character ) form is a 
formatted file: 






READ (KARDS, 100) (V(I),Z>1, 48) 
100||\FORMAT (._8F10.0  ) 

——_ 

EACH RESCAN 


BEGINS THE NEXT RECORD 

The length of each record is limited by the external medium € in this 
case the punched card ). It would be a mistake to write the format 
as (48F10.0) because one eighty-column card cannot contain 48 fields 
each of 10 columns. The FormAT statement is responsible for chopping 
up the I/0 list = ¢.e. the logical record = into manageable physical 
records by rescans or slashes or both. 





“A set of printed output is also a formatted file: 


WRITE (LINPRN, 200) (V(I), T= 1, 48) 
FORMAT (_1X, 6F 10.1, ) 
S= 


EACH RESCAN 
BEGINS THE NEXT RECORD 


30.0 40.0 
90-0 





eo09000000ee99 
eoeea0ce8 0 e0 





The length of each record is again limited by the external medium ¢ in 
this case the longest line the particular printer can print ). It would 
be a mistake to write the format as (1x,8Ffo.1) if the printer had a 
line length of 72 columns. Again the Formar statement is responsible 
for chopping the 1/0 list into manageable physical records by rescans 
or slashes or both. 


A formatted file may also reside on punched poper tape or magnetic 
tape or on a disk file. In such cases the I/o list in conjunction 
with the FoRMAT statement again determines the length of record = but 
the external medium does not impose any practical constraint on 
length. @ The computer does divide magnetic tape into “blocks” but this 
should not have fo concern the Fortran programmer directly.) In the 
following example the record length is chosen as eighty columns « in 
other words as a “card image”. 










INTEGER LINE (80) 
26 10 I=1,6 


i, 7 READ (KARDS, 300) LINE 
10\|WRITE(ITAPE , 300) LINE a7 
FORMAT (80AL) 





FORMATTED 
FILE 


Beware of writing records in one format and later trying to read 
them back in another. This is because many operating systems do 
not “pad out” short records with blanks to represent card images. 
Instead they keep account of the number of characters actually in 
the record. Below is shown an attempt to read twelve characters where 
only ten were written: all right with cards but not with disk files. 


WRITE (IFILE, 400) X ala READ (IFILE, 500) Y 
4OO|\|FORMAT ( F10-2) one 500} |FORMAT ( F 12.0) ~sp 
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NRO (FYES thot reeminoiocr 


“An unformatted file consists of logical records, each as long as the 
I/o list demands. @Some operating systems do impose a limit on length but 
it is usually quite big- 2 An unformatted file cannot be sent to a line 
printer; if it could the output would be unreadable because unformatted 
files are coded in binary form. 


A formatted file may be converted to an unformatted file on a_ suitable 
medium @say magnetic tape ) by means of a program: 
BINARY G4®S 


ll DIMENSION LINE (80) 
i! DO 10 I= 1, 6 “CARD | 
C READ (KARDS, 100) LINE ia e 
7 {0| \WRITE (ITAPE) LINE [o> | A g iy 
ee Lalit ¢ S0A1 ) UNFORMATTED FILE 


In the above example each unformatted record is the binary equivalent 
of eighty integers » a “binary card image” in which each integer 
stores one Hollerith character in its most significant bits. 



















FORMATTED FILE 


G@onversely an unformatted file may be converted to a formatted file 
by a program. To write this you have to know: 


@ the number of records in the unformatted file 
@ the number of items in each logical record 
(they could all be of different length / ) 

@ the type of each such item @ these could all 

be different too ) 
@ the limiting length of formatted record to be 
written 


Assuming the file of six logical records illustrated above « and assuming 
the data listed Opposite « the program could be as follows: 


Ps = £5 6 
((“ READ (ITAPE ) LINE 
—EE===_ oo 








© 10:0 200 30:0 40:0 2 
© 60.0 60:0 70.0 §0.0 ° 
° 90-0 . ! n.0 oO 


FORMATTED F/LE 
(OF TWELVE LINES ) 
be 5 he o 
© 410.0 420.0 430-0 4400 © 
© 430-0 460:0 4700 480-0 © 


Finally we illustrate a straightforward application of unformatted files 
used as “backing store” in the “number crunching” type of problem. 
A, B, € may be big multi-dimensional arrays accounting for very long 
records. 








READ (IFILE) A 
READ (JFILE) B 


\CALL CRUNCH( A, B, C) 
WRITE (KFILE) € 





PO you may do more than just read and write files serially; Fortran 66 
provides £NDF/LE, REWIND and BACKSPACE @ See overleaf ). These offer only 
limited control of files but the deficiencies in fortran 66 are overcome in 
Fortran 77. Here you may detect an end-of-file ; make access to any record 
of a file directly; inquire into the status of a file @e.g. whether it exists ) 
and take action if there is a misread. Consult your own Fortran user’s 
manual and check if any such facilities offered conform to Fortran 77 standards. 
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BND ZOE , REWIND. BACKSRASES Srarctens” 


As their individual names suggest, the auxiliary 1/0 statements were origin- 
ally designed for manipulating magnetic tapes. They are equally applicable 
to manipulating files stored on a disk or other modern secondary - storage 
device. But the auxiliary 1/0 statements would make no sense in controlling 
the more traditional devices: you can’t “rewind” or “backspace” a card 
punch, Cord reader, paper-tape punch or reader, or line printer. 


NOFA puts a dab of conceptual paint on a magnetic tape. This is called 
an end-of-file record (EOF). In bygone days the EOF was useful when 

the reel of tape was taken off the computer’s tape deck and mounted ona 
special deck for driving a line printer ¢“ off-line printing” in the jargon). The 
EOF ~ the dab of conceptual paint = was recognized by the printing installation 
which would throw a few blank pages before printing the next file. But 

Fortran 66 does not recognize an EOF record when reading a file. If an EOF is 
encountered during input different Fortrans do different things: some give up 
altogether. 







DATA MyYFILE (#/ 
WRITE (MYFILE) (VWI), I= 1, 800) 
ENDFILE MYFILE 





KECZIND positions a conceptual pointer at the first record of the nominated 
file. If this file happens to be already rewound then REwinp acts as a 
“do nothing * statement. @ There is danger using REwivd on magnetic 
tapes: some Fortrans rewind the entire tape rather than just one file. 
With disk files, however, the RE&wiwD statement should behave properly. ) 





WASIMOIAN SI: causes the conceptual pointer to move back over the 
previous record @ E0F record no exception ) of the nominated file. If this 
file happens to be rewound @€ the pointer at the first record ) then BSACKSPACE 


acts as a “do nothing” statement. 
RECORD 
a 















“Do NOTHING” 


Because of the different behaviour of different Fortrans it is best not to 


use BAcKSPACE with unformatted files. Here is an example using a formatted 
file. 


















DATA KFILE /[8/ 
DO 10 K= 1, 601, 200 
L=kK*+ 200-1 

10||WRITE (KFILE, 100) (A(t), I=K,L) 
100||FORMAT (200F 10.2) 


> 
| | {BacKsPAce — KFELE ies, THe EOF ecORD 

Se ow BACKSPACE | 
| BACKSPACE — KFILE 






D 2ND 
ru 3k ts7 
RECORD RECORD 




















SECOND 
BACKSPACE 


120 


Whe forms of the auxiliary input and output statements are as follows: 


ENDFILE unit 
REWIND unit 
BACKSPACE uait 


where: 

unit is either an integer constant @ digits only ) or the 
name of an integer variable @not an array element ). 
unit denotes a connection to a magnetic-tape deck 
or disk file. The connection is made by a 
command to the computer’s operating system ~ 
perhaps by Q PROGRAM command placed immediately 
before the main program. 


Phese are erecvtable statements and therefore may be labelled. 


[For the following example assume the installation shown. Units 5 and 6 
are connected fo the card reader and line printer according to the usual 
convention. Unit 4 is connected to a magnetic-tape deck or disk file. 
The object is: 


@ to read the waiting cards to a file on 
magnetic tape or disk 
G @ print the records of the file in reverse order. 





Whe end of the card deck is signified by a blank card. All other 
cards have non-zero integers right justified to columns 10 and 20. 


INTEGER READER, PRINTR, TAPDEK, NUMBER, I, M,N 
DATA READER, PRINTR, TAPDEK [ 5,6, 4 | 


REWIND TAPDEK 

NUMBER = -1 
NUMBER = NUMBER + 1 

READ (READER, 100) N,M 

IF ((N .£Q.0) .AND. (M.£Q:0)) GO TY 20 

WRITE (TAPDEK, 100) N,M 

GO TO 10 


IF (WUMBER .£0.0) STOP f£ 

ENDFILE TAPDEK 
DO 30 I= £, NUMBER 
BACKSPACE TAPDEK FIRST TIME ROUND 
BACKSPACE  TAPDEK THe EOF RECORD 
READ (TAPDEK, 100) N,M 
WRITE (PRINTR, 200) N,M 
CONTINUE 


TAP DEK 


(2Z10) 
(1X, 2IZ10) 
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23K SISIES CHAPTER 1 


14.1 Develop Q@ program for writing card images (either from punched 


cards or from lines of data fyped at a wu) to a file on 
magnetic tape or disk. Do this for writing both formatted and 
unformatted files and compare the sizes of files created from 
the same card images. Most, if not all, installations allow users 
to discover the sizes of files they have created. 


11.2 Develop a program for printing the contents of a magnetic -tape 


or disk file of card images. Your Fortran probably provides a 
means of detecting the EOF record. Use the facility provided 
but do discover whether or not it conforms to the Fortran 77 

standard. 


44.3 “Amend one or two of your favourite programs to read from magnetic- 
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tape or disk file = and write results to another magnetic-tape 
or disk file. Use the programs specified in 11.1 and 11.2 above 
to provide the data and print results. 


2 


VOR 
AXA APES 


LINEAR S/MULTANEOUS EQUATIONS 
SHORTEST ROUTE THROUGH A NETWORK 
REVERSE POLISH NOTATION 
EXERCISES 


DYNA SOO GRAN SOS EeMOAuie8 


(Somputers spend many hours solving sets of linear simultaneous equations. 
They crop up in engineering @ stress analysis of bridges, buildings, aircraft) 
and many other fields. The solution of a set of several thousand equations takes 
much ingenuity to devise because intermediate results have to be filed on disk 
and rapidly recalled. And there is the problem of accuracy; solution in a 
badly chosen order can yield useless answers or none at all. The example 
here illustrates just the bare bones of the problem. 


Gonsider these three simultaneous equations ¢ called linear because there is 
no z?, y3 etc.): 
sx + Oy + 5Z 


oe 2 
Feecall that we are allowed to 12x * My + BZ = 7.3 
multiply any equation all the way 6x + 36Z = 3.5 


through by a non-zero multiplying 

factor; also to subtract one equation term by term from another so as to 
replace either of these equations. Neither of these operations changes the 
solution unless we start getting small differences between large quantities 
hence introducing inaccuracy. The atm of these operations @ demonstrated 
below) is to create a new set of equations looking like this: 


@nce the equations are reduced 3x 4+ fOyY + SZ = G2 
to this triangulated form it is 16y + 4Z = 3.9% 
easy to solve them. Starting with 35z% = 2.805 


Zt: Z= 2.805+35 = 0.0801 . 

Knowing the value of z substitute in the second equation to find y: 
16y + 4x 0.0801 = 3.94 hence y = (3-94 - 4x 0.0801) + 16 = 0.2262 
Knowing the values of y and z _ substitute in the first equation to 
find xz: 5x + 10x0.2262 + 5x0.0801 = G2 hence x=(4.2-— 10x 0.2262 
- 5x 0.0801)+ 15 = 0.1025. In summary x= 0.1025, y= 0.2262, 
Z= 0.0801. This is called back substitution. But first we must 
triangulate. 


“fo get rid of the coefficient of x in the second equation subtract 
a multiple of the first equation term by term from the second = and 


so create Q new second equation. To do this the multiplying factor 
must obviously be /2+/5: 
2° equation: 12x + 24 + @z* 


y 7.3 
minus + (12415) x15x% + (12415)x 10y + (1271 X § Z = (12415) * 4:2 


new 2% eg. : Ox + 16y +4Z * 3.94 


Wo get rid of the coefficient of x in the third equation subtract a 

multiple of the first equation term by term from the third. This time the 

multiplying factor must be 6+/5: 

37 equation : 6x + 36Z = 3.5 
minus : (6415) x15 x # (64/5 )x Loy + (6+15)x5z = (6215) x G2 


new 37 eg, : Ox -4y 4$34Z = 1.82 


VeVe have now finished eliminating coefficients of x € and incidentally 
have finished using multiples of the first equation which remains 
unchanged ). The equations now look like this: 


GSH +#+10Y + SZ = 4-2 
16y + 4Z = 3-94 
- 4Y 434Z = 1£.82 
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‘yho get rid of the coefficient of y in the third equation subtract a 
multiple of the second equation from the third. The factor is -4+16. 





3° equation: -4y + 34Z = 1.82 
minus : (4+ 16)x16y + (-4Fl6)x GZ = (-4+ 16) x 3.94 
new 3 eg, : oy + 35Z = 2.805 


At last the equations are triangulated. Notice the second equation remains 
unchanged from above: 


15x + 10y ¢t 52 = 4.2 
and the solution may be l6y + 42 = 3.94 
found by back substitution 35Z = 2.805 


as already demonstrated. 


SPhis method works as long as you don’t divide by zero when forming the 
multiplying factor « and gives good results when coefficients on the 
diagonal @/5, 24,36 opposite ) are large compared to those off the diagonal. 


She following subroutine operates on array A(,) in which the numerical 
coefficients of the equations are stored « and vector B() in which the 
right-hand side is stored. The third parameter, WUMBER, specifies the 
actual number of equations fo be solved. The subroutine changes array 
A(,) to triangulated form and replaces the contents of vector 8() with 
the values of the “unknowns”. 


SUBROUTINE EQUNS (A, 8, NUMBER) 

INTEGER NUMBER, PENULT, I, J,K, NEXT 

REAL A(NUMBER, NUMBER), B(NUMBER), DIVISR, FACTOR 
Ir 


INATE COEFFICIENTS “ADJUSTABLE” DIMENS/ONS 
PENULT = NUMBER - 1 gg Codey 
Dé 20 I= tf, PENULT 
DIVISR = A(I,I) 


IF (ABS(DIVISR).LT. 0.000001) STOP 1 <3 
NEXT = I+ 1 
Dd 30 J = NEXT, NUMBER 
FACTOR = A(J,I)/ DIVISR 
DO 40 K= 1, NUMBER 
( A(J, K) =A(J,K) — FACTOR* A(I, K) 
CONTINUE 
BOT) = BIT) -— FACTOR *B(I) 


CONTINUE 
CONTINUE 


SUBSTITUTE 


D@ 50 I= 1, NUMBER SE T couwrs “up” 
IF ( K .&Q. NUMBER) GOTO FO 
NEAT = K#1 
DO 60 J= NEXT, NUMBER 
(2 = B(K) — B(J)* ACK, J) 
CONTINUE 
B(K) = B(K)/A(K,«K) 
CONTINUE 
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[Finding the shortest Gor longest ) route through a network is a problem 
that crops up in various disciplines — one of which is critical path 

analysis for the control and monitoring of construction projects. Given a 
network such as that below the problem is to find the shortest route 

from the node marked s7arrT to that marked F/w/sH. The journey must 
at all times be in the direction of the arrow. The number against each 
arrow shows the journey time for that arrow. 









ALL SUBSEQUENT 
RECORDS 


From WTO JOURNEY 
NODE ,» NODE TIME 
IN ANY ORDER 


Whe program starts by reading the first record: >) 







INTEGER NODES, ARROWS, START, FINISH, ROUTE (50), 
HEAD (50), LINK(20), TIP (120), KOUNT 

REAL TIME (120) , BESTIM(50), BETTER 

OGICAL SWITCH (50), ON, OFF 


\AS ALLOW 50 NODES 
ZZ1 AND 120 ARROWS 












DATA ON, OFF / .TRUE. , FALSE. / 





READ (5, 100) NODES, ARROWS, START, FINISH 
100| |FORMAT (425) 


“Phen four entities are established for each node as illustrated below: 


ot swirew(i){ ON] Bestin(i) [1000000.9] HEADCi) | O 4 roureci [OJ 
k ene) GR) 
where the switch @ explained later) is set to ov; the best time to this node is 
set impossibly high except for the s7agr node for which the best fime is 


invariably zero. The head of a chain linking all nodes running out of this 
node is set to zero; so is the link on the chain to be created along the best rovle, 


-? 
mee * 


DO 10 TF=1, NODES 
SWITCH (I) = ON 
BESTIM (2) = £6 


HEAD(I) = O 
ROUTE(I) = O 
CONTINUE 
BESTIM(START) = O 





— 


Now the remaining data are read. The node at the tip of each arrow is 
read into the vector named 7/P(); the journey time is read into the vector 
named 7/ME(); the node at the tail of each arrow is linked into a chain 
held in the vector z/vK() and with its head in vector wWEeAD(). 


DQ 20 J= 1, ARROWS 
READ (5,200) WN, TIPCJ), TIME(S) 
LINK (J) = HEAD(N) 


HEAD(N) = J 
CONTINUE 
FORMAT (225, F410.0) 
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he structure of data stored in the computer is now as depicted below: 





LINK TIP TIME 
SWITCH BEST/M ROUTE gz 
(1) [on | za 
(2) {ON | | 2 | 
START o> (3) |_ON | ER 
¢ 2 oe 
| S| 
rousuci>Ce) For] re] 
- 
Ep ae a 


Row for the working Perl of the program. Start at node s7ARkT and visit each 
node in turn: 3,4,5,6,.... going round again...,1,2,3,4,... until all switches are off. 


K = START 
GO TO 40 o> © 


At each node consult the switch: if off move to the next node: if on then run 
down the chain of arrows out of this node. For each arrow add its journey time to 
the best time so far achieved at the tail « and see if this sum gives a better 
route to fhe tip. If so switch on at the arrow’s tip and replace the previous 
best time with the new one. Also put the node number at the tail into the 
vector ROUTE() so as to build a chain along the best route. Having completed 
all this work at one node switch off that node. Keep variable xKounr for 
counting the number of cycles through the nodes = but reset to zero every time 
you meet a node which is 07. So when Kounr7 reaches 2 all switches are off. 












K=kK#tI 


IF (kK .LE. NODES) GO TO 35 

ke<f ue 
KOUNT = KOuNT +1 

IF (KOUNT .GT. 1) GO TO 85 

IF (.NOT. SWITCH(K)) Gt TR 


® Cee 49 |Kounr = 0 

N= HEAD (K) <EW HEAD oF {nEAD OF CHAIN, FOR NODE K ) 

IF (NW .£Q. 0) GO TO 80 ) 
BETTER = BESTIME(K) + TIME(N) 


L= TIP(N) 
IF ( BESTIM(L) .LT. BETTER) GO TO 70 


SWITCH (L) = ON 

BESTIM (L) = BETTER 

ROUTE (L) = 

N= LINK(N) 

GO 76 50 

SWITCH (K) = QFF Pe 


G® TO 30 











ALL SWITCHES 
ARE OFF 















WRITE (6, 300) START, FINISH, BESTIM( FINISH) 

FORMAT (16H BEST ROUTE FROM, Ig, 3H 70, TG, 6H TAKES, 
F6-1, 8H THROUGH ) 

NM = FINISH 

IF (N -.£Q. 0) STOP 

WRITE (6,400) N 


’ 


GO TO 90 















FORMAT (1X, I4&) 
END 







©sing the data opposite, this program should say the best route from 3 to 6 
takes 31.0 (units of time) through nodes 6,5,4, 2, 1,3. 


127 


RSS RIS] INO WAUOIN tse or sraces 


Aiigebraic expressions may be converted from the conventional form to a 
form without any parentheses. This notation is called reverse Polish. For 
example: 

At(B-C)*D-F/(GtH) becomes ABC-D¥+FGH +/- 


‘rhe new expression is easier to evaluate than might appear. for example 
let A=6.0, B<4.0, C=1.0, Dx2.0, F=3.0, G=7.0, H=5.0. With these 
values the expression to be evaluated is: 

6.0 4.0 10- 2.0 x + 3.0 7.0 5.0 + f - 
Work from left fo right pointing to each item in turn. Whenever you come to 
an operator apply that operator to the previous pair of terms. This 
procedure is illustrated below in steps: 

6.0 4:0 1.0 — <doreraror 

6.0 3.0 2.0 * <KhoPERATOR 

6.0 6.0 + <Qyoreraror 


12.0 3.0 7.0 5.0 + <khoreraror 
1N OTHER WORDS 
6.04(b-0= 1.0) ¥ 2.0 = 3.0/(7.0 45.0) 12.0 3.0 12.0 J SBaoreraror 
= 11.75 12.0 0.25 — ey OPERATOR 


11-75 <Kqresuc7 


‘he logic for converting an expression into reverse Polish notation is given 
below as a flow chart. Two stacks are used ond referred to as STACK Xx 
and stack Y. Operators have precedence: * and / have the highest 
(no exponentiation in this example ) + and - come next; finally space is 
considered an operator with lowest precedence. A space is used to terminate 
the expression. 


READ A RECORD OF 80 CHARACTERS 
SET POINTER TO ZERO 
(10) 


ADVANCE POINTER TO GET NEXT CHARACTER; C aaa: 
LETTER? » PUSH CON STACK X 
iy. PUSH C ON STACK Y 














IS PRECEDENCE OF ¢ 
GREATER THAN THAT OF X? 


PUSH K ON STACK X 








PUSH K ON STACK Y 









PUSH CON STACK Y 






NO 






PRINT CONTENTS 
OF STACK X 
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Phe following program assumes the existence of subroutines PUSHxX, POPX, PUSHY, 
Pory for managing two integer stacks. These may be modelled on PusH(EXPRN) 

and POP(TroP,ok) on page 79 but also declare common/sracxs/ as in the main pro- 
gram. The program also refers to a function subprogram named prEceD(r) for 
returning the precedence of an operator *# using values given to characters by 

function /WDEX(W) on page 89: 


INTEGER FUNCTION PRECED(Z) 

INTEGER DIV, MULT, MINUS, PLUS, SPACE 

DATA DIV, MULT, MINUS, PLUS, SPACE / 42, 41,40, 39,37 / 
IF ((1 -£Q. DIV) .OR, (I .EQ. MULT)) PRECED = 3 


IF (C1 .£Q. PLUS). OR. (I -£Q. MINUS )) PRECED = 2 
IF (ZI -£Q. SPACE) PRECED = 1 

RETURN 

END 





Were is the program itself: it reads a record containing an expression composed 
of letters, operators and brackets = then prints the expression in reverse Polish 
notation. The program returns to read and decode a new expression starting 
in column 1 but stops if it meets a dollar sign. 











LOGICAL OK 
INTEGER CARD(60), POINTR,C,K, LEFT, RIGHT, SPACE, PLUS, MINUS, 


i MULT, DIV, DOLLAR, A, Z, PRECED IPX & IPY SHOULD BE 
COMMON /SIGMA/ ICHR(4?) SET TO ZERO IN THE 
COMMON / STACKS / ISX(20), IPX, ISY(20), TPY <—TSLOCK DATA SUBPROGRAM 


DATA LEFT, RIGHT, SPACE, PLUS, MINUS,MULT, DIV, DOLLAR, 4,Z 
fi f 48); 66, SH 39, 60; Mf, 42 6%,01,36/ 


. 5 READ (5, 100) CARD AF EXPRESSION STARTS /N COL.1 
100||FORMAT (80AL) 
POINTR = O 
10||\POINTR = POINTR + 1 STOPPING CONVENTION 
IF (C.£Q. DOLLAR) STOP 
IF Se gees (C.LE.2)) G6 76 20 <> 
IF (C.&Q. LEFT GO 7 30 SY BRACKETS 
IF (C.£@. RIGHT) GO 76 40 HY aeackers) 
IF ((C.£@. PLUS) «OR. (C-E0. MINUS) )GO TO 50 
TF ((C.£Q:MULT) .OR. (C.EQ. DIV)) GO TO 50 
IF (C.EQ. SPACE) GO TO 50 


STOP f <KI ERROR: 
c Sek 
CALL PUSHX(C) = 


GO 76 10 
CALL PUSHY (Cc) 

GO 70 10 

CALL POPY (k, OK) 

IF GNOT. OK) STOP 2 
IF (K .£Q. LEFT) GO TO 10 
CALL PUSHX(K) 

GO TO 40 

CALL POPY(K, OK) 

IF (.NOT. OK) GO TO 60 
IF (CK .£Q. LEFT) GO 76 55 
IF (PRECED(C) .GT. PRECED(K)) GO 7 55 

CALL PUSHX(K) PRINT 
GO TO 50 STACK X 
CALL PUSHY (kK) 

CALL PusHy (Cc) 

IF (Cc .NE. SPACE) GO 70 10 













































: 
7 
















896 


FoRivpe+AGDr+~n 1 


Is 









CALL POPX (k, OK) 
F (.NOT. OK) GO 76 5 
WRITE (6, 200) ICHR(K) 
GO TO 1 

200! |FORMAT (1X, At) 










999000000 


Cc 


° 





| END 


EXERCISES) CHAPTER 12 


1%.1 Phe example on simultaneous equations is deficient because it cannot 


cope with zero divisors. Assuming the equations do have a solution this 
deficiency might be overcome by solving the equations in a different 
order. Amend the example so that if there is a zero divisor the 
program looks down the remaining equations to find one with a non- 
zero divisor. Exchange this equation with the one giving trouble and 
proceed as before. 


13.8 ‘Whe example on simultaneous equations deals with just one right-hand 


side but could easily be made fo deal with several = either in 
sequence or simultaneously. In particular, if the right-hand sides 


look like this: 
AS MANY RIGHT-HAND $/DES 
AS THERE ARE EQUATIONS 


then the array of resulting solutions is called the ¢nverse of the matrix 
of coefficients. If you are familiar with matrix algebra write a 
subroutine fo generate an inverse matrix. See Exercise 7.5. 


oo ™ 
—“—O— Oo 
~Oo oO 


12.8 Omprove the shortest roule program so that it makes some checks on the 


validity of input data. For example ensure there is no more than one 
arrow between any two nodes; at least one arrow pointing out of 
the start node; at least one pointing into the finish node. Provide 
error messages to be printed when such checks fail. 


12.44 @hange the program on reverse Polish notation so that if reads real 


12.5 


130 


numbers rather than letters. For example instead of reading: 
A +(8-C)*D-F/(GHH) 
Change stack x fo type &£AL and make the program capable of 
reading: 
6-5 4(3.76 -2)#7.9 — 13-3 /( 3-14 + 2-8) 


tb do this you may find it convenient to modify and use the 
free formaf input routine on page 115. Make if read unsigned 
numbers and ferminate on meeting a space, operator or bracket. 


Laodify the above program so that it has the form of a subroutine: 
SUBROUTINE READER (V, OK) 


where V returns the next item to be read; logical variable Ox reports 
the success or failure of any call to the subroutine. 


Shange the logic dealing with stacks. When an operator is about 
to be pushed on STACK x: 

@ pop two numbers from stTAck x 

@ apply the operator that was 

about to be pushed 

@ push the resulting value on sTAck x 
When one complete expression has been dealt with there should be 
just one value on stack x. This value should be returned by par- 
ameter V of subroutine READER (V, 0K). 


his subroutine permits anyone who writes numerical data for your 
programs ¢ those employing READER) to include expressions instead of 


numbers wherever he or she chooses. 
[ 2*3.14* 27-6 3.9# (1.1-0-03)/ (6% (7-25 +4)) 


BASLIOARAP EY 


Whe Fortran described in this book ¢ Fortran 66) is formally defined in: 


AMERICAN STANDARD FORTRAN 

ANSI X3.9-1966 

American National Standards Institute Inc. 
New York, N.Y. , U.S.A 


A\n amplified description of Fortran 66 including many examples of syntax 
may be found in: 


STANDARD FORTRAN PROGRAMMING MANUAL 
SECOND EDITION 1972 

National Computing Centre Limited 
Manchester, U.K. 

ISBN O 85012 063 2 


A commentary on the above # defining a subset of Fortran 66 for writing 
portable programs * has been made. This commentary is the fruit of many 
years experience in making large programs run on ao range of different com- 
puters in Britain and abroad: 


HECB PROGRAMMING INSTRUCTION MANUAL 

VOLUME II - FORTRAN 1979 

Highway Engineering Computer Branch 
Department of Transport, St Christopher House, 
Southwark Street, London SE1_ , U.K. 


Another study of Fortran for writing portable programs * again drawing upon 
much experience * has been published: 


COMPATIBLE FORTRAN 
A. Colin Day , 1978 
Cambridge University Press 
ISBN O 521 22027 O 


‘Whe full Fortran 77 language is formally defined in: 


AMERICAN NATIONAL STANDARD PROGRAMMING LANGUAGE FORTRAN 

ANSI X3.9- 1978 

American National Standards Institute Inc. 

New York, N.Y. , U.S.A 
and it is well to check whether the facilities offered by the Fortran you are 
using » and which fall outside the scope of Fortran 66 = conform to this 
definition. This particularly concerns the use of files; an area in which 
Fortran 66 is deficient and different Fortrans employ different solutions. 


LAany useful tricks of the programmer’s art are beautifully described in: 


FORTRAN TECHNIQUES (WITH SPECIAL REFERENCE 

TO NON-NUMERICAL APPLICATIONS ) 

A. Colin Day, 1972 

Cambridge University Press 

ISBN O 521 086549 7 & O 521 09713 3 (paperback) 
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WNS4 


A 


adjustable dimensions, 69, 73,77, 97 
arguments 
actual, 56, 58,61,66-7, 69-71 
dummy, 60-1, 66-71, 80,86 
input and output, 69 
arithmetic 
expressions, 26-7, 61 
operators, 26 
arithmetic IF, 42, 45 
array elements, 23, 30 
not allowed, 39,40,43,50, 61,94, 121 
arroys, 48-9 
as arguments, 66-9, 78 
in 1/0 lists, 96 
initialization of, 86-91 
assigned GO TO, 43 
assignments, 30-1, 5,25 
in FUNCTION subprograms, 67 
ASSIGN, 43 
association, 71, 76-8 
asterisk, 26-7, 20 
in DATA statements, 86 
auxiliary 1/0 statements, 


B) 


BACKSPACE, 120-1 
BASIC, vill, 66 
basic externa] functions, 58-9, 62 
using names of, 61, 70 
binar 
digits, 8, 43 
files, 118-21; see also “records” 
blank 
COMMON, 76-7, 87 
record, 100 
blanks (spaces), 20, 12-15, 25, 102, 113 
in data items, 105, 89, 113, 114-15 
BLOCK DATA, 87, 16-17 
uses of, 79, 88 
Boolean values, 28-9, 36, {if 
brackets ( parentheses), 20 
in expressions, 5, 26-9 
in COMPLEX constants, 25 
in FORMAT statements, 99 
in I/o lists, 96-7 


C’ 


120-1 


CALL, 68 

card 
image, 89, 94, 106, 114, 118-19 
punched, 12-13, 6-9 


reader, 5, 9, 94-5 


carriage control, 101, 104 
chains, 82-3, 126-7 
characters 


descriptors (A,H), 2-13, 5 
initialization of, 88 

in run-time formats, 102 
manipulation of, 88-9, 23 
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PRINCIPAL REFERENCES, WHERE SUCH 
EXIST, ARE PUT FIRST IN THE LIST 


standard set of, 
columns 

blank; see “blanks” 

of arrays, 49, 9% 

of punched cards, 
commands, 9, 94 
comment lines (C), 14, éx, 20 
COMMON , 76-9, 16-17 

blank , 76-7 

initialization of, 87, 71 

order of elements in, 76-7 

with EQUIVALENCE, 8! 
compilation, separate, 70 
compiler, Fortran, 8,9 
COMPLEX , 22-3, 17 

arguments, 56, 58 

assignment, 30 

constants, 25 

descriptors (F,E,D,G,P), 106-10 

FORMAT for, 101 

size and precision of, 22, 78 
computed GOTO, 39 
constants, 24 

as arguments, 69 

in DATA statements, 86 
continuation line, 14, 17 
CONTINUE, 39, 41 
control statements, 


D 


data 
deck, 6-9, 13 
files, 7, 118-21 
numbers in, 105 
tape, 6-8 
DATA, 86-7, 
descriptors, 
digits, 20 
binary, 8,43 
in symbolic names, 21 
DIMENSION, 48-9, 17 
DO loop, 40-{ 
implied, 96-7 
DOUBLE PRECISION, 22-3, 17 
arguments, 56, 58, 67 
assignment, 30 
constants, 24 
descriptor (D), 108, 
expressions, 27 
functions, 56,58, 66 
items of data, 105 
size and precision of, 22, 78 


[s 


ENDFILE, 120-1 

END line, 14, 5,66,68,87 

EQUIVALENCE, 80-1, 17 

executable statements, 15, 17 
labels on, 15 

exponent form, 24, 105, 107-9 


20, 14,99 


12-14, 4 
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17, 102 
98-9, 101, 104-13 


104-5 


expressions 
arithmetic, 26-7 
logical, 28-9 
extended range, 41 
EXTERNAL, 70, 17 
fields, 99, 4-5, 13, 104-5 
files, 118-21; see also “records” 
formats 
fixed point, 106, 109 
floating point, 107-9 
free, 88-9, 114-15 
nested, 99 
run-time, 102 
FORMAT, 98-101, 4-5 


label on statement, 15 
formatted (readable) 


files, 118-21 
records, 94-101 
Fortran, viii, 4 
functions, 5, 27 
basic external, 58-9 
intrinsic, 56-7 
statement, 60-1 
FUNCTION , 66-7, 16-17 
GoTo, 38 
graph plot, 103 
Hollerith, Hermann, 12, 5, 25 
Hollerith 
constants, 25, Ti, 102 
items (data), 13, 14, 20, 100 
literals, 113, 25, 102, 104 


Q 


implicit typing, 
implied 
decimal point, 
DO loop, 96-7 
multiplication 
initialization, 86-91 
of characters, 886-9 
1/0 list, 96-7, 94-5 
with descripfors , 
with files, 118-19 
INTEGER, 22-3, 3-4, 17 
arguments, 56,58, 61 
arrays, 48 
assignment, 30 
constants, 24 
descriptor (I), 
expressions, 26 
functions, 56,58, 66-7 
size and precision of , 22-3,15,78,88 
intrinsic functions, 56-7, 62 
using names of, 61, 70 
invocation, 56, 58, 61, 66-7, 68, 71 


J 


23, 48-9, 61, 66 


105-6 


106-13 
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job control language (JcL), 9 


KS 


keywords, 21 
label, 15, 14 


in transfers, 38,39, 42, 43 
of FORMAT statements, 15, 98 
terminating a DO loop, 40-1 
letters (AtoZ), 20-1 
line printer, 101, 8,5, 95, 100 
lines of a program, 14 
little boxes, 3,4, 23, 48, 69 
LOGICAL, 22-3, 17 
arrays, 78 
assignment, 30 
constants, 24 
descriptor (L), 1M 
expressions, 28-9, 36, 61 
functions, 63, 88 
operators , 28-9 
logical IF, 36, 26-9 


GA 


magnetic tape, 7,9, 
matrix, 51,54, 74 
mismatching, 100 
mixed mode, 27 


IN 


name; see 

non-integral power, 

non-standard 
characters, 20 
Fortran, 20, 38, 86 


Co 


octal numbers, 38 
operating systems, 9 


118-20 


“symbolic name” 
31 


operators 

arithmetic, 26 

logical, 29 
order 

of characters, 20, 88-9 

of statements, 17 
parameters 

of DO loops, 40-1! 

of statement functions, 61 
parentheses; see “brackets” 
PAUSE , 38 
peripheral devices, 4-9, 94-5 


portability, wiii-ix, 20,21,23,27,43,78, 114 
precision, 22 
program 
concept of , 3 
deck, 6,8,9 
illustration in Fortran, 
labels in, 45 
main, 16-17 
units, 16-17 


4 
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eek CONTINUED ») 


READ, Qh, & symbolic names, 21 

REAL, 23, 3-4, 17 in basic external functions, 58-9 
arguments, 56, 58, 6! in FUNCTION subprograms, 66-7 
arrays, 48 in intrinsic functions, 56-7 
assignment, 30 in SUBROUTINE subprograms, 87 
constants, 24 in statement functions, 60-1 
descriptors (F,E,D,G,P), 106-7, 109-10 of arrays, 48 
expressions, 26 of BLOCK DATA subprograms, 87 
functions, 56, 58, 61, 66 of comMMoNn blocks, 177 
Size and precision, 22, 78 Ty 

records 
blank, 100-1 tape 
end-of-file (EOF), 120-1 blocks of, 94-5, 118 
physical, 94-5, 118 Poper, 6, 9, 95 

recursion, 16, 96, 98 type, 23, 17 

relational expressions, 28 association of, 76-8 

repeat count character, 88 
in DATA statements, 86 of constants, 24-5 
in FORMAT statements (2), 96, 104 of expressions, 26-7 

rescan, 99, 101, 118-19 of functions, 56-62, 66 

RETURN , 66-7, 686, 71 of variables, 22-3 


reverse Polish notation, 128-9 
REWIND, 120-1 


Roman numerals, 90-1 unconditional transfer, 38 
run-time formot, 102-3 undefined 
Grrays, 49 
input, 100 
scale factor (P), 110, 98, 104 variables, 23, 71, 80 
shapes (structures), 34-5 unformatted (binary) 
shortest route, 126-7 files, 119-214 
sign records, 94-7 
bit, 88 unit number, 4-5, 9, 101 
equals, 20, 5, 30 unsigned integers, 89 
plus and minus, 26, 20, 24 Vi 
simultaneous equations, 46, 124-5 
slash, 26, 20 variables, 22 
in FORMAT statements, 100, 118-19 control, 41 
sorting, 52-3, 54 in 1/0 lists, 96 
spaces; see “bianks” initialization of, 86-91 
stacks, 79, 83, 128-9 symbolic names of, 21 
statement functions, 60-3 vector, 50-1, 67 
state tables, 90-1, 114-15 visual display unit (VDU), 7,9 
STOP, 38, 5 wW 
Storage units, 22,78 
subprograms words, 22, 78, 88 
EXTERNAL, 70 WRITE , 95, 5 
FUNCTION, 66-7 
horrors with, 7! 74 
names of, 214 zero 
SUBROUTINE, 68-9 as a blank in data, 105 
SUBROUTINE, 68-9, 16-17 in column six, 14 
subscripts, 50-2, 80, ~86 in labels, 15 
generation of, 96-7 with arithmetic IF, 42 
ane END 
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[Fortran is @ computer language that has been around for a quarter of 
a century, and is still much used despite predictions throughout its life 
that it would be replaced by more elegant languages. But Fortran is 
still the only language in which it is possible = with care = to write 
truly portable programs. 


Osing his own special formula = original and readable style together with 
unique calligraphy = Donald Alcock now does for Fortran what he did for 
BASIC in his earlier book, Jilustrating BASIC . The Fortran described and 
illustrated is that defined by the American National Standards Institute in 1966 
(Fortran 66) but with allusions to the same Institute’s standard published 
in 1978 @& Fortran 77). 


[For both the person who is new to computers and programming, and the 
person who has met only BASIC before, this book will give an excellent 
introduction to standard Fortran; it also provides a reference manval for 
the language < one that emphasizes the self-discipline needed to write 
portable programs. The book also illustrates a few useful tricks of the 
programmer’s frade. 


“Also by the author 


DSRNA EASTS 


‘ 

Z\\cock takes the reader through the elements of the language, pausing now 
and again to pass on useful tips about list processing, matrix algebra, 
symbol state tables and network analysis ~ indeed, for these tips alone, the 
book makes worthwhile reading for any programmer whether he is familiar 
with Basic or not.’ 


Computer Bulletin 


€ 

Jumane, interesting, comprehensive, and « in paperback form, at least = 
excellent value for money. Congratulations to the author. Very highly 
recommended.’ 
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